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ABSTRACT 


A  computer  code  for  prediction  of  the  marine  atmospner ic 
boundary  layer  is  developed.  The  code  is  used  to  predict 
changes  of  the  capping  inversion  height,  the  strength  of  the 
specific  humidity  and  potential  temperature  'jumps 'at  the 
inversion,  the  well-mixed  relative  humidity  in  the  layer  and 
the  lifting  condensation  level  for  possible  cloud/ fog 
formation.  The  code  is  based  on  recently  formulated 
integrated'  models  for  the  clear  or  cloudy  marine  planetary 
boundary  layer  capped  by  an  inversion.  The  initialization 
is  based  on  radiosonde  data  and,  as  such,  the  code  was 
developed  to  be  used  with  the  Integrated  Refractive  Effects 
Prediction  System  (IREPS)  assessment  code.  It  has  been 
extended  to  be  used  with  any  single  platform  having 
atmospheric  sensing  capabilities,  radiosonde  or  dropsonde. 
IREPS  and  hence  the  MABL  prediction  code  are  under 
consideration  for  inclusion  in  Tactical  Environmental 
Support  System  (TESS). 
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I_. _ INTRODUCTION 

Modern  warfare  has  become  critically  dependent  on  the 
entire  electromagnetic  spectrum  for  command  and  control 
communications,  for  weapons  guidance,  for  electronic  warfare 
support  and  for  countermeasures.  Tactically  essential 
systems  are  highly  affected  by  the  environment  even  when 
conditions  are  not  severe  in  tne  historical  sense. 
Enhancement  or  degradation  of  the  operational  performances 
of  EM/EO  systems  has  become  a  primary  concern  of  tasK  force 
commanders.  As  such,  the  deployment  of  resources  and  the 
modification  of  tactics  based  upon  environmental  factors  in 
EM/EO  propagation  will,  to  a  very  large  extent,  determine 
the  effectiveness  of  sensor,  weapon  and  communications 
systems. 

Environmental  effects  on  EO/EM  system  performance  can  be 
grouped  in  the  general  categories  of  refraction  (EM),  wave 
front  distortion  (EO) ,  extinction  (EO) ,  and  dispersion 
(EM/EO) .  Environmental  factors  contributing  to  these 
effects  are  the  vertical  gradients  of  pressure,  temperature 
and  humidity  for  refraction,  small  scale  inhomogenities 
(turbulence)  of  the  index  of  refraction  for  wave  front 
distortion,  concentration  of  water  vapor  and  size 
distributions  of  aerosols  for  extinction  and  turbulent 
transport,  for  dispersion.  All  of  these  factors  are 
multi- var iable  in  terms  of  their  dependence  on  the  routinely 


measured  and  predicted  meteorological  variables:  pressure, 
wind,  temperature  and  moisture.  The  gradient  of  the  index 
of  refraction  is  the  only  factor  for  wnich  the  required 
accuracy  exists  for  direct  description  by  existing  or 
forseen  measurement  capabilities.  Turbulence  and  aerosol 
descriptions  will  have  to  be  obtained  by  indirect  metnods 
due  to  measurement  complexities  which  preclude  measurement 
on  operational  snips. 

There  are  requirements  to  describe  all  of  these  factors 
for  the  atmospheric  region  extending  from  the  surface  to 
0.5-1. 5  km  above  the  surface,  spanning  the  Marine 
Atmospheric  Boundary  Layer  (MABL) .  The  MABL  is  cooler  and 
more  moist  than  the  overlying  air  and  is  capped  by  a  layer 
(inversion)  ,  50-100  meters  thick,  in  whicn  temperature 
increases  and  humidity  decreases  with  height.  Critical 
values  of  the  index  of  refraction  gradients  leading  to 
anamalous  EM  propagation  can  be  present  within  the  shallow 
capping  inversion  because  of  the  humidity  and  temperature 
gradients.  The  entire  affected  region  (duct)  extends  below 
the  inversion,  and  determining  its  lower  boundary  is 
essential.  High  turbulence  can  be  present  in  the  inversion 
as  well  as  in  the  surface  layer  and  this  can  cause 
degradation  of  slant  path  optical  propagation.  Turbulence 
intensities  in  the  inversion  are  usually  one  to  two  orders 
of  magnitude  greater  than  they  are  immediately  below, and  ar 
two  to  three  orders  of  magnitude  smaller  immediately  above 


the  inversion.  Tactically  significant  degradation  of 
infrared  energy  is  due  to  water  vapor  absorption  and  to 
scattering  by  marine  aerosols,  and  is  thus  usually 
restricted  to  moist  MABL. 

All  of  the  above  features  depend  critically  on  local 
detailed  vertical  structure  of  the  atmosphere.  Therefore 
regional  climatologies  are  not  useful  for  operational 
predictions,  because  the  significant  vertical  structures  are 
lost  in  averaging,  and  also  because  most  nistorical 
measurements  do  not  relate  to  specific  EM/EO  requirements. 
Furthermore,  the  gradients,  and  hence,  the  affected  regions 
cannot  be  explained  on  the  basis  of  large  scale  atmospheric 
flows  since  the  nature  of  the  gradients  are  due  to  dynamic 
processes  which  are  controlled  by  near  surface  as  well  as  by 
larger  scale  features. 

We  believe  a  'gap'  has  existed  in  past  efforts  to 
characterize  the  tactical  environmental  conditions.  The  gap 
is  between  two  extreme  approaches  of  relating  conditions  to 
1)  near  surface  observations,  and  2)  to  larger  scale 
predicted  synoptic  patterns.  Clearly,  to  assess  the  above 
features,  local  measurements  are  desirable.  Measurements 
(radiosondes)  are  made  infrequently  and,  as  time  progresses, 
the  initial  point  measurement  becomes  less  applicable  and  a 
predictive  scneme  is  needed.  One  must  consider  a  transition 
to  climatology,  large  scale  numerical  analyses  and 
predictions,  or  dynamic  models  based  on  the  initial 
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soundings  which  are  available  at  the  operational  location. 
Climatology  has  already  been  addressed  and  discounted. 
Ruggles  (1975)  has  examined  the  capabilities  of  large  scale 
numerical  procedures  at  remote  facilities  and  has  argued 
convincingly  that  they  are  not  sufficient. 

Needed  are  characterizations  of  coupled  local  oceanic 
and  atmospheric  mixed  layer  features  for  time  scales  from  12 
to  18  hours  and  for  spatial  scales  from  50  to  300  km. 
Significant  cnanges  occur  in  both  mixed  layers  over  these 
temporal  and  spatial  scales.  Predicting  these  changes  is 
now  possible  using  available  measurements  and  reasonable 
physical  models  which  have  recently  been  formulated. 

This  paper  will  describe  results  from  observational  and 
model  evaluation/formulation  efforts  with  these  models  which 
provide  evidence  that  tactically  relevant  forecasts  can  be 
made  over  the  required  scales.  The  MABL  is  emphasized  in 
this  discussion  but,  as  will  be  seen,  the  MABL  model  must 
eventually  be  coupled  with  an-  Ocean  Boundary  Layer  (OBL) 
model  to  adequately  describe  the  system.  This  is  the  object 
of  separate  investigation  by  the  author  and  faculty  at  NPS. 

The  long  term  objective  of  work  in  these  areas  is  to 
make  reliable  5-18  hour  forecasts  of  properties  within  the 
MABL  and  OBL.  The  near  term  objectives  are  1)  to  verify  and 
reformulate  existing  models  for  the  responsible  physical 
processes  and  2)  to  develop  measurement  procedures  which  are 
suitable  for  routine  observations  from  a  single  ship  or 
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aircraft.  For  the  finalized  model  it  will  be  necessary  to 
predict  such  properties  as  the  height  of  the  inversion,  the 
strength  of  the  gradients  and  fluxes  at  the  inversion  and  at 
the  surface,  and  the  values  of  well  mixed  parameters  in  the 
region  between  the  surface  and  the  inversion. 
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II.  BACKGROUND 


This  section  will  provide  a  general  outline  of  the 
environmental  descriptions  required  and  the  models  which  can 
used  to  predict  them.  The  possible  input  parameters  and 
models  are  delineated  in  Figure  1.  Tactical  forecasts  to 
which  the  improved  understandings  would  apply  appear  on  the 
extreme  right  hand  side.  Single  station  assessment 
constraints  associated  with  the  local  atmospheric  and  ocean 
models  are  timely  and  compatible  with  the  concept  of  the 
Tactical  Environment  Support  System  (TESS)  as  recently 
described  by  NEPRF  scientific  personnel  (NEPRF  METRO  Report, 
July  1981) . 

As  stated,  we  will  consider  only  the  local  atmospheric 
mixed  layer  model  and  associated  descriptions.  The  local 
oceanic  predictions  would  be  made  by  an  oceanic  mixed-layer 
model  (e.g.  Garwood,  1977)  and  the  regional  predictions 
would  be  made  by  a  three  dimensional  atmospheric  (e.g. 
NORAPS)  and  a  regional  oceanic  (e.g.  Tactical  Ocean 
Prediction  System-TOPS)  numerical  prediction  model.  The 
ocean  mixed  layer  model  and  the  regional  numerical 
prediction  schemes  are  objectives  of  other  basic  and 
exploratory  Navy  sponsored  research  efforts. 

A.  METEOROLOGICAL  DESCRIPTIONS  AND  MODELS 

From  a  local  perspective,  let  us  consider  an 
idealization  of  the  oceanic-atmospheric  system.  The  sea-air 


interface  is  bordered  by  oceanic  and  atmospheric  turbulent 
mixed  layers  which  effectively  insulate  the  bulk  ocean  and 
atmospheric  regions.  The  primary  sources  of  the  turbulence 
within  the  layers  are  the  velocity  (current)  and  buoyancy 
(density)  gradients  at  the  interface.  Even  under  conditions 
where  the  water  is  slightly  cooler  than  the  air,  buoyancy 
forced  velocity  fluctuations  within  the  layer  can  be  quite 
large  and  mix  the  entire  MABL  from  the  surface  to  the 
inversion.  The  large  vertical  mixing  yields  constant 
(well-mixed)  wind,  temperature  and  humidity  profiles  above 
the  surface.  At  the  top  of  the  atmospheric  mixed  layer 
there  is -the  inversion,  a  thin  transition  region.  Examples 
of  observed  well-mixed  profiles  and  the  inversion  are  shown 
in  Figure  2(a). 

The  atmospheric  mixed  layer  interacts  with  the  free 
atmosphere  at  this  interface  by  means  of  turbulence  forced 
entrainment.  Entrainment  brings  dry,  warm  air  into  the 
mixed  layer  and  also  increases  the  surface  winds  if  there 
are  higher  winds  aloft.  Stratus  clouds  form  within  the 
layer  if  entrainment  causes  the  moist  mixed  layer  to  extend 
above  the  lifting  condensation  level. 
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Fig.  2.  Simplified  meteorological  and  tactical  feature*  of  the  *BL  and 
overlying  region  including  (a)  humidity  mixing  ratio,  Q/o,  and 
potential  tmopsrature,  a,  profiles  and  subregions,  (b)  modified 
index  of  refraction,  N,  profile,  (cl  optical  turbulence  inten¬ 
sity,  C^2,  profiles,  (d)  total  extinction  coefficient,  > 
(absorption  and  scattering) ,  profile. 


Over  a  time  scale  of  24-36  hours,  physical  processes 
within  the  mixed  layer  are  also  controlled  by  the  large 
scale,  non- stationary,  synoptic  flows.  Synoptic  scale  mean 
vertical  motion  and  advection  will  have  to  be  prescribed  or 
predicted  over  the  forecast  period. 

The  previously  inferred  well-mixed  nature  of  the 
convective  MABL  has  implications  for  the  vertical 
distributions  of  mean  values  and  the  vertical  flux  scales  of 
wind-(J*,  tempera ture-T*  and  humidity-Q*.  One  implication 
already  illustrated  is  that  properties  wnich  are  conserved 
during  mixing  can  be  treated  as  being  constant  with  height 
within  the  MABL.  These  parameters  are  the  specific  humidity 
and  potential  temperature  for  a  clear  MABL  and  equivalent 
potential  temperature  and  total  (vapor  plus  liquid)  specific 
humidity  for  a  cloudy  MABL.  A  second  implication  is  that 
vertical  fluxes  of  the  well-mixed  parameters  vary  linearly 
with  height. 

These  implications  enable  predictions  of  MABL  evolution 
to  be  based  solely  on  fluxes  at  the  upper  (inversion)  and 
lower  (surface)  boundaries  and  the  large  scale  subsidence 
and  advection.  They  form  the  basis  of  recent  model 
formulations  by  Deardorff  (1976)  and  Stage  and  Businger 
(1981).  Fluxes  at  both  boundaries  are  due  to  buoyant  and 
mechanically  generated  turbulence  for  the  clear  case.  The 
linear  height  variations  of  the  fluxes  allows  buoyant  fluxes 
at  the  inversion  to  be  more  readily  estimated  on  the  basis 


-20- 


of  general  cloud  features.  Methods  exist  for  estimating 
synoptic  scale  forcing  from  single  station  measurements  but 
further  efforts  are  required  to  achieve  the  accuracies 
required  in  MABL  predictions. 

In  general,  existing  models  are  quite  good  for  the  clear 
sky  MABL  and  fair  to  good  for  the  cloudy  MABL.  Considerable 
effort  is  now  being  directed  to  improving  the  models  for  the 
cloudy  MABL.  The  existence  of  cloud  layers  in  the  MABL  is 
very  sensitive  to  both  the  sea  surface  temperature  and  the 
entrainment  rate  and,  in  turn,  clouds  have  a  profound  effect 
on  the  short  and  long  wave  radiation  budget  at  the  surface. 
Fairall  et  al  (1981)  provide  a  discussion  on  the  approach 
and  the  status  of  abilities  to  predict  the  evolution  of  the 
MABL. 

B.  TACTICAL  DESCRIPTIONS 

Figures  2b-d  illustrate  tactical  descriptions  discussed 
in  the  Introduction  (appearing  in  the  right  hand  column  of 
Figurel) .  Meteorological  processes  and  features  relevant  to 
these  tactical  descriptions  are  those  that  were  discussed 
above  (Figure  2a  and  column  4  of  Figure  1).  Several 
meteorological  factors  affect  each  tactical  description  as 
indicated  by  the  numbers,  1-4,  in  Figure  1.  A  few  of  the 
significant  properties  and  values  in  the  tactical 
descriptions  are  described  in  the  following  paragrapns. 
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The  M  profile  in  Figure  2b  describes  tne  refracted  EM 
ray  radius  of  curvature,  relative  to  the  earth's  radius  of 
curvature.  If  M  decreases  with  height,  a  trapping  layer  is 
formed  because  the  EM  ray  will  bend  downward  relative  to  the 
earth.  This  causes  the  formation  of  a  duct.  The  upper 
boundary  of  the  duct  is  at  the  M  minimum;  the  lower  boundary 
is  at  the  height  where  this  same  value  of  :i  occurs  below  the 
trapping  layer,  as  shown  in  Figure  2b.  The  M  profile  is 
determined  by  pressure,  temperature  and  humidity  profiles, 
so  values  of  well  mixed  temperature  and  numidity  determine 
the  lower  boundary  of  the  duct.  A  decrease  of  mixed  layer  M 
value  with  time  could  cause  a  surface  based  duct  to  become 
an  elavated  duct  as  shown  in  Figure  2b.  This  could  occur 
if  the  mixed  layer  became  warmer  and  drier  due  to 
entrainment  of  overlying  air.  Warming  and  drying  of  the 
mixed  layer  by  entrainment  would  also  increase  the  height  of 
the  evaporation  duct,  a  ducting  layer  immediately  adjacent 
to  the  surface.  Examples  of  extended  ranges  and  holes 
occurring  with  a  duct  appear  in  Figure  3  (Hitney,  1979). 

The  Cn2  profile  in  Figure  2c  shows  the  vertical 
variation  of  turbulence  which  would  affect  optical  wave 
front  distortion.  We  see  that  Cn2  is  largest  near  the 
surface  and  in  the  inversion.  Values  in  these  regions  are 
1-2  orders  of  magnitude  larger  than  those  in  the  mixed  layer 
and  above  the  inversion.  Cn2  values  near  10 (-14)  m  -2/3  are 
representative  for  these  two  regions.  Tne  importance  of 
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this  value  is  illustrated  in  Figure  4  with  a  simulated  image 
of  a  remotely  piloted  vehicle  as  viewed  througn  Cn2  regions 
of  0  (no  turbulence),  10{-15),  and  10 (— 14 )  m  -2/3  (Kearns 
and  Walter,  1978) . 

The  extinction  coefficient,  B,  profile  is  shown  in 
Figure  2d.  B  is  inversely  proportional  to  the  range  of  IR 
systems  and  primarily  depends  on  the  absolute  humidity  and 
the  aerosol  size  distribution.  Aerosol  size  distributions 
depend  on  the  generation  rate  of  sea  salt  particles  and  on 
the  relative  humidity  as  determined  by  both  the  moisture  and 
temperature.  Therefore,  IR  extinction  is  very  dependent  on 
the  entrainment  of  relatively  clean,  dry  and  warm  air  from 
above  the  inversion  and  on  the  surface  wind  generation  of 
the  sea  salt  particles  (Davidson  et  al,  1982).  As 
illustrated,  the  extinction  coefficient  is  large  within  tne 
mixed  layer  and  increases  with  relative  humidity  which 
increases  with  height.  The  increase  of  extinction  with 
height  is  important  in  3lant  path  range  considerations. 


Fig.  3.  Ex  an  pies  of  ascended  rangna  and  holes  for  EM  discs  for  (a)  surface- 
basnd  duct  and  a  shipboard  air-anarch  radar  and  (b)  elevated  duct  end 
an  airbomn  sarly-waminq  radar  (froa  Hitnny,  1979). 


Fig.  4.  Canpucar  a  inula  tad  iaaga  raiocaly  piloted  vehicle  iRtV)  as  viewed  with 
3-12  a  Si  sensor  (range  2  ka>  visibility  IS  <a,  threshold  4.  5  rU) 
through  ionosphere  with  C^2  values  of  (a)  >  (no  turbulence),  (b)  3.7  x 
ia'iS  a_2/3  and  (c)  1.  0  x  10“i4  a~2/^  (from  Kearns  and  Walter,  1978). 
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III.  MABL  OBSERVATIONS  AND  PREDICTIONS 
A  recently  obtained  data  set  will  be  used  to  demonstrate 
the  prediction  of  changes  in  the  MABL.  The  demonstration 
will  show: 

a)  the  status  of  available  data  sets.- 

b)  that  the  assessment  of  relative  roles  of  dynamic 
processes  in  the  MABL  can  be  based  on  a  relatively  simple 
physical  model  which  includes  specifiable  inputs, 

c)  that,  at  the  present  time,  we  only  have  limited 
evidence  that  changes  can  be  predicted.  The  rigorous 
specification  or  description  will  require  further 
interpretive  efforts  and,  perhaps,  improved  formulation 
of  the  separate  models. 

A.  MABL  OBSERVATIONS 

The  data  were  obtained  during  the  Cooperative  Experiment 
on  West  Coast  Oceanography  and  Meteorology  (CEWCOM-78) 
conducted  west  of  San  Nicolas  Island,  CA  during  May  of  1978. 
Observations  of  the  oceanic  and  atmospheric  mixed  layers 
were  made  from  the  R/V  ACANIA;  radiosonde  observations  were 
also  made  at  surrounding  shore  stations.  The  data  to  be 
shown  are  from  a  48-hour  period,  5/19/1200  to  5/21/1200  PST. 
The  R/V  ACANIA  was  cruising  slowly  (2-3  knots)  into  the  wind 
and  returned  to  an  initial  point  approximately  every  12 
hours.  The  general  location  of  the  R/v  ACANIA  during  the 


5/19  to  5/22  periods  and  locations  of  surrounding  shoreline 


radiosonde  sites  appear  in  Figure  5. 

The  period  was  one  of  steady  onshore  flow  caused  by  the 
combined  effects  of  intensification  of  the  Eastern  Pacific 
High  and  the  persistence  of  the  Mexican  thermal  low.  The 
only  apparent  change  in  synoptic  scale  forcing  was  an 
increase  in  the  offshore  pressure  gradient.  Advection  in 
the  atmosphere  was  moderate.  However,  MABL  evolutions  were 
primarily  determined  by  subsidence,  surface  fluxes,  and 
entrainment  at  the  inversion.  Satellite  imagery  showed 
uniformly  increasing  stratus  coverage  (thin  to  heavy)  during 
the  period  with  a  cellular  (broken)  coverage  occurring  late 
on  the  21st. 

Acoustic  sounder  returns,  mean  surface  layer  parameters, 
and  sea  surface  temperature  during  the  19-21  May  period  were 
all  measured  from  the  R/V  ACANIA  and  are  shown  in  Figures 
6b-d.  Composite  potential  temperature  profiles  from 
shipboard  and  shore  station  radiosonde  and  temperature 
profiles  from  ship  deployed  XBT's  appear  in  Figures  6a  and 
6e,  respectively.  Although  these  descriptions  of  MABL  and 
OBL  changes  were  obtained  from  an  instrumented  research 
vessel,  all  can  be  obtained  from  operational  ships.  This 
includes  the  acoustic  sounder  record. 
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Tract  of  R/V  Acania  and  location  of  shoreline  radiosonde  sites  during 
CEWCCM-78.  General  location  of  R/V  Acania  during  19-21  May  period  is 
indicated  by  hatched  area  N-NW  of  San  Nicolas  Island  (SNZ). 


Changes  of  atmospheric  features  which  would  have  been 
tactically  significant  in  view  of  the  time  scale  and  their 
magnitudes  were: 

a)  The  MABL  depth  increased  from  250  m  to  750  m  over 
the  48  hour  period.  The  changes  occurred  over 
relatively  short  intervals,  from  20/00  to  20/04, 
and  from  21/00  to  21/08.  The  MA3L  depth  remained 
nearly  constant  during  the  intervening  12  to  18 
hour  periods. 

b)  The  surface  layer  temperature  increase  from  20/00 
to  20/15  is  indicative  of  entrainment  of  overlying 
warm  air.  The  entrainment  was  presumably  also  a 
factor  in  the  wind  speed  increase  over  the  same 
period,  which  caused  the  warm  shallow  ocean  layer 
to  be  deepened  by  mixing,  lowering  the  sea  surface 
temperature  after  20/1800  (Figure  6d) . 

In  conjunction  with  changes  in  MABL  structures  and 
parameter  values  depicted  in  Figure  6,  changes  occurred 
which  are  important  to  operational  systems.  Some  of  these 
changes  can  be  determined  directly  from  parameters  measured 
in  this  scientific  observational  effort.  M  profile  and 
cloud  base  heights  and  condensation  levels  can  be  calculated 
from  the  measured  composite  humidity  and  temperature 
profiles.  Surface  layer  extinction  coefficients  have  to  be 
determined  from  the  measured  humidity  and  aerosol  size 
distributions.  The  method  for  determining  extinction 
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coefficients  from  aerosol  data  has  been  described  by 
Schacher  et  al  (1981) .  Surface  layer  Cn2  values  and 
evaporation  duct  heignt  can  be  determined  with  considerable 
certainty  from  measured  surface  layer  wind,  temperature, 
humidity  and  surface  temperature  values  using  expressions 
described  by  Davidson  et  al  (1981)  and  Fairall  et  al  (1978). 
Extinction  and  Cn2  profiles  have  to  be  based  on  more  recent 
and,  hence,  less  substantiated  expressions  such  as  those  by 
Wells  et  al  (1977)  and  Wyngaard  and  Lemone  (1980) . 

Results  from  these  determinations  appear  in  Figures  7 
and  8  along  with  the  composite  potential  temperature,  0,  and 
water  vapor  mixing  ratio,  q,  profiles  (Figure  7a)  and  the 
inversion  height  (Figure  8a) .  The  observed  results  indicate 
that  for  the  observation  period: 

a)  The  EM  duct  associated  with  the  inversion  gradients 
evolved  from  being  surface  based  to  being  elevated 
(Figure  7b)  and  the  evaporation  duct  height  ranged 
from  below  5  m  to  above  8  m  (Figure  8b), 

b)  Optical  turbulence,  Cn2 ,  in  the  surface  layer 
varied  from  10 (-15)  to  10  (-14)  m  -2/3  over  the 
period  (Figure  3c)  and  values  in  the  inversion 
increased  from  10 (-16)  to  10 (—15)  m  -2/3  (Figure 
8c)  , 

c)  Surface  layer  aerosol  extinction  in  the  8-12  urn  IR 
region  increased  during  the  first  6  hours  (5/19 
1200-1800)  because  the  relative  numidity  increased. 
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It  decreased  during  the  last  18  hours  (5/20/1800  to 
5/21/1200)  because  of  relative  humidity  and  wind 
decreases  (Figure  3d) , 

d)  Stratus  clouds  persisted  through  the  period.  The 
base  lifted  from  about  100  m  at  5/19/1800  to  600  m 
at  5/21/1200.  This  is  evident  in  the  extinction 
profile  (Figure  7d) . 
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Fig  .  7 

Observed/computed  and  predicted  ^A8L  profiles  for  5/19/1200  to 

5/21/1200  CEWCCM-78  period  including  (a)  Q  (g  m”^)  and 
2  -2/3 

(b)  M,  (c)  CM  (ta  )  and  (d)  total  extinction  coefficient,  3 
-1  ^ 

).  Solid  lines  ccr-espond  to  :bse”ved  and  dashed  to  model 
prediction. 
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Fig .  8 

Observed/ computed  and  predicted  surface  layer  values  for 
5/19/1200  to  5/21/1200  CSWCOI-78  period  including  (a)  inversion 
height  and  lifting  condensation  level,  (b)  evaporation  duct, 

(c|  Cj,2  m'2/3  and  (d)  aerosol  scattering  extinction  coefficient 
3.  X's  are  observed  and  dashed  lines  are  model  predictions. 


B.  MOOEL  PREDICTIONS 

During  the  past  £ive  years,  investigations  were 
performed  on  surface  layer  and  mixed-layer  scaling  of  small 
scale  turbulence  and  aerosol  properties  (Davidson  et  al, 
1978;  Fairall  et  al,  1980;  Scnacher  et  al,  1981a  and  1981b). 
The  investigations  led  to  the  application  of  a  time 
dependent  MABL  models  based  on  entrainment  energetics  and  on 
cloud  radiative  transfers  suggested  by  both  Deardorff  (1976) 
and  Stage  and  Businger  (1981).  The  model  requires  as  input 
an  initial  atmospheric  sounding  (radiosonde) ,  the  mean  winds 
at  a  level  within  the  surface  layer  (10-30  meters)  and  the 
surface  temperature.  Well-mixed  temperature  and  humidity  are 
predicted  so  the  surface  wind  and  wind  shear  at  an  inversion 
are  the  only  local  atmospheric  variables  which  have  to  be 
prescribed.  The  larger  scale  subsidence  and  advection  must 
also  be  prescribed  for  the  forecast  period. 

The  steps  in  the  prediction  computation  are  shown  in 
Figure  9,  where  it  is  noted  that  procedures  are  the  same  for 
clear  and  cloudy  cases  except  for  entrainment  computation 
and  estimating  cloud  top  cooling.  Fairall  et  al  (1981) 
include  comprehensive  discussions  of  how  applicable  physical 
processes  are  teated  in  the  model  and  how  available  data  can 
be  used  to  estimate  the  processes.  Because  of  the 
simplified  physical  model,  the  computations  do  not  require 
numerical  integration  on  a  vertical  grid  so  the  computer 
storage  requirements  can  be  satisfied  by  available  shipboard 
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microcomputers.  In  fact,  the  use  of  a  hand-held  calculator 
is  being  considered  for  these  procedures.  Particular 
computations  of  interest  are  those  for  the  evaporation  duct, 
the  surface  flux,  flux  of  total  heat  and  momentum  and  the 
depth  change  of  the  mixed  layer.  Routines  for  computing  and 
graphing  the  M-profile  can  be  system  are  used  on  the 
HP-41C/V  and  printer  hand-held  system. 

Predictions  of  MABL  changes  were  made  with  the 
simplified  model  for  the  48  hour  period  (5/19/1200  to 
5/21/1200)  corresponding  to  observed  results  in  Figure  6a-e. 
The  predictions  were  made  for  two  separate  24  hour  periods 
starting  at  19/1200  and  at  20/1200.  The  initial  profile  for 
the  19/1200  starting  time  was  that  obtained  at  San  Clemente 
Island  (SCI  in  Figure  5)  at  approximately  1000.  Initial 
profiles  for  the  20/1200  starting  time  were  averages  of  the 
composites  at  20/0500  and  20/1700.  The  sea  surface 
temperatures  were  those  at  each  start  time  as  shown  in 
Figure  6d  and  were  15  and  14.5  C,  respectively.  The  surface 
wind  observed  during  the  forecast  period  was  specified  to 
increase  linearly  from  5  m/sec  at  19/1200  to  10  m/sec  at 
20/2100  and  then  to  decrease  linearly  to  7  m/sec  at  21/1200. 

Larger  scale  synoptic  advection,  was  estimated  to  be 
zero  from  the  thermal  wind.  Subsidence  for  the  second 
period,  starting  at  20/1200,  was  adjusted  on  the  basis  of 
that  value  required  for  agreement  between  predicted  and 
observed  mixed  layer  depths  during  the  first  period.  Hence, 
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the  model  was  also  used  to  estimate  the  "most  probable” 
synoptic  scale  subsidence.  The  tuning  of  the  subsidence 
value,  based  on  comparisons  of  the  latest  sounding  and 
predictions  from  a  previous  sounding,  has  been  described  as 
"hindcasting" . 
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The  predictions  are  illustrated  in  Figures  7  and  8  along 
with  the  observed  values.  These  predictions  were  made  using 
the  Deardorff  (1976)  cloud  topped  entrainment 
parameterization  and  Davidson  et  al  (1978)  surface  layer 
parameterization.  Results  for  the  same  period  using  the 
Stage  and  Businger  entrainment  parameterization  appear  in 
Appendix  B.  The  agreement  near  20/1200  in  Figure  8  occurs 
because  the  model  was  reinitialized  at  that  point. 
Significant  outcomes  from  comparisons  of  observed  and 
predicted  results  are: 

a)  The  mixed  layer  depth  (inversion  height,  Figure  8a) 
was  predicted  quite  well.  This  increase  during  the 
later  part  of  both  periods  was  caused  by  cloud  top 
entrainment  because  the  predicted  surface  layer  became 
stable  and,  hence,  there  existed  no  surface  buoyancy 
flux  forcing. 

b)  The  cloud  cover  was  accurately  predicted  to  persist 
throughout  the  first  period  and  for  the  first  15  hours 
of  the  second  period,  as  indicated  by  the  condensation 
level  being  below  the  inversion  in  Figure  8a.  However, 
the  predicted  cloud  thicxness  was  less  than  that  which 
actually  occurred.  As  indicated  in  preceding 
discussions,  the  cloud  cover  was  observed  to  breaX  up  on 
the  21st  and  the  prediction  after  21/0600  agres  witn 


that , 


c)  The  well  mixed  temperature  and  humidity  and  tne  jump 
predictions  (Figure  7a)  are  reasonable  except  at 
20/0500.  The  predicted  well  mixed  temperature  was  too 
warm  by  1-2  C  at  all  other  comparison  times.  This  led  to 
the  surface  layer  becoming  stable  when  in  actuality  it 
remained  unstable  throughout  both  periods  (i.e.  surface 
temperature  higher  than  10  m  temperature . )  It  is 
believe  that  the  cooling  associated  with  the  clouds  was 
not  sufficient.  It  will  be  seen  that  tne  surface  layer 
was  too  warm  and  stable  during  the  night  hours. 

d)  The  predicted  M  profiles  (Figure  7b)  were  very 
accurate  in  terms  of  the  surface  based  duct  evolving  to 
an  elevated  duct.  It  is  interesting  that  even  wnen  the 
predicted  well  mixed  temperature  and  humidity  were  much 
too  high  (20/0500)  the  M  profile  was  verified.  This 
occurred  because  their  respective  affects  on  M  were 
compensated . 

e)  The  predicted  Cn2  profiles  (Figure  7c)  gave  good 
values  for  tne  inversion  region  and  reasonable  values 
for  the  mixed  layer,  except  for  20/0500.  Agreement 
within  a  factor  of  three  is  considered  reasonable  for 
Cn2.  The  inversion  region  values  were  cased  on  the 
convective  mixing  velocity,  W*  (Wyngaard  and  Lemone, 

1980)  for  unstable  conditions.  Tne  friction  velocity  , 
U*,  was  used  instead  of  W* ,  wnen  a  stable  surface  layer 
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existed  (20/0500).  Values  above  the  inversion  were  set 
to  10 (—13)  for  both  ooserved  and  predicted  profiles. 

f)  The  total  predicted  extinction  at  3-12  urn  (Figure  7d) 
was  quite  good  outside  of  cloud  layers.  Again,  values 
above  the  inversion  were  set  to  the  same  value  for  both 
observed  and  predicted  profiles.  The  sensitivity  of 
aerosol  size  to  relative  humidity  Decomes  larger  wnen 
the  relative  humidity  is  above  90%.  Hence  we  expect 
this  comparison  to  be  poor  in  tne  cloud  region  wnere  tne 
relative  humidity  approaches  100%.  However,  the  good 
agreement  outside  of  clouds  is  encouraging. 

g)  The  evaporation  duct  height  (Z*)  prediction  (Figure 
3b)  is  not  very  good.  The  disagreement  occurs  because 
the  predicted  surface  layer  was  more  stable  than  was 
observed.  Z*  increases  almost  linearly  with  increasing 
stability.  Very  large  Z*  values  were  predicted  at  the 
end  of  the  second  period  because  it  was  stable  and 
because  the.  predicted  well  mixed  humidity  was  too  low. 
Figure  7a. 

n)  The  predicted  surface  layer  Cn2  values  (Figure  8c.) 
were  always  less  than  those  observed  Decause  of  the 
previously  mentioned  inaccuracy  in  surface  layer 
stability.  Minima  at  19/2300,  20/1300  and  20/0500 
correspond  to  the  surface  layer  oeing  neutral  as  it 
passed  from  unstable  to  staDle,  unstable  to  stable  and 
stable  to  unstaDle  conditions,  respectively. 
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i)  The  predicted  surface  layer  aerosol  scattering 
extinction  for  8-12  um  (Figure  3d)  was  very  good  during 
tne  second  period.  The  observed  and  predicted  decrease 
from  20/1800  to  21/1200  was  associated  with  a  decrease 
in  relative  humidity  and  wind  speed.  The  wind  speed 
used  for  the  prediction  was,  of  course,  prescribed  on 
the  basis  of  observed  results.  Hence,  the  role  of 
relative  numidity  on  the  extinction  is  tnat  oeing 
evaluated. 

C.  PREDICTION  INTERPRETATION  AND  APPLICATION 

The  demonstrated  capability  of  this  model,  which  is 
based  on  routinely  observed  input  data,  should  be  of 
assistance  to  the  geophysics  officer  in  estimating  possible 
changes  in  the  MASL  and  in  the  inversion  and  the  potential 
impact  upon  tactically  important  parameters.  Exploitation  of 
the  existing  ducts  is  presently  possible  with  IREPS. 
Specifically,  the  duct  changes  between  elevated  and  surface 
based  nature  are  of  considerable  concern  to  placement  of 
platforms  to  ennance  or  to  minimize  communication  linxs. 
Alerting  commanders  and  airborne  assets  in  advance  to 
possible  changes  is  a  capability  not  previously  availaole  in 
on  board  assessment  systems.  Of  additional  use  is  the 
computation  of  the  lifting  condensation  level  to  indicate 
the  occurence  of  clouds. 
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W.  CONCLUSIONS 


The  importance  or  and  an  approacn  for  predicting 
evolutions  of  the  marine  atmospneric  boundary  layer  nave 
been  described.  Cnanges  in  the  mixed  layer  deptn  strengtns 
of  overlying  inversion  and  well-mixed  properties  lead  to 
significant  cnanges  in  the  weapons/systems  tactical 
environment.  These  include  £.4  duct  regions,  optical 
turbulence  ana  optical  extinction,  DOtn  within  and  at  tne 
ooundaries  of  the  inversion  capped  mixed  layer. 

Altnougn  the  tactical  descriptions  are  multi-var iable 
and  require  detailed  vertical  descriptions,  a  simplified 
model  based  on  routinely  observed  data  appears  to  oe  quite 
useful  for  predicting  cnanges  over  12-18  hours.  For  a  cioua 
topped  period,  it  was  demonstrated  that  tne  refraction, 
optical  turbulence  and  optical  extinction  profiles  were 
predicted  quite  well,  as  were  cloud  coverage  cnanges. 

The  primary  differences  between  observation  and 
predictions  were  for  tne  near  surface  layer  and  pertain  to 
tne  evaporation  duct  height  and  the  optical  turbulence.  The 
predicted  evaporation  duct  was  mucn  too  nign  and  tne 
predicted  optical  turbulence  was  much  too  low.  In  botn 
cases,  tne  predicted  air  temperature  was  too  nign  (2-1  C) 
wnicn  resulted  in  a  neutral  to  stable  condition.  Tne  cause 
for  this  appears  to  oe  due  to  improper  specification  of 


cloud  cooling  witnin  cne  mixed  layer.  This  is  an  area  in 
wnicn  researcn  is  still  oeing  performed. 

Two  scnemes  are  presently  being  evaluated  for 
specification  of  cloud  top  cooling  and  entrainment.  Stage 
and  Businger  (1981)  and  Dearcorff  (1976).  Both  address 
entrainment  and  cloud  top  cooling  for  snort  and  long  wave 
radiation,  out  large  disparities  exist  in  actual 
parameterization.  More  study  in  the  variance  of  absorption 
witn  insolation  and  cloud  buoyancy  should  resolve  tne 
difference.  Both  of  these  effects  may  nave  great  impact  on 
tne  predicted  results  oecause  of  the  bulk  nature  of  tne 
model  employed. 

Tactical  decisions  can  only  be  made  witn  timely, 
pertinent  and  accurate  information.  The  specific 
information  presented  by  this  model  can  oe  extremely 
important  to  temporal  and  spatial  variation  of  ducting  and 
refraction.  As  a  bulk  model,  it  is  a  gross  approximation  to 
actual  conditions,  out  affords  decision  makers  an 
estimationof  most  probaole  relevant  cnanges  in  the 
meteorological  conditions. 
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MODEL  USAGE  PROCEDURES 

1.  Model  Components  and  Flowcnact 

Using  this  code  requires  minimal  knowledge  of 
microcomputer  systems.  Inputs  required  are  an 
initial  radiosonde  (IREPS),  forecast  winds  (up  to 
30  hours),  and  subsidence  (divergence)  and 
thermal  advection.  Because  of  forecasting,  some 
ability  to  estimate  large  scale  subsidence  and 
advection  is  required. 

1.1  Model  Components. 

A  block  diagram  and  flow  chart  for  the  program  appear  in 
Fig.  A-l  and  A-2,  respectively.  The  block  diagram  (Fig  A-l) 
gives  an  overview  of  the  program  operation  and  flow,  as 
broad  as  possible  without  specifics.  The  flow  chart  can  be 
useful  to  keep  track  of  position  within  the  execution  of  the 
program.  The  program  flow  (Fig  A-2)  was  designed  to  oe  as 
autonomous  as  possible,  yet  enable  the  user  to  interface 
with  the  precedure  as  deemed  necessary.  The  two  critical 
points  in  program  flow  are  the  input  of  subsidence  and 
selecting  an  initial,  simplified  model  for  the  atmospheric 
structure  (digitizing). 

The  following  sections  include  separate  discussion  of 
program  operation,  initial  parameter  input,  selection  of 
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output,  and  procedures  for  adjusting  inputs  to  produce 
meaningful  output. 

2.2  Overview  of  Model  Flow 

The  first  list  of  items  for  selection  is  tne  DIRECTION 
MENU  from  which  the  user  selects  the  operation  to  be 
performed.  The  operation  can  be  to: 

1)  Predict  the  change  of  the  inversion, 

2)  Display  the  predicted  results,  graphics  or  tabular 
form, 

3)  End  prediction  model  and  rewind  tapes, 

4)  Return  to  IREPS, 

5)  Review  or  Change  any  initial  conditions,  then  rerun 
the  prediction  model. 

Any  data  set  on  an  IREPS  data  tape  can  oe  called.  The 
prediction  would  normally  be  the  first  cnoice  in  a 
procedure.  In  the  prediction  procedure  (choice  1) ,  the 
computation  loop  steps  are: 

1)  Choose  the  IREPS  data  set, 

2)  Enter  forecast  winds  and  other  initial  conditions, 

3)  Approximate  a  stick  structure  to  the  temperature  and 
specific  humidity  profiles  from  the  selected  IREPS  data 
set. 

Do  not  use  IREPS  bailout  keys  to  correct  errors. 
EXCEPTION:  User  Defined  Key  1  will  taxe  the  user  back 

to  the  DIRECTION  MENU. 
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Once  all  the  computations  have  been  performed  for  a 
forecast  period,  the  program  will  return  to  the  DIRECTION 
MENU.  If  the  display  (2.3)  is  selected  at  tnis  point,  the 
user  may  select  displays  of  various  computed  values. 
Selection  may  be  of  the  inversion  heignt  cnange  and  'm' 
profiles  at  spaced  intervals  or  the  inversion  strength 
'jumps'  and  well-mixed  temperature  and  specific  humidity  or 
a  tabular  print  out  of  the  computed  parameters.  After  eacn 
of  these  displays  are  completed,  the  program  will  return  to 
the  DISPLAY  MENU  for  a  choice  of  different  display  or  return 
to  the  DIRECTION  MENU. 

The  user  may  also  use  the  DIRECTION  MENU  to  cnoose  a 
different  value  for  subsidence  or  other  parameters  by  going 
through  the  REVIEW/CHANGE  MENU.  The  prediction  model  must 
be  run  after  changes  are  made  to  produce  new  outputs. 

In  the  following  section  we  will  consider  tne  program 
steps,  the  block  diagram  should  be  used  to  follow  the 
discussion. 
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PREDICTION  FLOW  CHPRT 


Fig.  A-2 


PROGRAM  FLOWCHART 


2.  PROGRAM  EXECUTION  STEPS 


2.1.  Starting  Procedures 

2.1.1  Loading  and  Starting  Program 

The  prediction  model  program  is  loaded  oy  typing: 

LOAD  "MIXED2 : xxxM , 1  then  EXECUTE 

Where  xxx  is  the  proper  mass  storage  device  (e.g.  T14) 
containing  the  model  program. 

NOTE 

In  default,  the  program  will  look  for  the  data 
tape  to  be  in  the  left  nand  tape  drive  (T14) . 
Questions  requiring  user  response  will  be  set  off  by 
asterisks,  the  selection  menus  can  be  viewed  on  the  flow 
diagrams  or  at  the  beginning  of  each  major  section. 

2.1.2  Enter  Execution  Date 

************************************************************ 

ENTER  Execution  Date  as  date  of  forecast  run,  just  press 
CONT  if  the  same 

************************************************************ 

The  entered  date  will  be  printed  on  all  displays  to 
indicate  the  date  of  run,  this  is  not  necessarily  the  date 
of  the  data  set.  If  it  is  left  blank,  the  date  will  default 
to  that  of  the  IREPS  data  set  (i.e.  IREPS  date  if  IREPS  has 
just  been  run,  blank  if  this  program  is  run  alone.)  The 
computer  will  then  read  the  IREPS  data  table  and  display  the 
DIRECTION  MENU. 
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DIRECTION  MENU 


2.2  DIRECTION  MENU 


Cnoices  of  procedures  are: 

1)  COMPUTE  THE  LAYER  CHANGE 

2)  DISPLAY  THE  OUTPUTS 

3)  END  THE  PREDICTION  MODEL 

4)  RETURN  TO  IREPS  PROGRAM 

5)  REVIEW/CHANGE  ANY  INITIAL  CONDITIONS 

ft*********************************************************** 

ENTER  option  (1  to  5)  then  press  CONT 

************************************************************ 

2.2.1.  COMPUTE  THE  LAYER  CHANGE. 

The  list  of  data  sets  on  the  IREPS  data  tape  is  presented 
and  the  selection  is  made  by  entering  the  corresponding  line 
number,  followed  by  CONT. 

2. 2. 1.1.  Data  Set  Choice 

A*********************************************************** 

ENTER  THE  NUMBER  OF  THE  DATA  SET  DESIRED 

************************************************************ 

If  the  desired  data  set  is  not  on  the  present  data  tape, 
hence  is  not  displayed,  SFK  3  (kl — top  right  of  keyboard) 
should  be  used  to  return  to  the  DIRECTION  MENU,  or  type 
'CONT  Option',  then  EXECUTE.  If  a  different  data  tape  is  to 
be  used,  insert  it  in  the  left  tape  drive  and  restart  the 
program  with  STOP  and  RUN.  Otherwise,  the  program  will 
retrieve  the  entered  choice  of  data  set  and  proceed  toward 
the  next  procedure  wnicn  is  digitizing.  The  program  will 
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ask  for  winds  (2. 2. 2. 4.),  subsidence  (2. 2. 1.5.)  and  tbe  sea 
surface  temperature  (2. 2. 1.6.). 

2.2. 1.2  Latitude  and  Julian  Day 

The  program  uses  these  two  inputs  for  calculating  the 
zenith  angle  of  the  sun  and  thus  the  amount  of  insolation  at 
the  top  of  the  inversion  or  cloud  layer. 
************************************************************ 

ENTER  Latitude  of  the  prediction  area  [-90  to  90J  [use  - 
for  South] 

a*********************************************************** 

ENTER  the  Julian  day  of  the  year  [maximum  1  digits] 

ft*********************************************************** 

2. 2. 1.3  Beginning  Time 

a*********************************************************** 

ENTER  THE  BEGINNING  TIME  [4  digits] 

***************  **  **  ***************************************** 

The  beginning  time  will  be  displayed  on  all  the  output. 
The  entered  time  snould  be  fairly  close  to  tne  start  time 
since  it  determines  the  beginning  of  the  30  hour  forecast 
period.  The  entered  time  must  have  four  characters  (e.g. 
1900) .  The  program  displays  will  only  be  in  whole  hour 
segments  for  easier  graphic  interpretation. 

2. 2. 1.4  wind  Forecast 

*********************************************************** 

IS  THIS  rfIND  FROM  IREPS  OKAY  (Y/N)  ? 

*********************************************************** 
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The  program  will  display  tne  initial  time  and  tne  wind 
speed  as  previously  entered  in  tne  I  REPS  data  set.  If  tne 
wind  is  to  scay  tne  same,  press  COST  wit n  no  entry,  tne 
response  will  derault  to  the  first  of  tne  options.  If  tne 
wind  is  to  be  cnanyed,  enter  'N ' ,  press  CONT,  in  wnicn  case 
the  following  question  will  as*  for  a  new  beginning  wind 
speed. 

************************************************************ 

ENTER  dEJ INNING  SIS D  SPEED  [KNOTS!  at  time  XXXX 


Winds  are  requested  in  KNOTS.  The  code  performs 
computations  with  wind  in  m/s,  a  conversion  is  performed, 
(c.f.  Appendix  0.)  The  display  is  now  set  up  to  enter 
forecast  winds  for  tne  period.  (Fig  3-1) 


flhmm,ww  FROM  THE  FORECAST,  ENTER s  TIME, WIND  (KNOTS) 

************************************************************ 

Tne  maximum  numoer  of  winds  J3ED  is  10.  When  done 
entering  winds,  enter  'e'  or  'END'.  Winds  will  be  linearly 
interpolated  oetween  input  times.  This  is  a  gross 
approximation,  but  worxs  well  with  tne  degree  of  accuracy 
used  in  forecasting  winds.  To  get  a  smoothing  closer  to  tne 
forecast,  enter  more  winds  at  the  times  of  rapid  cnange.  If 
tne  fine  structure  is  availaole  or  needed  only  for  tne  first 
part  of  the  period,  enter  the  winds  as  desired.  If  tn-e  last 
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wina  entered  is  ror  a  time  oefore  tne  end  or  tne  oeriod  it 


will  oe  neld  constant  to  tne  enu  of  tne  period. 

If  the  initial  wind  is  to  oe  used  througnout  tne  period, 
enter  '2 '  or  'END'  as  cne  first  response  to  tnis  question. 

As  a  prompt,  the  time  of  tne  last  possiole  input  period,  30 
nours  after  tne  oeginning  time,  is  displayed.  All  times 
must  be  entered  sequentially,-  times  beyond  tne  last  period 
will  not  be  considered. 


00  TuESE  WINDS  LOOK  OKAY  ( Y/N ) ? 


If  different  winds  are  desired,  the  computer  will  as* 
for  correction  of  entered  winds,  one  line  at  a  time.  If  tne 
displayed  winds  are  correct,  press  CQNT,  if  not,  enter  .1  and 
refer  to  the  previous  qestion.  Tne  next  requests  will  oe 
for  suosidence,  and  a  cnecx  of  sea  surface  temperature.  The 
program  flow  will  tnen  proceed  to  digitizing. 

2.2. I. 5;  Subsidence  Forecast 

************************************************************ 


ENTER  SUBSIDENCE  (m/sec)  [First  value*-. 003  ],  or  CQNT 
to  use  same  value? 

************************************************************ 

Enter  a  value  of  subsidence  in  meters  per  second.  This 
can  oe  calculated  from  tne  divergence  field,  determined 
explicitly  from  tne  previously  used  value  or  by  manipulating 
tnis  model  oetween  successive  soundings. 


2. 2. 1.6.  Sea  Surface  Temperature  Specification 

******#*iH**iHHt**iH****ifr****i»******#**iHt»******llr****i«*#**##*#* 

ENTER  THE  SEA  SURFACE  TEMP  IN  'C  or  just  CO NT  for  same 
value 

ft*********************************************************** 

The  sea  surface  temperature  from  IREPS,  wnicn  is 
displayed,  may  need  refinement  from  tne  initially  entered 
value.  This  would  arise  if  the  temperature  used  in  tne 
IREPS  data  set  was  the  sea  water  injection  temperature.  To 
use  the  existing  value,  just  press  CONT  without  entering 
anything.  To  cnange  tne  value,  enter  tne  new  temperature  in 
degrees  centigrade,  and  press  CONT. 

2. 2.1.1  Profile  Initialization 

The  program  will  now  go  to  the  digitizing,  the 
determination  of  the  vertical  structure  of  potential 
temperature  and  specific  humidity.  The  steps  establisn: 
Inversion  neight 
Potential  temperature 

(well-mixed,  jump  and  lapse  rate  values) 

Specific  humidity 

(Nell-mixed,  jump  ana  lapse  rate  values) 

Previous  inputs  will  be  disiayed  on  the  plot.  One 
display  with  data  and  digitized  results  appears  in  Fig.B-2. 
At  the  end  the  user  will  be  asxed  to  specify  wnether  tne 
plotting  is  OKA 'i  or  if  he  wants  to  do  it  AGAIN,  if  an  error 
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is  made  during  the  digitizing  routine,  continue  to  tnis 
question,  then  redigitize. 

Tne  program  will  draw  out  tne  digitizing  graph  and  its 
labels,  then  enter  tne  digitizing  mode,  rne  first  structure 
to  plot  will  be  the  potential  temperature,  tne  second  will 
be  the  specific  humidity.  The  profile  data,  taxen  from  tne 
IREPS  data  tape,  will  be  displayed  to  2400  meters  or  the 
maximum  neignt  of  entered  data  if  less  than  2400  meters. 

Nnen  digitizing  can  begin,  the  computer  will  peep  and  a  full 
screen  cursor  will  appear.  To  position  the  cursor,  use  the 
arrow  keys  at  the  top,  center  of  the  keyboard,  tnen  CONT. 

NOTE 

Finer  structure  is  achieved  by  using  snift  when 
positioning  the  cursor.  This  will  move  it  one 
pixel,  or  picture  element,  at  a  time. 

2. 2. 1.7.1.  POTENTIAL  TEMPERATURE:  3 

la)  Inversion  Height:  position  the  horizontal 
line  at  a  representative  inversion  neight 

(usually  near  the  middle  of  the  sharp 
temperature  'jump'.) 

lb)  well-mixed  Temperature:  position  the 
vertical  line  at  a  representative  value  for 
the  well-mixed  temperature  below  the 
inversion,  then  CONT.  POINT  3 

2)  Inversion  'jump':  horizontally  position 
the  small  cursor  at  the  point  wnere  tne 
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extension  of  the  profile  above  the  layer  meets 
tne  norizontal  'jump'  line  at  the  inversion 
neignt  (use  precision  eyeball),  then  CONT. 

POINT  2 

3)  Gradient  Above  the  Inversion:  position  the 
small  cursor  at  the  top  of  the  screen  or 
heignest  data  point  to  get  the  best  linear 
approximation  of  potential  temperature  aoove 
inversion,  determining  lapse  rate,  then  press 
CONT.  POINT  3 

2.2. 1.7. 2. _ SPECIFIC  HUMIDITY:  J 

la)  The  inversion  height  determined  for  the 
potential  temperature  sequence  will  be  the  same 
for  the  specific  humidity  procedure.  POINT  4 

lb)  well-mixed  Specific  Humidity:  position  the 
vertical  line  at  a  representative  value  for  the 
structure  below  the  inversion, then  CONT.  POINT  4 

2)  Inversion  'jump':  position  the  cursor  where 
tne  extension  of  the  profile  above  meets  the 
norizontal  'jump'  at  the  inversion  (again,  this 
is  very  subjective) ,  then  CONT.  POINT  5 

3)  Gradient  Above  the  Inversion:  position  tne 
small  cursor  at  the  top  of  the  screen  or  highest 
data  point  to  get  the  best  linear  approximation 
of  the  specific  humidity  above  the  inversion, 
determining  lapse  rate,  then  press  CONT.  POINT  6 
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Fig.  A-4.  EXAMPLE  OF  DIGITIZED  VERTICAL  'STRUCTURE 
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Fid.  A-5.  EXAMPLE  OF  DIGITIZED  '/EKTICAL  STRUCTURE 
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If  the  digitized  plot  is  not  as  desired,  move  the  cursor 


to  AGAIN,  press  CONT,  and  begin  with  POTENTIAL  TEMPERATURE 
(2.2.1. 6.1).  If  tne  digitized  plot  is  agreeaole,  move  tne 
cursor  to  the  OKA*  and  press  GO NT.  This  will  start  tne 
program  into  the  computing  stage.  A  picture  of  tne  layer 
computatutions  is  progressively  displayed. 

At  the  completion  of  the  run  the  program  will  return  to 
tne  DIRECTION  ME  MU .  Tne  running  prediction  may  be  PAUSEd  at 
any  time,  then  kl  or  'CONT  Option',  'EXECUTE'  will  taxe  tne 
program  bacx  to  the  DIRECTION  MENU. 

2. 2. 1.8.  Change  Subsidence  or  Other  Initial  Parameters 

A  decision  point  exists  nere  as  to  the  reasonableness  of 
the  layer  movement.  If  it  looks  good,  choose  to  display  tne 
change  (2.2.2.).  If  the  results  of  tne  computation  do  not 
go  as  expected,  cnoose  to  re-compute  using  a  new  value  of 
subsidence  or  other  initial  data,  or  redigitize  the  vertical 
structure.  3oth  can  be  performed, by  cnoosing  RE VIEW /CHANGE 
MENU  (2.2.5.).  Use  of  the  flow  cnart  can  be  aelpful  wnen 
running  the  program. 

2.2.2.  Display  tne  Outputs 

Selection  of  this  option  will  take  the  user  to  the 
DI3PLAT  MENU  for  more  defined  and  varied  displays  of 
computed  inversion  layer  movement: 

1)  overall  display  of  inversion  movement  and  .l-orofiles 

2)  tne  inversion  jumps  of  temperature  and  humidity 


o  2 
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3)  tabulation  of  selected  values  tnrougn  tne  period 
Refer  to  paragraph  2.3  DISPLAY  .-l£NU 

2.2.3.  End  Prediction  .lodel  and  Rewind  the  Tapes 

This  choice  will  simply  rewind  the  tapes  in  tne  computer 
and  end  the  prediction  model  program.  Remember  to  use  tne 
EJECT  bar  to  remove  tapes. 

2.2.4.  Return  to  IREPS 

This  choice  will  taKe  the  user  back  to  tne  IREPS 
program.  Tne  IREPS  driver  tape  must  be  in  the  right  nand 
tape  drive,  before  pressing  CONT,  to  load  successfully. 

2.2.5.  REVIEW/ CHANGE  ME MU 

This  option  is  for  rerunning  the  program  witn  different 
values  for  initial  data.  This  can  also  be  used  for  going 
through  the  digitizing  routine  and  retaining  the  winds  tnat 
were  input.  Wnen  the  end  of  the  digitizing  routine  is 
reacned,  the  user  will  again  be  given  the  choice  of  changing 
the  digitized  points.  If  all  looks  satisfactory,  press  CONT 
whicn  will  return  the  program  to  this  menu.  This  procedure 
can  be  continued  until  a  good  approximation  is  generated. 
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Fig.  A-S.  DISPLAY  '*ENU  FLOWCHART 
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2.3.  DISPLAY  MENU 


Tne  following  options  will  oe  listed: 

1  DISPLAY  OP  INVERSION  AND  *M'  PROFILES  grapnics 

2  'JUMP'  STRENGTH  AT  INVERSION  graphics 

3  LIST  OF  COMPUTED  VALUES  taoular 

4  RETURN  TO  THE  DIRECTION  MENU 
*********************************************************** 

ENTER  display  choice  tnen  press  CONT 

************************************************************ 

The  user  will  be  asxed  of  he  wants  naru  copy,  after  tne 
grapnics  displays  are  drawn  and  before  the  tabular  listing 
is  run.  Option  4  will  take  the  program  nack  to  tne  DISPLAY 
MENU  as  will  SFK  1  at  any  time.  Further  discussion  and 
interpretation  of  the  various  aspects  of  the  displays  occurs 
in  Appendix  3. 

2.3.1.  Graphics  Display  of  Inversion  and  'M '  Profiles 
As  the  main  graphic  product  of  the  prediction,  this 
display  depicts  the  cnanges,  over  tne  30  hour  period,  in  the 
inversion  height,  the  lifting  condensation  level  (for  cloud 
formation),  the  relative  humidity  and  the  'A '  profile 
(modified  index  of  refraction)  indicating  elevated  and 
surface  based  duct  conditions. 

2.3.2  Grapnic  Display  of  Jumps  and  Nell-Mixed  Values  of 
Potential  Temperature  and  Specific  Humidity 

To  illustrate  the  change  in  inversion  strength,  tnis 
display  graphically  snows  the  variations  in  tne  well-mixed 
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values  next  to  the  'jump'  values  of  ooth  potential 
temperature  and  specific  humidity.  Also  displayed  is  tne 
relative  humidity  as  it  evolved  through  the  forecast  period. 
2.3.1  Tabular  Display  of  Selected  Parameters 

This  option  lists  tne  values  that  are  stored  during 
computation  to  be  used  for  generating  grapnic  displays.  They 
are  useful  to  more  fully  understand  tne  movement  of  the 
inversion,  the  cnange  in  tne  duct  area  and  tne  inversion 
streng  th. 

2.3.4  Return  to  tne  DIRECTION  MENU 

This  option  will  return  tne  program  to  the  DIRECTION 
MENU  for  further  cnoices  to  compute  again,  return  to  tne 
IREPS  program  or  just  end  the  prediction  model. 
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? ij.  A- 6.  R2VIEW/CHANGE  MENU 
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2.4  REVIEW/CHANGE  MENU 

Select  one  o 1  the  following  to  cnange  a  parameter  after 
whicn  the  prediction  needs  to  ran  again: 

1)  List  Current  Initial  Values 

2)  Change  Wind 

3)  Change  Subsidence 

4)  Change  Sea  Surface  Temperature 

5)  Change  Latitude  or  Julian  Day 

6)  Redigitize  tne  Vertical  Structure 

7)  Select  a  different  data  set 
3)  Return  to  the  DIRECTION  MENU 

Once  this  option  is  exercised,  the  model  must  run 
througn  another  prediction  before  drawing  out  any  displays. 
*********************************************************** 

ENTER  choice  from  RE VI EW/CJANGE  .4ENCJ  (1-8) 

*********************************************************** 

2.4.1.  List  Current  initial  Values 

This  will  list  the  values  of  the  indicated  variaoles 
that  were  last  used  or  modified.  The  user  will  be  preompted 
to  get  a  hard  copy.  Default  of  the  question  will  be  to  tne 
CRT,  press  CONT  wnen  ready. 

2.4.2.  Change  Wind 

Refer  to  procedure  2. 2. 1.3  and  2. 2. 1. 4  and  the  return  to 


REV IEN/CHANGE  .1ENJ 


2.4.3 


nange  Subsidence 


Refer  to  procedure  2. 2. 1.5,  return  to  REV ISW/CHANGE  .  1EMG 

2.4.4.  _ C.nange  Sea  Surface  remperatur e 

Refer  to  procedure  2. 2.!. 5,  return  to  REV  I EW/ CHANGE  M£:iU 

2.4.5.  Change  Latitude  or  Julian  Day 

Refer  to  procedure  2. 2.!. 2,  return  to  R E V I E ;i /CHA. !G E  . IE l< J 

2.4.5.  Redigitize  the  Vertical  Structure 

Refer  to  procedure  2. 2. 1.7,  return  to  RE VI EW /CHANGE  MENU 

2.4.7,  Select  a  Mew  Data  Set 

Choosing  a  new  data  set  will  taxe  the  program  tnrougn 
all  the  initial  steps  to  set  up  the  necessary  parameters  to 
make  a  prediction  model  run. 

2.4.8.  Return  to  PI RECTI CM  MENU 

This  option  will  return  tne  program  to  tne  DIRECTION 
MENU  for  further  choices  to  compute  the  forecast  again, 
return  to  the  IREPS  program  or  end  the  prediction  nodel. 
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APPENDIX  B 


PREDICTION  OUTPUT 

This  appendix  contains  examples  of  the  prediction  output 
wnich  are  obtained  by  making  selections  from  tne  DISPLAY 
MENU  for  output.  Two  different  prediction  runs  are  snown, 
eacn  for  a  different  observational  period.  The  two  periods 
were  a  cloud  covered  case,  previously  discussed  in  the  body 
of  the  thesis,  and  a  clear  sky  case,  described  by  Davidson 
et  al  (to  be  submitted,  JAM,  1982). 

Entrainment  energetics  for  the  cloud  covered  case  in 
this  appendix  were  based  on  fomulations  by  Stage  and 
Businger  (1931).  Results  for  the  same  period  were  shown  in 
the  body  of  the  thesis,  but  were  based  on  a  formulation  by 
Deardorff  (1976).  The  scope  of  this  thesis  was  not  to 
compare  the  differences  between  the  entrainment 
parameter izat ions.  Such  a  comparison  will  be  performed 
separately  by  Naval  Postgraduate  School  investigators. 

Prediction  output  examples  appear  in  two  sets,  one  for 
each  prediction  period.  Each  set  contains  a  listing  of 
several  computed  variables  and  two  pages  with  graphic  plots 
of  several  computed  variables  and/or  profiles.  The  plotted 
variables  do  not  encompass  all  those  listed.  All  selected 
outputs  automatically  appear  on  the  CRT,  and  the  printed 
copies  (hard  copy)  may  be  obtained  by  selecting  the 
appropriate  option  after  CRT  display. 


The  listing  of  predicted  results  (Fig.  B-3  and  3-6)  are 
self-explanatory.  The  listed  values  are  for  every  30  minute 
period  (the  prediction  time  step)  and  the  units  are: 
velocity-m/s,  temperature-Kelvin,  neight-ra  and 
humidity-g/kg . 

The  graphic  plots  contain  information  about  the 
predicted  vertical  structure  (Fig.  B-l  and  B-4)  and  tne 
predicted  values  in  the  mixed-layer  and  at  tne  inversion 
(F  ig  .  B-2  and  3-5 )  . 

The  predicted  vertical  structure  is  presented  in  two 
panels  in  Fig.  B-2  and  Fig.  B-4.  The  top  panel  shows 
predicted  M  profiles  (Modified  index  of  refraction)  at  6 
hour  intervals,  including  possible  surface  based  duct 
conditions.  It  is  determined  from  the  predicted  well-mixed  T 
and  q  and  the  predicted  jumps  at  and  lapse  rates  above  the 
inversion.  The  bottom  panel  shows  the  top  of  the  well  mixed 

layer  ( - ) ,  the  top  of  the  inversion  zone  (-  -  -) 

assuring  an  inversion  thickness  of  .1  layer  height  and  the 
lifting  condensation  level  (.....).  The  lifting  condensation 
level  was  computed  from  the  well  mixed  temperature  and 
specific  humidity  values.  A  description  of  duct 
determination  appears  at  the  bottom  of  the  hard  copy  page. 

The  graphic  display  of  predicted  values  (Fig.  B-2  and 
Fig.  B-5)  in  the  mixed  layer  and  at  the  inversion  are  of 
several  variables  appearing  in  the  listing  option.  The 
solid  lines  in  the  top  panel  represent  the  well  mixed  values 


in  the  layer  appropriate  for  tne  10  meter  level.  The  dashed 
lines  are  the  'jump'  of  that  parameter  (T  or  q)  at  the 
inversion.  The  'jump'  magnitudes  are  important  in 
determining  the  strength  of  the  duct  and  its  depth.  The 
bottom  panel  shows  the  relative  humidity  at  the  10-meter 
level.  It  can  be  used  to  determine  the  proximity  of  fog  or 
clouds  and  would  be  important  in  aerosol  equilibrium  size 
determination. 
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DUC  T  HE  I GHT  (meters' 


This  i*  i  SIMPLIFICATION  of  *.hf  real  |tr«jc*yrt, 

Thi»  display  i  s  divided  into  VPPgR  and  LC'UgP  vindovs: 

LONER  UINDOU  display!  the  top  and  middle  of  ♦  he  elevated  l  kvir  and  ms 
forecast  continued  i»out»tnt  for  4  thirty  hour  ptnod  at  *  er  beginning. 
At  the  bottom  is  the  Relative  Humidity  it  *4Ch  su  hour  period. 

The  tightly  dotted  line  is  the  lifting  •:  ondensat  i  on  level. 

CMith  enough  moisture  near  the  inversion,  this  can  be  used 
4S  4  flsg  for  possible  cloud  for«4tion,  or  for  high  mixed- 
l  4ver  humidity,  near  the  surf  see,  to  forecut  fog.  1 
JPPgR  HlNDQM  picks  out  fl-vslue  structure  using  only  4  points:’ 
st  each  six  hour  period  snd  will  indicate  4  surf  see  bssed  duct 
>n  I  y  if  the  el  ousted  Pl-uslut  is  less  than  or  uithm  ?  of 
the  surf  see  M-value.  The  sampled  times  sre  displayed  st  top. 


Fi q.  a-l.  Vertical  Profile  and 


Profiles-  Cloudy  Sky 


gff J»ggfl  TU*g  >nd  HUW1PITV  for  NPS  33  29 /T3. 0*30 


Thu  1*  k  1  i  *».  of  plotted  HP? 

HOUR  Zi  Z  J  ■: !  R»IHu«  HIND  Th 
0500  326  330  70.7  3.7  12. i 

0330  330  332  30.3  3.3  12.0 

0600  334  330  31.4  3.3  12.0 

0630  53?  312  31.3  7.3  U.a 

0700  544  301  32.1  7.!  u.  t 

0730  34?  2?3  32.1  7.2  12.1 

0800  333  237  31.3  7.3  12.2 

0830  560  301  81.3  7.4  12.4 

0000  336  310  31.1  7.3  12.3 

8930  371  320  30.6  7.3  12.7 

1000  377  333  08.1  7.7  12.0 

1030  382  346  79.3  7.3  13.0 

1100  387  339  79.0  0.0  13.2 

1130  591  371  70.3  3.1  13.3 

1200  393  382  79.2  9.2  13.3 

1230  300  392  77.9  3.4  10.6 

1300  303  400  77.6  3.3  13.7 

1330  609  406  77.3  3.3  13.3 

1400  614  410  77.4  0.0  13,3 

1430  318  412  77.4  3.3  14. 0 

1500  323  411  77,5  9.8  14.0 

1330  527  409  77.7  9.1  14. 0 

1680  532  403  79.0  9.2  14.0 

1530  638  399  79.3  9.0  14.0 

1730  343  392  73.5  3.7  14.0 

1730  649  383  79. 0  3.4  13.  9 

1800  633  374  79.4  3.  1  13.3 

1330  662  366  79.7  7.3  13.8 

1900  568  333  79.9  7.?  ^3.7 

1930  675  332  30.1  7.2  1  2.5 

2000  632  348  30.2  5.9  13.3 

2030  689  343  38.3  6.6  13.3 

2100  696  343  30.4  6.3  13.3 

2130  702  341  30,4  6.0  13.4 

2200  709  340  30.4  3.7  13.4 

2230  713  339  30.4  3.4  13.4 

2300  721  339  30.4  5.2  13.3 

2330  727  339  30.4  3.2  13.3 

0000  733  333  30.3  3.2  1  3.2 

3030  739  333  30.3  5.2  13.2 

0100  743  338  30.3  3.2  13.1 

0130  731  337  30.3  3.2  12.1 

0200  735  337  30.3  3.2  12.1 

0220  762  337  30.  3  3.2  1  3.0 
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The  pom  bi  1  1 1 v  of  a  *Mrf  ic«  based  duct. '  oBD'  is  indicated  and  only 
due  to  the  'M'  v«lu*  greater  at-  t  h*  jurfacf  than  at  th*  inversion 
height.  Variation  in  *  h*  vertical  structure  nuit  be  tMohui:*d 
••/hen  knifing  t  his  forec  ast  m  that  a  oBD  wav  not  exist. 

This  ia  a  •SIMPLIFICATION  of  *h*  real  structure. 

This  display  is  divided  into  UPPER  and  LQUER  vindous: 

LONER  NINDQM  displavs  th*  top  and  middle  of  the  *  legated  U"*r  an©  it* 
forecast  continued  movement  for  a  thirty  hour  period  after  3**5 inning. 
At  the  bottom  1  j  the  Relative  Hu»i*jit"  at  each  si  -  hour  period. 

The  lightly  dotted  lin*  13  * n*  lifting  condensation  e**’. 

CUith  enough  moisture  near  ’h*  inversion,  this  can  be  used 
as  a  flag  for  possible  cloud  formation,  or  for  high  i»i  <#4- 
1  aver  humidity,  near  the  surface,  to  forecast  fog.  I 
UPPER  NINDOM  picks  Out  M-value  structure  using  onl*>  4  points) 
at  each  six  hour  period  and  will  indicate  a  surface  based  duct 
only  if  the  elevated  M-val  ue  <s  less  than  or  vithin  5  of 
the  surface  M-val  ue.  The  sampled  times  are  displayed  at  *op. 


ig.  3-4.  Vertical  Profile  and  'cl''  Profiles-  Clear  SKy 
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TEMPERATURE  lod  HUMIDITY  iMr  SCI  Id.  03^?6-  14i30^SCI 


Mixed  Layer  and  'Jump'  \/alues-Clear  Sky 
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Fig.  3-6 


APPENDIX  C: 


TERMS  USED  IN  THE  DISCUSSION  AND  PROGRAM: 


Cn2  ...  Index  of  refraction  structure  function  parameter 
[m  -2/3] 

Extinction  ...  loss  of  EM/EO  energy  due  to  absorption/ 

scattering  by  species  (molecules/aerosols)  in  the  path. 

Lifting  Condensation  Level  ...  Level  at  wnicn  a  parcel  will 
reach  saturation  if  cooled  while  undergoing  lifting. 

Mixed  Layer  ...  That  part  of  the  boundary  layer  wnere 
turbulent  mixing  has  destroyed  the  gradients  of 
conservative  quantities. 

Mixing  Ratio  ...  the  ratio  of  the  water- vapor  density  to  the 
density  of  the  dry  air. 

Potential  Temperature  ...  the  temperature  that  an  air  parcel 
would  assume  if  its  pressure  were  changed  by  a 
dry-adiabatic  process  to  some  standard  pressure. 

Relative  Humidity  ...  the  ratio  of  vapor  pressure  to 
saturation  vapor  pressure  times  100.  Vapor  in 
equilibrium  with  sea  water  is  assumed  to  have  a 
relative  humidity  of  93  per  cent. 

Specific  Humidity  ...  the  ratio  of  water-vapor  density  to 
total  density  and  i~:  usually  designated  by  q. 

Subsidence  ...  a  sinning  of  air  from  nigh  to  low  levels.  The 
causes  of  its  occurrence  are  dynamic,  such  as  tne 
frictional  outflow  of  air  near  the  surface  associated 
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C-l 


1 


witn  nign-pressure  systems,  the  divergence  of  air  from 
near  the  surface  nign  pressure  centers,  and  the 
divergence  of  air  due  to  change  in  latitude  as  an  air 
mass  moves  from  north  to  soutn.  (conservation  of 
absolute  vorticity) 

Thermal  Wind  ...  vector  wind  cnange  with  height  due  to 
horizontal  temperature  gradient.  This  implies 
advection  above  the  inversion. 

Conversion  ...  KNOTS  *  (m/s)  *  1.94 

1.94  *  ft  *  Nm  *  3600 sec 


.  3048m  6076.1ft  Hour 


5  i  *****************  ************************************  * 

10  i  ***  FOR  TIE  TO  IREPS  DATA  PREDICTION  *** 

15  !  ***  LT  DAVID  A.  BROWER,  USN  ,  NPS,  MONTEREY ,  CA .  *** 

20  !  ****************************************************** 
25  i  BASIC  DAVIDSON  MODEL  ".MIXED 2  '  MARCH  30 

30  !  REVISED  CHRIS  FAIRALL  ENTRAINMENT  COMPUTATION  11JAN32 
3!  !  NOW  INCLUDING  LONG  AND  SHORT  WAVE  RADIATION 
1  3  J  AN  8  2 

35  EXIT  GRAPHICS  !  UPDATED:  19  MAR  82 

40  GCLEAR 

45  OPTION  BASE  l 

50  NORMAL 

55  OFF  ERROR 

60  PRINTER  IS  0 

65  PRINT  CHRS (27) &"&132T" 

70  PRINTER  IS  16 

75  Mass$=":T! 4"  J  CHANGE  HERE  FOR  A  NEW  MASS  STORAGE 

80  GOTO  155 

35  !  THE  FIRST  LINE  FOR  'COM'  NEEDED  ONLY  FOR  2  ..0 

90  !  COM  I r pr og$ , I r dataS , Rev$ , Date  S  >  Pr int  s$ , Per f or a tea , 

Date  wanted 

95  COM  Names [24] ,LocS[24] , TimeS [24] ,TyoeS [ 1] , Heights [1 ] , 

Evao$ [1 ] 

100  COM  WmoS (30) [11] ,WmontS [ 5J 

105  COM  SHORT  Presur ( 30 ) , Temper ( 30 ) , Relhum ( 30 ) 

110  COM  SHORT  Height  (  30 )  ,Muni  ts(  30  )  ,Nuni  ts  (  30  ) 

115  COM  SHORT  Wind ,  Sea  ,A  ir  , Relhm  ,H  tzer o,  Pr  zero ,  Nmax ,  De  1  ta , 

Change, Recenv 

120  COM  SHORT  Coml  (  32 , 4)  ,Envsq  ( 1 6)  , Syssg (  3 2)  , Lossq  {  32) 

] 25  COM  Envnam$( 16) [24] ,SysnamS(32) [24] ,Losnam$(32) [24] 

130  !  COM  Over lay$[ 24 ], Options [ 80] 

135  !  COM  Nrec ,Unorotected 
140  !  DIM  DummyS [ 160] ,Errm$[80j 

145  !  DIM  Dlabel$[ 80] , Defaults [ 80] ,Doptions$ [80] 

150  !  DIM  Next  overlay$[24] 

155  DIM  WindTn (1 0, 2) , Sticks ( 3) , Frame  locS [24] , Frame_date$ [ 
24] 

160  INTEGER  J 

165  DIM  Pot  temo( 30) , Spec  hura( 30) , Zsnd ( 30) , LS [ 1 60] 

170  Julday*T50 

175  Lat»33 

180  SHORT  G (63 ,7) 

185  RESTORE  190 

190  READ  Z  ,Zot,Kkk,  Alpha_t,Deltim, Sticks  (*)  !  Alpha  t  is 
the  ratio  of  heat”to 

195  DATA  10 ,5E-5  ,  .35, 1  .  35 ,  1800  !  momentum  transfer 

at  Psi*0  [from  BULK] 

200  DATA  "Gradient  above" ," Inver sion  jump","  well- 
mixed" 

205  Wsl=WS*-.003  i TO  SET  A  TYPICAL  DEFAULT 

210  PRINT  PAGE, "THIS  IS  'MYREP3 '  WHICH  MODELS  THE  CHANGE 
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IN  THE  MIXED  LAYER" 

215  PRINT  TAB (1 8) , "ABOVE  THE  SURFACE";  ! 

REV: 18NOV80 

220  Start:  Tday$*Date$ 

225  IF  (TdayS*"")  OR  (TdayS*"  ")  THEN  Tday$=“  6MAR82" 

230  DISP  CHR$  ( 7)  ;  "ENTER  Today's  Date,  THEN  Dress 
'CONT'"  ; 

235  INPUT  TdayS 

240  PRINT  TdayS 

245  GOSUB  Ireps  data 

250  I  MINIMUM  FREQ  TRAPPED  IN  HERTZ*  3 . 603 3E2 1 *a < - 3/2 )  d 
DUCT  THICKNESS  IN  m. 

255  Option:  !  TO  WORK  WITH  THE  IREPS  KEYS 
260  Menu:  PRINTER  IS  16 
26  5  Aa=Q 

270  PRINT  USING  "g  ,  K  ,  2/,  K "  ;  "This  is  the  DIRECTION  MENU", 
"Select  option  below,  ENTER  NUMBER  and  PRESS 
'CONT'. " 

275  PRINT  USING  "#,2/,2  (!0X,K,2/)  1  COMPUTE  NEW- 

LAYER  (displayed  as  computed  on  CRT)", "2  PLOT  of 
LAYER  CHANGES  " 

280  PRINT  USING  " 3 ( 1  OX, K, 2/) , K" ; " 3  END  'MIXED2  ' 

Program”," 4  RETURN  TO  IREPS  OPTION  LIST  PROGRAM  ( 
with  IREPS  driver  in  T15)"," 5  REVIEW/CHANGE  MENU" 
285  INPUT  “Enter  Option (1  to  5)  then  press  CONT",Aa 

290  IF  NOT  (  (Aa*l)  OR  (Aa=2)  OR  (Aa=3)  OR  (Aa*4 )  OR  (Aa=5) 

)  THEN  Menu 

295  IF  Aa*4  THEN  Menu  !  LOAD  " IREPS : T15 ", Input 

1  PUT  IN  LATER 

300  ON  Aa  GOTO  Start_out, Picture , Bye , Menu ,Revcnange 
(Subsidence 
305  GOTO  Menu 

310  Start_out:  IF  Flag_rev  THEN  Run 

315  Start_out_again:  GOSUB  Disp_data 

320  GOSUB  Julia t 

325  GOSUB  Enter_winds 

330  GOSUB  Subsidence 

335  GOSUB  Sea  temo 

340  GOSUB  Pit” struct 

345  GOTO  Run 

350  Revchange:  PRINTER  IS  16  !  TO  ENTER  ANY  CHANGES 

OR  REVIEW  DATA 
355  Rc$="0" 

360  IF  NOT  Flag  data  THEN  Menu  !  RETURN  TO  DIRECTION 
MENU  IF  NO  DATA  SET 

365  Flag  rev*]  J  TO  FLAG  A  RUN  THROUGH  THE  REVISION 

ROUTINE 

370  PRINT  USING  "3,K ,/" ; "This  is  the  REVI EW /CHANGE  MENU" 

375  PRINT  USING  " I,/, 4(1  OX, K, 2/)  " ; " 1  LIST  CURRENT 

INITIAL  VALUES", "2  CHANGE  WIND", "3  CHANGE 
SUBSIDENCE", "4  CHANGE  SEA  SURFACE  TEMPERATURE" 
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330  PRINT  USING  M t , / , 3 ( 1  OX , K , 2/ ) H ; " 5  CHANGE  JULIAN  DA* 

or  LATITUDE", "6  REDIGITIZE  VERTICAL  STRUCTURE" 7 
SELECT  NEW  DATA  SET" 

385  PRINT  USING  "10X,K";"8  RETURN  TO  DIRECTION  MENU" 

390  INPUT  "Enter  Option  (I  to  8)  tften  press  CONT",RcS 
395  PRINT  PAGE 

400  IF  (NUM(RC$) <49)  OR  (NUM (Rc$) >56 )  OR  ( LEN (  Rc$)  < >1 ) 

THEN  Revcnange 

405  IF  VAL(Rc$)»7  THEN  Start  out  again 
410  IF  VAL(Rc$)*8  THEN  Menu 

415  ON  VAL ( Rc$)  GOSUB  List_c ur r , Enter_winds , Subsidence, 

Sea  temp ,  Julia  t  ,P  Instruct,  Disp_data,Menu 
420  GOTO  Revchange 
425  List  curr:  GOSUB  Hard_output 

430  IF~Hd$="*"  THEN  PRINTER  IS  0 

435  PRINT  "THESE  ARE  THE  LATEST  VALUES  FOR  THE  PREDICTION 
MODEL " 

440  PRINT  USING  "/,K"  ;  "DATA  SET:  "S,Loc$&"  "STimeS 
445  PRINT  USING  "/K , X 3D, 2 X9A  ,M2D , 2/ , K ,M. 5D, 8X , K , 2D" ; 

"DATE :  "&Tday$& "  JULIAN  DA* : " , Julday , 

"LATITUDE:" , Lat," SUBSIDENCE : " ,Wsl , "SURFACE  TEMP: ", 
Sea 

450  PRINT  USING  "/ ,K ,/"; "FORECAST  WINDS  FOR  THE  PERIOD  [ 
KNOTS  J " 

4  55  FOR  Win=»l  TO  Wn 

460  PRINT  USING  " ZZ , 2A, 2X, 3D. D" ? N indin ( W in , 1 ) /2  MOD  24, 
"00",  wind  in  (Win  ,2)  *  1 . 94 
465  NEXT  Win 

470  IF  Hd$=*"Y*  THEN  RETURN 

475  DISP  "PRESS  CONT  WHEN  READ*" 

480  PAUSE 
485  RETURN 

490  Jullat:  PRINT  USING  "@,K,/,K,  3D"  ;  "FOR  DATA  SET:  "S.Loc$& 
"  "SrTimeS  , "CURRENT  JULIAN  DA*:  ",  Julday 
495  INPUT  "ENTER  NEW  JULIAN  DA* ( 1  TO  366),  PRESS  CONT 
" , Julday 

500  PRINT  USING  " 3 , K ,M2D . 4D" ; "CURRENT  LATITUDE  IS  :";Lat 
505  INPUT  "ENTER  LATITUDE (DEG  -90  TO  90),  PRESS  CONT 
",  Lat 

510  RETURN 

515  Enter_winds:  PRINTER  IS  16IKNOTS  =  METERS  *1.94  I 
1  .*94*  ft  *  Nm  *  3600sec 

520  MAT  W indin*ZER  !  SECOND  I 

.3048m  6076.1ft  Hr 
525  ON  ERROR  GOTO  530 

530  PRINT  USING  "@,K,  4X  ,K" ; Loc$ , Time$ 

535  INPUT  "ENTER  THE  BEGINNING  TIME  (  TO  NEAREST  HOUR 
e.g  .  0600) "  ,lnll$ 

540  IF  (LEN  ( InllS)  04 )  OR  (VAL ( I  nl  1$  ( 1 ;  2 J  )  >24 )  THEN  530 
545  IF  Inll$[l?2J= "00 "  THEN  InllS [ 1 ; 2) *"24" 

550  REDIM  G(VAL(InllS [1 ; 2) ) * 2-1 : VAL ( I nl 1$ (1;2) ) *2+61,7) 
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555  OFF  ERROR 
560  Wn*l 

565  Windin(l ,1) =Nr 11»VAL( Inll$ [1 ; 2] )*2  1  PUT  PERIOD  IN 

WIND  HOLDING  ARRAY 
57  0  y$  =  "Y" 

575  Wind  Deginl*Wind  !  USE  THE  VALUE 

FROM  I REPS 

560  PRINT  USING  ■ 3 ,KX , 3D . D, 2XK " ; " WIND  FROM  I REPS  IS  "; 

W ind*l. 94, "KNOTS" 

585  INPUT  "IS  THIS  WIND  OKAY?  PRESS  CONT  TO  USE  THIS 
WIND  (Y/N)",Q$ 

590  IF  Q$[l;  1]  *"Y"  THEN  610 

595  DISP  "ENTER  BEGINNING  WIND  SPEED  [KNOTS]  at  time  "; 
InllS ; 

600  INPUT  Wind  beginl 

605  Wind_beg  inT=Wind  beginl/1.94  ! CHANGE  INTO 

METERS  PER  SECOND 

610  PRINT  USING  " 3 , K , 4A, K , 3D. D, 6A " ; "AT  BEGINNING  TIME:  "; 

I  nil? [1 ,2] i" 00";"  THE  WIND  SPEED  IS  " ; wind_Deginl * 
1.94;"  KNOTS" 

615  PRINT  USING  " 2/ ,K , 2 (/ , K) " ; "Enter  the  nour  and  minutes 
of  eacn  prediction,  then  the  wind. ", "Separate  time 
from  wind  by  a  comma" , "DISPLAYED  winds  in  KNOTS" 
620  PRINT  USING  " / ,K , 2Z  ,16A  ,/ ,K" ; "Tne  last  period  is  30 
nours  after  the  beginning  (i.e.  ",  (Nrll/2+30)  MOO 
24, "00  the  next  day )", "PERIOD  TIME  WIND" 

625  PRINT  USING  "X2D , 5X, 4A  ,5X , 2D .D" ; 1 , InllS [ 1 ; 4 ] , 
wind  beg  in  1*1. 94 

630  G (Nrll+bO , 3) *W indin (1,2) *MAX (Wind_beg ini , .0001) 

635  Q$»"Y" 

640  FOR  Win*2  TO  10  I  MAX  OF  10  WIND  INPUTS .... F IRST 
ONE  FROM  ABOVE 

645  One  wind:  L INPUT  "hhmm,WW  FROM  THE  FORECAST, 

ENTER  TIME,  WIND(  KNOTS)  [END  when  done]",Inl2$ 

650  IF  (Inl2$[l;2]  =  "99")  OR  (UPCS (Ini  2$ [1 ;1] ) ="E")  THEN 
Print  wind 

655  IF  <LEN(Inl2S) <5)  OR  (Inl2$«"")  THEN  645 
660  IF  NOT  ( Inl2$ [ 5 ;  1] =" , " )  OR  (VAL ( Inl2$ [1 ; 2 ] ) >24 )  OR  ( 
Inl2S [ 6 ;1 ] * " " )  THEN  645 

665  PRINT  USING  "X2D , 5X, 4A  , 5X , K " ; W in , Ini  2$ [ 1 ; 4 ]  , I nl2 $ [ 6 ] 

670  Wind  in  (Win,  1)  *VAL(  Inl2S  [1;  2]  )  *2  1  TIMES  2  TO  GET 

i  OF  PERIOD 

6  75  G (N rll+6 0 , 3) *w indin (Win  ,2) =MAX ( VAL ( Inl2S ( 6] ) /l . 94 , 
.0001) 

680  IF  UPCS (0$ [If  1] )  =  "N"  THEN  RETURN 

685  W n*w in  !  COUNT  THE  TOTAL 

WINDS  ENTERED 
690  NEXT  Win 

695  Print  wind:  Q$  =  "Y"  !  HERE  TO  PUT  WIND  IN  'G  '  ARRAY 
700  PRINT  PAGE, "WINDS  THROUGH  FORECAST  PERIOD" 

705  FOR  Win*l  TO  wn 
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710  PRINT  USING  " 5X , 2D , 2X , ZZ , 2A ,  2X, 1D.D" ;W in ,rfind in (Win , 1 ) / 
2  MOD  24, " 00", Wind  in (Win, 2) *1.94 
715  NEXT  Win 
720  0$="Y" 

725  INPUT  "DO  THESE  WINDS  LOOK  OKAY  (Y/N)?”,Q$ 

710  IF  UPC$(Q$[1,H  )*"Y"  THEN  745 
715  IF  UPC$ (Q$(l ,1] ) *"N"  THEN  Fixwind 
740  GOTO  Print  wind 

745  G (Nr 11-1 , lI*G(Nr 11 , 1) =Windin( 1 ,2) 

750  Winl=0 

755  Iaayl=Nr 11 

76  0  FOR  W  in*2  TO  Wn 

765  IF  Windin( Win, 1) <Idayl  THEN  Winl=Ninl+48 

770  G (MIN (Wind  in ( Win, 1 ) +Winl,Nr 11+61) , 1) =W indin( Win, 2 ) 

(TO  PREVENT  TIME  BEYOND 
775  Idayl*Windin(Win,l) 

780  G(Nrll+60,3)«W indin  (Win, 2)  !  AUTOMATIC  PUT  THE  LAST 

WIND  INTO  END  PERIOD 
785  NEXT  Win 
790  RETURN 

795  Fixwind:  INPUT  "ENTER  PERIOD  NUMBER  TO  CHANGE", Win 

800  IF  Win>Wn  THEN  Fixwind 

805  GOSUB  One_wind 

810  GOTO  Print  wind 

815  Subsidence:  “l  COME  HERE  NORMALLY  OR  FOR  RETRY  WITH 
DIFFERENT  SUBSIDENCE 
820  FIXED  4 

325  DISP  "ENTER  SUBSIDENCE  (m/sec)  [Prior  value*" ;Wsl ;"] , 
or  CONT  to  use  same  value" ?CHR$ ( 7) ; 

810  INPUT  Wsl 
815  RETURN 

840  Sea_temp:  PRINT  "PRESENT  VALUE  OF  SEA  SURFACE 
TEMPERATRURE  IS  ";Sea 

845  INPUT  "ENTER  THE  SEA  SURFACE  TEMP  IN  'C  or  just 
CONT  for  same  value"  ,Sea 
350  Tsf c*Sea+271 . 16 

855  RETURN 

860  Assig_error:  IF  ERRN*80  THEN  No  data 
365  DISP  "ERROR  # " ; ERRN ; "  in  line  T ; E RRL ; "  PRESS 

CONT  when  ready"  ,CHR$ ( 7) 

870  PAUSE 
375  GOTO  Start 

380  No  data;  PRINT  PAGE, "DATA  FILE  NOT  AVAILABLE  IN  LEFT 
"TAPE  DRIVE. ",  LIN (2) /’EITHER  NO  TAPE  or  NOT  A  DATA 
TAPE. ", LIN (2) 

885  PRINT  "INSERT  OR  MAKE  OTHER  CHOICE" ,  LIN  (  2)  ,  "PRESS 
CONT  TO  TRY  AGAIN  WITH  DATA  TAPE  IN  LEFT  SIDE 
TAPE  DRIVE" 

890  BEEP 
395  PAUSE 

900  Ireos  data:  !  TO  DISPLAY  THE  FILE  NAMES  OF  DATA  ON 
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I  REPS  DATA  TAPE 

905  ON  ERROR  GOTO  Assig  error 

910  ASSIGN  #5  TO  "ENVIR^StMassS  ,  Ier 

915  IF  Ier  THEN  No  data 

920  ASSIGN  #6  TO  " SEQNCE " &Mass$ , I er 

925  IF  Ier  THEN  No_data 

930  READ  #6;Envsq(*) ,Syssq( *) , Lossq( *) 

935  ON  END  #5  GOTO  All  out 

940  FOR  1*1  TO  16 

945  READ  *5 , I ; EnvnaraS ( I ) 

950  NEXT  I 

955  GOTO  All  full 

960  All_OUt:  ~  FOR  J*I  TO  16  !  I  HAS  SEEN  INCREMENTED  1 

ABOVE  LAST  READ  I 
965  Envsq(J)*0 
970  NEXT  J 

975  All_full:  RETURN  <  ARRAY  FILLED  WITH  NAMES  and 
REMAINDER  WITH  0 's 


980  Disp_data:  PRINTER  IS  3 6  !  TO  PUT  THE  DATA  SETS  UP  FOR 
CHOICE 

985  PRINT  USING  " 3 , K , 5X , K, 2/ , 5X , 2D  ,  2X, K" ;CHR$ ( 7)  , 

"EXISTING  ENVIRONMENTAL  DATA  SETS  ON  IREPS  TAPE:", 

0 , "NEW  or  RETURN  TO  MENU  SELECTION" 

990  FOR  1*1  TO  16 
995  Rec*Envsq(I) 

1000  IF  NOT  Rec  THEN  Pick_data 

1005  PRINT  USING  " 5X, DD,  2X,X";  1  ,Envnam$  (Rec)  [1 ,2  3] 

1010  NEXT  I 
1015  1*17 

1020  Pick  data:  Nrec*I-l 

1025  OFF  ERROR 

1030  IF  Nrec<*0  THEN  Menu 

1035  INPUT  "ENTER  THE  NUMBER  OF  THE  DATA  SET  DESIRED", 
Dummy 

1040  IF  Dummy*0  THEN  Menu 

1045  IF  NOT  ((Dummy>*l)  AND  ( Dummy <=Nrec) )  THEN  Disp_data 
1050  Rec*Oummy 

1055  D ISP  "PROGRAM  IS  WORKING  RECORD  NUMBER  ",Rec,";", 

TRIMS (Envnam$(Envsq (Rec) )[ 1,23] ), "  IS  BEING  READ" 
1060  READ  #5 ,Envsq ( Rec) ;Name$,Loc$, Times , Types ,HigntS, 
EvapS,Wmo$ (* ) , NmohtS 

1065  READ  #5 ; Presur (*) , Temper ( *) ,Relnum( *) ,He ignt (*) , 
Munits(*) , Nunits (*) 

1070  READ  #5;'Wind  , Sea , Air ,  ReliTm,Htzero  ,Przero,Nmax  , Delta , 
Change 

1075  FOR  1*1  TO  Nmax-1 

1080  IF  (Presur ( I)  *0)  OR  (Height(Nmax-I+l) >10000)  OR  (1*2!) 
THEN  GOTO  1110 

1085  Spec  hum ( I ) =FNUva 1 ue ( Temper ( I ) , Relnumf I )  ) 
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10*0  Pot  temp{  I)  “Temper  ( I) +. 0098*Height  (Nmax-I) 

1095  ZsncJ(  I)  *He  ight  {Nmax-I )  *1 

1100  Digit  top*I 

1105  NEXT  I 

1110  Press*Przero 

1115  IF  Press*0  THEN  Press *10] 2  1  DEFAULT  TO  STANDARD  IF 

NO  OTHER  VALUE 
1120  MAT  G*ZER 

1125  Flag  rev*0  3  TO  FLAG  A  NEW  DATA  SET 

1110  Flag~data*l  !  TO  FLAG  A  DATA  SET  CHOSEN 

1135  RETURN  3  GET  THE  WINDS  THEN  DO  VERTICAL 

STRUCTURE 


1140  Pit  struct:  1 

1] 45  PLOTTER  IS  1 3 , " GRAPH  ICS " 

1150  GRAPHICS  i **********laBEL 

VERTICAL  PLOT********** 

12  55  PEN  1 

1160  LINE  TYPE  1 

1165  SCALE  -3,30,-200,2400 

12  70  CLIP  0,30  ,0,2400 

1175  AXES  1,500,0,0,5,2 

1180  UNCLIP 

1185  LORG  2 

2190  MOVE  -3,-6  5 

1195  LA3EL  USING  "K"?"SPEC.  HUM.  Q" 

1200  MOVE  -3,-130 

1205  LABEL  USING  "K";"POT.  TEMP.  "C" 

1210  LORG  8 

1225  MOVE  30,2285  3  ************DRAW  WIND  ON 

PLOT**  ******* 

1220  LABEL  USING  " 8AX" FORECAST " 

1225  MOVE  30,2225 

2230  LABEL  USING  "K"  ; "TIME  WIND" 

1235  FOR  L*1  TO  Wn 
1240  MOVE  30 , 2200-L*  80 
1245  BEEP 

1250  LABEL  USING  " ZZ ,2A, 2X, 3D. D" ; Windin (L ,1 ) /2  MOD  24, "00” 
Windin  (L  ,2)  *1. 94 
1255  NEXT  L 
1260  LORG  5 
1265  FOR  Ll-1  TO  2 
1270  FOR  L* 5  TO  25  STEP  5 
1275  MOVE  L,-65*L1 

1280  LABEL  USING  "M3D" ; L*L1- ( Ll-1 ) *  1 0 

2285  NEXT  L 
1290  NEXT  LI 

1295  FOR  L=500  TO  2000  STEP  500  3  ********labEL  SIDE 

NUMBERS******* 


130  0  .10  VE  -1 ,  L 

no 5  LABEL  USING  "4D";L 

1310  NEXT  L 

13! 5  DEG 

1320  LOIR  90 

1325  LGRG  2 

1330  MOVE  -2.5,0 

1  335  LABEL  USING  “  X" ;  Loc$& ■  “&Tinte$&"  "&"HEIGUT  (meters) 
"&Tday$ 

1340  LDIR  -30 
1345  MOVE  3,2400 

1350  LA3EL  USING  "XK" ; "SPECIFIC  HUMIDITY  0" 

1355  LORG  3 
1360  LDIR  30 

1365  MOVE  22,2400 

1370  LABEL  USING  "X" POTENTIAL  TEMP  0" 

1375  MOVE  22,2400 

1380  LABEL  USING  "A";"-" 

1385  LDIR  0 

1390  LORG  5 

1395  MOVE  25,500  !**********LABEL 

LEGEND********** 

1400  LABEL  USING  "15A";"  V  Q  I  Theta  /" 

1405  LORG  3 

1410  MOVE  30,750 

1415  LABEL  USING  " 1 OA ,2X ,MZ. 4D" ; "S ubs idence" , rtsl 
1420  MOVE  25,600 

1425  LABEL  USING  " 1 3 A" ; "MI XED  LAYER  HEIGdT" 

1430  FOR  L-l  TO  3 

1435  MOVE  25 ,500-L*100 

1440  LABEL  USING  " 14 A, 7 X” ; Stick? (L) 

1445  NEXT  L 

1450  SCALE  -16,50,-200,2400 

1455  MOVE  Pot  temp (1) , Height (Nmax)  !**********DRAW  THETA  * 

*  *  *  *  *  h 

1460  FOR  L»2  TO  MIN  (Nmax- 1 , D Ig  i  t_top) 

1465  DRAW  Pot  temp  (L) , Height  (Nmax-L) 

1470  NEXT  L 

1475  Zi*Zi3 

1480  POINTER  Pot  temp( 2) ,Height(Nmax-3) ,1  !  FULL  SCREEN 

CURSOR  ” 

1485  BEEP  j **********3ET  THETA*** 

***★★★* 

1490  DIGITIZE  Mixtemo,Zi 
1495  MOVE  40,600 

1500  LORG  2 

1505  LABEL  USING  "X4D , 2A" ; Zi j "m . " 

1510  MOVE  40,200 

1515  LABEL  USING  "XM2D. D "  jMixtemo 

1520  BEEP 

1525  POINTER  MixtemD+5 , Z i , 2 
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1510  DIGITIZE  Top  inv  temp,Mmra 
1515  HOVE  40,100  ~ 

1540  LABEL  USING  "X2D.D" ;Top  inv  temp-Mixtemp 

1545  POINTER  Top  inv  temo+1072 35ff,2 

1  550  BEEP 

1555  DIGITIZE  Top  teniD,Top  temo  alt 
1560  MOVE  40,400 

1565  LABEL  USING  "XZ. 4D" ; (Top_temp-Top_inv  temp)/( 
Top  temp  alt-Zi)  1  ""  ’’ 


1570  LINE  TYPE  1  1 *  * *** *****DRAW  STICK  THETA******* 

*  * 

1575  MOVE  Mixtemp, 0 

1580  DRAW  Mixtemp,Zi 

1585  DRAW  Top  inv  temp , Zi 

1590  DRAW  Top” temp, Top  temo_alt 

1595  LINE  TYPE~1 

1600  LORG  8  1 WRITE  UNDER  'Q '  SIDE 

1605  SCALE  -1,10,-200,2400 

1610  BEEP 

1615  LINE  TYPE  8 

1620  MOVE  Spec_hum( 1 ) , Height (Nmax)  i**********draw  Q  *** 
*****★“■ 

1625  FOR  L»2  TO  MIN (Nmax-1 , Dig it_top) 

1610  ORAW  Spec_hum(L) ,  He  ight  (Nmax-L) 

1635  NEXT  L 

1640  LINE  TYPE  1 

1645  POINTER  Spec_hum( 2) ,Zi,l 

1650  DIGITIZE  Mi xq ,Mmm  !  USING  THIS  TERM  SO  SCALE  STAYS 
THE  SAME 

1655  MOVE  25,200 

1660  LABEL  USING  "M2D . 2DX" ;Mixq 

1665  POINTER  Mixq/2 , Zi , 2 

1670  BEEP 

1675  DIGITIZE  Top  inv  q,Mmm  !  USE  FOR  SAME  SCALE  FACTOR 
GRAPHICS 
1680  MOVE  25,100 

1685  LABEL  USING  "M2D . 2DX" ;Top  inv  q-Mixq 
1690  POINTER  Top  inv  q/2, 2350,1 
1695  BEEP 

1700  DIGITIZE  Top  q,Top  q  alt 
1705  MOVE  25,400 

1710  LABEL  USING  "M. 4DX" ; (Top  q-Top  inv  q)/(Top  q  alt-Zi) 
1715  Zi3-Zi 

1720  LINE  TYPE  4  ! *****#****DRAW  STICK  0  ****** 

#** 

1725  MOVE  Mixq  ,0 

1730  DRAW  Mixq,Zi 

1735  DRAW  Top_inv_q,Zi 


1740 

1745 

1750 

1755 

1760 

1765 

1770 

1775 

1780 

1785 

1790 

1795 

1800 

1805 

1810 

1815 

1820 

1825 

law 

1835 

1840 

1845 

1850 


1855 

1860 

1865 


1870 

1875 

1880 

1885 

1890 

1895 

1900 

1905 

1910 

1915 

1920 

1925 

1930 

1935 

1940 

1945 

1950 


DRAW  Top  q,Top  q  alt 
LINE  TYPE  1 

SCALE  -16,50,-200,2400 
MOVE  40,100 

LORG  5  i **********  FOR  DECISION  OF  A  GOOD  FIT  OR 

NOT********** 

LABEL  USING  "17A";"  OKAY  <-|->  AGAIN" 

POINTER  40,100,2 
BEEP 

DIGITIZE  D,Mmn> 

IF  D*40  THEN  1770 

IF  D >40  THEN  Pit  struct 

PRINTER  IS  16 

PRINT  PAGE 

MOVE  40,100 

DISP  "CLEARING  PICTURE" 

EXIT  GRAPHICS 
PEN  -1 

LABEL  USING  "17A"?"  OKAY  <-|->  AGAIN" 

GOSUB  Hard  output 
IF  HdS*"N"_THEN  1885 
PRINTER  IS  0 

PRINT  USING  "  1X,K,2(6X,K) ,/"; "ENVIRON.  DATA  LIST, 

FROM  MIXED2" ,Loc$,Time$ 

PRINT  USING  "3(10A,X4D.D,9X18A,M3D.D,/) ", •"WIND", Wind, 
"AIR  TEMP" , Air , "SEA  TEMP" , Sea, "REL  HUM* ,Relhm, "SFC 
PRESS" ,Przero,"IREPS  EVAP  DUCT  HT", Delta 
PRINT  "LEVEL  PRESS  TEMP  RH  %  METERS  M  UNITS 
POTENT  TEMP  SPEC  HUM" 

FOR  1*1  TO  MIN (18, Digit  top) 

PRINT  USING  "3D,3X,4D.D72X,M2D.D,X,3D.D,X,6D.D,4X,4D,2 ( 
6X,  M3D.D)  "  ;  I  , Presur  ( I)  ,Temper  ( I )  ,Relhum(  I)  ,  He  ight  ( 
Nmax-I )  ,Munits(Nmax-l )  ,Pot  temp  (I), Spec  num(I) 

NEXT  I 

DUMP  GRAPHICS 
PRINT  USING  "<§" 

PRINTER  IS  16 
EXIT  GRAPHICS 
RETURN 

Run:  Tl0=Mixterap  i *** *******prep  DATA  FOR  RUN  WITH 

IREPS********** 

Zi=Zi3 
Qp*M  ixq 

Dqp*Top_inv_q-Mi xq 

Dqd z  u*  ( T op_q  -Top_  in  v_q )  /  { To p_q_ a  1 1-  Z  i ) 

Dth*Top_inv_temp -Mix temp 

D td z u* ( Top_ temp -To p_i n  v  te  mp )  /  ( T  op_  te rap_a 1 1- Z i ) 

Ns*WSl  ~ 

Frame__loc$*Loc$ 

Frarae”date$*Time$ 

EXIT  GRAPHICS 
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/ 


1955 


1960 

1965 

1970 

1975 

1980 

1985 

1990 


1995 

2000 

2005 

2010 

2015 

2020 

2025 

2030 

2035 

2040 

2045 

2050 

2055 

2060 

2065 

2070 

2075 

2080 

2085 

2090 

2095 

2100 

2105 

2110 

2115 

2120 

2125 

2130 


CALL  3 ky2 ( Zi , Umax* l , Pot_temp ( *) , Zsnd (*) , Spec_hum( *) , 
Tsky,Ftop) 

Mhr*6  0  i  TO  SET  >1HR*  NUMBER  Of  HALF  HOUR 

PERIODS 
OFF  ERROR 

PRINT  PAGE , TAB (9) , "THIS  IS  'MIXED2',  MODELING  THE 
CHANGE  IN  THE  MIXED  LAYER" 

PRINT  TAB(IO) , "ABOVE  THE  OCEAN  SURFACE.  Today  is  "& 
Tday$ 

Numhr*Nrll+Mhr 
PRINTER  IS  16 

PRINT  USING  "/,15X,K,/,20X, 4Z , 2XK ,/19X,4D,K,/"; "DATA 
PARAMETERS  CHOSEN" , Nr  11/2  MOD  2 4* 1 00  , "BEGINNING 
TIME" , Zi, "m  CAPPING  INVERSION  HEIGHT" 

Vd  t*0 

!  INPUT  "ENTER  A  VALUE  FOR  'Vdte'  (  'C/12hr)  ,  or  PRESS 

CONT  FOR  DEFAULT  TO  0  ?",Vdt 
Vd  te*Vd  t/2  4 

!  PRINT  "THE  VALUE  FOR  'Vdte  '  IS  "  ; Vd t;  "  'C/1 2hr  or 
Vdte;"  'C  per  TIME  STEP" 

Th*Tl0+273 .16 
JaNrll-1 
G (Nrll-1 , 1) -Zi 
GOSUB  Bulk 
Nrl-Nrll 

I*Nrl  1  FIND  FIRST  INPUT  DATA  LINE 

1*1+1  1  INCREMENT  TO  NEXT  LINE 

IF  G(I,3)*0  THEN  20451  IF  NO  DATA,  INCREMENT 
Nd-I-Nrl  i  FIGURE  TIME  SPREAD 

X5»(G(I,3)-G(Nrl,3))  /Nd  !U10 
FOR  J*1  TO  Nd-I  !  STEP  THROUGH  EMPTY  POINTS 

G(J+Nrl,3)*G(Nrl,3)+X5*J 
NEXT  J 

gradient  above  inv. 

Nrl*I 

temperature  jump 
IF  I <N umhr  THEN  204  5 

gradient  above  inversion 
T*Th-.0098*  Zi*. 5-273 .16 
humidity  jump 
LCD*  (5  96. 73-.  601*'.  )/.  24 
ootential  temperature 
T  t*T- . 009 8* Zi*. 5  +  273. 16 
soecific  humidity 
Te*Th  *  ( 1+Lcp/Tt  *Qp/1000) 

D  te*D  th+Th/Tt  *Lcp*  Dqp /l 00  0 
between  solutions 
Dtedzu*Dtdzu+Th/Tt*Lcp*  Dqdzu/1000 
Dqw*Dqp/1000 
Qw»Qp/1000 

Aws*W s/Zi  1  AVERAGE  NS  THRU  THE  MIXED 


!  DQDZU*  Specific  hum 

!  DTH=  Potential 

!  DTDZU*  Temp 

!  DQW*  Specific 

J  TH*  Well  mixed 

!  Qp*  Well  mixed 

!  ZI*  Inversion  heignt 
!  DELTIM=Time  step 


91 


D-  11 


LAYER 

2135  Zws*Zi  !  REMEMBER  HEIGHT  OF 

BEG  INN  I  MG  Ns 

2140  PLOTTER  13  1 3 , "GRAPH ICS" 

2145  SCALE  Nr  1 1-20 , Nr  11  +  80 , -100 , 2000 
2150  CLIP  Nrll , Nr  11+61, 0,1200 
2155  AXES  4, 200  , Nrll,  0,3, 5 
2160  FRAME 
2165  GRAPHICS  1 


2170  RunS3"* " 

2175  Looper:  FOR  J=Nrll  TO  Nrll+Mhr 

2180  Zil*Zi  l  HOLD  ONTO  PREVIOUS  HEIGHT 

2185  PLOT  J,Zi  !  FOR  THE  QUICK  PLOT 

2190  T10=Th-273 . 16  !  T10  CARRIED  BACK  FROM  'BulK '  AS  Th 

2195  Q10*Qp  i  NEW  VALUE  OF  Qp 

2200  FIXED  2 

220  5  D  tedzu*Dtdzu+Th/Tt  *Lcp*Dqdzu/10Q0 

2210  IF  J>Nr  11  +  56  THEN  RunS-T* 

2215  DISP  "PERIOD  BEGIN  :  "  ; TAB  ( 14+  ( J-Nr  11 )  /2 )  ;  Run$  ;  TAB  ( 

45);": END  Layer  at" ;G( J-l ,1) 

2220  Q2-Dqdzu/1Q00  !  TO  CHANGE  G/KG  TO  G/ 

G 

2225  Timesec»J* 1800 

2230  CALL  Zenith  (Tijuesec ,Julday,Lat, Theta) 

2235  CALL  M2 (T10,Q10, Zlcl , Zi , Ustar , Tstar ,Qstar ,Tstv ,ws, Delr r 
02,Dte  ,Dqw  f  Dhdt  ,Ddtedt  ,Ddqwdt,  Dtedt, Dqwdt,  Wtv , 
Dtedzu,Tsfc,Theta,Tsky) 

2240  G(J,2)=Zlcl 

2245  G(  J,  1)  *Zi»Zi+Dhdt*Deltim 

2250  Qw*Qw+Dqwdt*De  lv.im  !  GRAMS  PER  GRAM 

2255  G(J,  4)=‘Q10=Qp=Qv;i1000  i  GRAMS  PER  KILOGRAM 

2260  Te*Te+Dtedt*Deltim  l  EQUIVALENT  POTENTIAL 

TEMP 

2265  Dqw*Dqw+Ddqwdt*De ltim  !  GRAMS  PER  GRAM 

2270  G  ( J  /  6)  »Dqp*Dqw*  1000  !  GRAMS  PER  KILOGRAM 

2275  Dte*Dte+Dd tedt*Deltim+Vdte  l  ADDING  THERMAL  WIND 

2280  Ws=Aws*Zi  !  RESETTING  LAYER 

SUBSIDENCE 

2285  Zws»Zws+Aws* Zws*Deltim 

2290  Tt-Tt-,0098* (Zi-Zil)  I  Tt  IN  KELVIN 

2295  T*Tt+. 0098*Zi*.5-273.2 

2300  Lcp* (596. 73-. 601 *T)/. 24 

2305  G( J,5) »Th»Te/{l+Lcp*Qw/Tt)  !  Th  in  KELVIN 
2310  G  ( J,  7)  »Dth*Dte-Th/Tt*Lcp*Dqw 

2315  !  G(  J ,2)  *MIN ( 100 ,  FNRelhunt  ( Qp, Th )  )  !  FIGURED  AT  10 

METER  LEVEL 
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2320 


2325 

2330 

2335 


2  340 

2345 

2350 

2355 

2360 

2365 

2  370 

2375 

2380 

2385 

2390 

2395 

2400 

2405 

2410 

2415 

2420 

2425 

2430 

2435 

2440 

2445 

2450 

2455 

2460 

2465 

2470 

2475 

2480 

2485 

2490 

2495 

2500 

2505 


JOS U3  Bulk  J  COMPUTE  STAR  VALUES  FOR  NEXT  TIME 
PERIOD 

NEXT  J 

EXIT  GRAPHICS 
GOTO  Menu  ! 


Bulk:  !  TO  COMPUTE  THE  NEXT  SET  OF  VALUES  FOR  U*,T*,Q*, 
T*v 

U10-MAX ( .2,G( J, 3) )  J  WIND  CHANGES  THROUGH  THE  PERIOD 
IF  U10>*2 . 2  THEN  Cdn*.  789*U10*  .  259/1000  !  FIGURE 
DRAG  FROM  WIND  CURVE  FIT 

IF  U10<2. 2  THEN  Cdn=U 1 0" (- . 1 5) *1 . 08/1 000  ![  DAB 
SEP  80] 

RAD 

Tdelta*Th-Tsfc  i  DIFFERNCE  IN  TEMPERATURE  AND  Q- 
value 

Qdelta*Qp~FNQvalue (Sea, 98)  !  FROM  lOmeters  TO  THE 
SURFACE 

Zo*Z*EXP (-Kkk/SQR (Cdn) ) 

Ctn_sqrt*Alpha  t*Kkk/LQG  { Z/Zot) 

Sl*S*So*Kkk*9 . 8*Z/Th*  (Ctn_sqr  t/Cdn)  *  (Tdelta+ .00061* 
Th*Qdelta) /U10~2 
IF  ScKO  THEN  24  50 
IF  So<  2  THEN  GOTO  2420 
S=10 

Psil*-47 
Psi2*-64 
GOTO  Stars 
Psil—SM.7 
Psi2*-S*6 . 4 

3»So*Alpha_t/Kkk*  ( (Kkk-SQR(Cdn)  *Psil)  ~2/(Alpha_t* 
Kkk-Ctn”sqrt*Psi2) ) 

IF  ABS(S-Sl) <7001  THEN  Stars 
Sl*S 

GOTO  2420 

XI* (1-15*S) ".25 
X2* ( 1- 9*S ) * . 5 
Psi2*2*LOG ( ( 1+X2) /2) 

PS  il*2*LOG ( (1+Xl) /2) +LOG( (1+X1~2) /2) -2*ATN (XI ) + 
PI/2 

S*So*Alpha_t/Kkk* ( ( Kkk-SQR (Cdn) *  Psil ) "2/ (Alpna_t* 
Kkk-Ctn”sqrt*  Psi2) ) 

IF  ABS(S-Sl) <7001  THEN  Stars 
S1*S 

GOTO  2450 

Stars:  Mult*Alpna  t* Kkk/ ( LOG ( Z/Zot) -Ps i2 ) 

Ustar*Kkk*U10/ (LOG (  Z/Zo)  -Psil) 
Tstar*Tdelta*Mult 
Qstar*Cdelta*Mult 
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25 TO  Ts tv=Tstar+. 61*Qstar*Th/l  000 

2515  RETURN  i  VALUES  BACK  FOR  BEGINNING  AND  FOP.  DURING 
RUN 


2520  Picture:  EXIT  GRAPHICS 
2525  IF  NOT  Flag  data  THEN  Menu 
2530  PRINTER  IS  16 

2535  OFF  ERROR 

2540  PRINT  USING  “ ? ,K , 2/ , 3 (1  OX, X , 2/) " ; "TYPE  OF  PLOT:'V'l 

PLOT  THE  'BIG  PICTURE 2  PLOT  OF  'o',  'THETA' 

AND  JUMPS", "3  LIST  OF  THE  PLOTTED  VALUES" 

2545  PRINT  "  4  RETURN  TO  DIRECTION  MENU" 

2550  ON  ERROR  GOTO  Picture 

2555  INPUT  "ENTER  MODE  of  output?", Vss 

2560  PRINT  PAGE, "PRESS  CO NT  wnen  picture  is  done", LIN ( 
2),"  (  k 3  TO  DISPLAY  GRAPHICS)" 

2565  IF  Vss*l  THEN  CALL  Fr ame (G ( * )  , Nr  1 1 , TdayS , Dqd zu , 
Dtdzu, Press , Frame_loc$ , Frame_date$ ,Mhr ) 

2570  ON  Vss  GOTO  Picture ,Mix_plot ,Gval , Menu 

2575  £  ON  Vss  GOTO  Picture , Picture ,3 val , Menu 
2530  GOTO  Picture 

2585  Hard  output:  Hd$="N" 

2590  INPUT  "DO  YOU  WANT  HARD  COPY  OUTPUT  (N/Y)",HdS 
2595  Hd$*UPC$ (Hd${  1  ,1]  ) 

2600  IF  (Hd$= "Y")  OR  (Hd$="N")  THEN  RETURN 

2605  GOTO  Hard  output  £ 


2610  Gval : !  PRINT  OUT  THE  VALUES  FOR  THE  ENTIRE  PERIOD 
2615  GOSUB  Hard  output 
2520  OFF  ERROR 

2625  IF  Hd$="Y"  THEN  PRINTER  IS  0 

2630  LS="This  is  a  list  of  plotted  values.  "iLocSSt"  "& 
TimeSS"  "&Tday$ 

2635  PRINT  L$[!;78) ,LIN(3) 

2640  PRINT  "HOUR  Zi  Zlcl  RelHum  WIND  Th  Qp 

Thv  Dtn  Dtnv  Dqp" 

264  5  FOR  Gl-drll  TO  Nrll+Mhr 

2650  G2=INT(G1  MOD  43/2) *1 00+FRACT (G1  MOD  48/2) *60 

2655  33*3(32 ,5) * ( !+. 0006 !*G (G1 , 4) ) i  Thv 

2560  G4=G(G1  ,7) +.00061*3 (Gl ,5) *G  (  Gl ,6)  !  Dthv 

2565  G5  =  FNRelnum((G(GI  ,4))  ,  (G(G!  ,5))  ) 

2670  PRINT  USING  " 4 Z 2X, 2 ( 4DX) , 2 (M 3D . DX) , 6 ( 3D. D2X) " : G2 , G ( 
32 ,1) ,G(G1 ,2)  ,35,3  (Gl  ,3)  ,G(G!  ,5) -273.1  6,G(G1,4)  , 
-j 3-273. 1o,G(j1  ,7)  ,G4,G(Gl  ,5) 

2575  NEXT  Gl 

2530  PRINT  "HOUR  Zi  Zlcl  RelHum  WIND  Th  Co 

Thv  Dth  Dtnv  Dap" 

2535  0 ISP  "PAUSING  FOR  2  SECONDS  THEN  RETURNING  TO  THE 

START” 
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2690  IF  Hd$="N"  THEN  WAIT  2000 

2695  !  WHEN  YOU  ARE  READY ,  JUST  PRESS  CONT 

2700  PRINT  PAGE 

2705  GOTO  Picture! 


2710  Mix  clot:!  BEGIN  PLOT  OF  Qd, Dqo, Th, Dth ,RH 

2  715  OFF  ERROR 

2720  PLOTTER  IS  13 , "GRAPHICS" 

2725  GPAPHICS 
2730  FRAME 

2735  SCALE  Nr ! 1-5 , Nr  1 1+6 0,-8 , 3 
2740  CLIP  Mr  17 ,Nr 1 1 +60  ,-8 , 3 
2  74  5  AXES  2 , 2 ,M r 1 1 , 0 , 3 , 1 

2750  LORG  2 

2755  MOVE  Nr  11, 7. 5 

2760  LABEL  USING  " K" ; "  POTENTIAL  TEMPERATURE  [  'C ] " 

2765  MOVE  Nr!!  ,7 

2770  LABEL  USING  "K";"  POT  TEMP  jump  -  - 
2775  MOVE  Nr!!, -.5 

2780  LABEL  USING  "K" SPECIFIC  HUMIDITY  [ g/kg ] " 

2785  MOVE  Nr! 1,-1 

2790  LABEL  USING  " K" ; "  SPEC  HUM  jump  -  - 
2795  FOR  K=0  TO  2  STEP  2 
2800  LORG  3 

2805  SCALE  Nr  11-5 ,Nr !! +60 ,INT(G (Nr !! ,4+K) ) -4 , INT (G (Nr  1 1  , 
4+K) ) +12 

2810  FOR  1=0  TO  ! 

2815  P—2 

2320  LINE  TYPE  1 

2825  FOR  L*  INT  (G  (Mr  1!  ,  4+K+I)  -3-273  .  IS*  ((  K=0)  AND  (1=1))) 

TO  INT(G (Nr !!, 4+K+I) +3-273. 16* ((K=0)  AND  (1=1))) 
28  30  MOVE  Nr  1 1 +K*  30  ,L 

28  35  LABEL  USING  ,,MDDX";L 

2840  NEXT  L 

2845  LINE  TYPE  !+3*(K=2) 

2350  FOR  J=Nr 1!  TO  Nr  11+60 

2355  PLOT  J,G (J, 4  +  I+K)-( ( 1  =  1)  AND  ( K= 0)  ) *  273 . 1 6  ,  ? 

2360  P=-l 

2865  NEXT  J 

28  70  SCALE  Nr ! 1-5 ,Mr 11+60 , INT (G (Nr ! 1 , 5+X) -!  2-2 73 . ! 6* (K=0) 

) , INT (G (Nr  11, 5+K) +4-273. 16* (X=0)  ) 

2875  NEXT  I 

2380  LINE  TYPE  4 

2885  NEXT  K 

2390  G0SU3  Hard  cutout 

2395  IF  Hd$  =  "YM”*THEN  PRINTER  IS  0 

2900  PRINT  USING  "X" ; "TEMPERATURE  and  HUMIDITY  for  "&Loc$& 
"  " &  T ime  3 

2905  IF  Hd$="Y"  THEN  DUMP  GRAPHICS 
2910  GCLEAR 
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29!  5 
2920 
2925 
2930 
2935 
2940 
2945 
2950 
2955 
2960 
2965 
2970 

2975 
2980 
2985 
2990 
2995 
3000 
3005 
3010 
3015 
3020 
3025 
30  30 
3035 
3040 
3045 
30  50 


30  55 
3060 


3065 

3070 
3075 
3080 
3085 
3090 
3095 
3100 
3105 
313  0 


SCALE  Nr!!-5, Nr  11+6  0,-20,120 
LINE  TYPE  1 

CLIP  Nr!l-5,Nr!l+60, 40,100 
FRAME 

AXES  2, 10, Nr  11 ,40 , 3,2 

MOVE  Sri! , FNRelhum{ (G(Sr!l,4) ) , (G (Nr 2 ! ,5) ) ) 

FOR  I*Nr  2 1  +  1  TO  Nr  11 +60 

DRAW  I , FNRelhum( (G ( I , 4) ) , (G ( I ,5)) ) 

NEXT  I 
LORG  2 

MOVE  Nr  11 ,105 

LABEL  USING  "K";"  RELATIVE  HUMIDITY  at  10  meter 
level" 

LORG  7 

MOVE  Nr  11, 80 

LABEL  USING  " 3 AX" ; " 30% " 

MOVE  Nr  12, 6  0 

LABEL  USING  "3AX";“60%" 

FOR  Cc=Nrl!  TO  Nr  2 1  +  6  0  STEP  12 
MOVE  Cc,110 

LABEL  USING  "ZZ ,ZZ" ; INT (Cc/2)  MOD  24 ,FRACT (Cc/2) *60 
NEXT  Cc 

IF  Hd$»"Y"  THEN  DUMP  GRAPHICS  40,115 
BEEP 

IF  Hd$  =  "N"  THEN  PAUSE 
EXIT  GRAPHICS 
IF  Hd$*"N"  THEN  Picture 
PRINTER  IS  0 

PRINT  USING  "2 (/, 9XX) "30-HOUR  PREDICTION  OF 

POTENTIAL  TEMPERATURE  AND  HUMIDITY" , "The  well- 
mixed  value  scale  is  along  the  left,  the  'jump' 
value" 

PRINT  USING  "8XK , ; "scale  is  along  the  right.  The 
third  frame  is  Relative  Humidity." 

GOTO  Picture  1 


Bye:  PRINT  PAGE , "THIS  PROGRAM  IS  ENDING  AND  REWINDING 

THE  TAPE" 

PRINT  LIN ( 1)  , "REMOVE  TAPE  BY  PRESSING  THE  EJECT  BAR" 
ON  ERROR  GOTO  3085 
REWIND  " :T15" 

ON  ERROR  GOTO  3095 
REWIND  " :T14" 

PRINT  LIN(6) ,TAB( 36) ,"Good  bye" 

BEEP 

OFF  ERROR 
En:  END! 
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3135  SUB  M2 (T3 0 ,Q10 ,Zlcl , Zi ,Ust ,Tst ,Qst ,Tstv, Ws, Delr , Dqdzu , 
Dte, Dqw  ,Dhdt,  Ddtedt  ,Ddqwdt,  Dtedt,  Dqwdt  ,Wtv,  Dtedzu, 
Tsfc, Theta, Tsky) 

3120  Tt-TlO+273. 2-,0098*Zi/2 

3125  Thet»(Tt+Zi*.0098)/Tt 

3130  Cl»<596.73-.601*(Tt-273.2))/(.24*Tt)  !  Cl*L/(Cp*T) 

3135  Wt— Tst*Ust 

3140  Wtv»-Tstv*Ust 

3145  Wq*-Qst*Us  t/1000 

3150  Wte-Wt+Thet*Cl*Tt*Wq 

3155  Te«T10+2.43*Ql0 

3160  CALL  ParcelO (Te ,Q10 , Zi ,Rh ,Thr ,Qv,Ql , Zlcl) 

3165  Qs»Qv/1000*Rh/100 

3170  Qliq*Ql/1000 
3175  !  WAIT  1000 

3180  CALL  Ne2 ( Zi, Zlcl, Dte ,Dqw ,Qs ,Thr ,Qliq,T 10 ,Tsf c,Ust,Tst, 
Qst ,Tstv ,We2 ,Delr , Tsky)  !  NEED  TO  CALL  PARCEL 

FIRST 

3185  CALL  Shrtwv(Theta,Ql,Zlcl,Zi,QlO,Dtemp) 

3190  We»We2 

3195  Dhdt»We+Ws  1  Rate  of  cnange  of  mixed  layer  deptn, 
Ws  is  subsidence 

3200  Ddtedt=Dtedzu*We+ (-We*Dte-Wte+Thet*Delr) /Zi-Dtemp 

3205  Ddqwdt=Dqdzu*We+ (-We*Dqw-Wq) /Zi 

3210  Dtedt* (Wte-Thet*Delr+We*Dte) /Zi+Dtemp 

3215  Dqwdt* (Wq+We*Dqw) /Zi 

3220  SUBEND  !  END  OF  'M2'  SUBPROGRAM 


3225  DEF  FNQvalue ( SHORT  Tx,R)  1  FIGURE  THE  '0'  VALUE 
GIVEN: 

3230  Qvb».053542*EXP(-5399.286*(l/(Tx+273.1o) -3. 5905E-3) ) *R 

3235  RETURN  Qvb  1TAKE  THE  Q  value  BACK 
3240  FNEND1 


3245  Frame:SUB  Frame(SHORT  G ( *) ,REAL  Nrl 1 ,Tday$ , Dqdzu, Dtdzu  , 
Press, Frame  loc$, Frame  dates, Mhr) 

3250  OPTION  BASE  l" 

3255  DIM  Zg (4) 

3260  PLOTTER  IS  13 ," GRAPHICS" 

3265  GRAPHICS 

3270  SCALE  Nr  1 1-6 ,Nr 11+6 0 , -20 0 , 2 400 

3275  CLIP  Nrll  ,Nrll +60, 0,2400 

3280  AXES  2, 100, Nrll, 0,6, 2, 3 

3235  MOVE  Nrll+26,-110 

3290  LORG  5 

3295  DEG 

3300  LABEL  USING  "X" ; "A tmospher ic  Mixed  Layer  Movement  30 
Hours  "&TdayS 
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3305  MOVE  Nrll+26,-170 

3310  LABEL  USING  "K"; Frame  loc$&"  "&Frame  date$&"  LT  0. 

BROWER  NPS" 

3315  LDIR  90 

3  320  MOVE  Nr  11-4. 5, 600 

3325  LABEL  USING  "K";"Mixed  Layer  ( meter  s)  ■ 

3330  MOVE  Nrll-4.5, 1800 

3335  LABEL  USING  "X";"DUCT  HEIGHT  (meters)" 

3340  CSIZE  2.7, .6 
3345  LDIR  0 

3350  MOVE  Nrll+30 , 1050 

3355  LABEL  USING  "K" ; "Ver tical  Profile  \M-value'  (every 
SIX  nours)  " 

3360  FOR  Cc*Nrll  TO  Nrll+Mhr  STEP  12 

3365  MOVE  Cc,-30  (Print  out  time  at  bottom 

3370  Bb*INT(Cc/2) 

3375  Bbb*FRACT(Cc/2) *  60 

3380  Bb*8b  MOD  24 

3335  LABEL  USING  " ZZ , ZZ" ; Bb, 3Qb 

3390  Bb=(3b+3)  MOD  24 

3395  MOVE  Cc+6,2370  i  PRINT  OUT  TIME  AT  TOP 
3400  LABEL  USING  "ZZ, ZZ" ;Bb,3bb 
3405  NEXT  Cc 

3410  FOR  Cv»l  TO  5  !  Label  side 

3415  MOVE  Nrll-2 ,Cv*200 
3420  LABEL  USING  "4D";Cv*20Q 
3425  NEXT  Cv 
3430  MOVE  Nr  11  ,G  (Nr  11 , 1) 

3435  GRAPHICS 

3440  FOR  V=Nrll+l  TO  Nrll+Mnr!  Draw  out  the  max  'a'  height 
at  the  inversion 
3445  DRAW  V,G(V,1) 

3450  NEXT  V 

3455  LINE  TYPE  4 

3460  MOVE  Nrll  ,G(Nrll  ,1)  *1.1 

3465  FOR  V»Nrll+l  TO  Nr 11+Mhr i Draw  out  the  min  'M'  neignt 
above  the  inversion 
3470  DRAW  V,G(V,1)*1.1 

3475  NEXT  V 

3480  LINE  TYPE  2 

3485  MOVE  Nrll  ,G(Nrll  ,2) 

3490  FOR  V*Nrll+l  TO  Nrll+Mnr!  DRAW  THE  LIFTING 
CONDENSATION  LEVEL 
3495  DRAW  V,G(V,2) 

3500  NEXT  V 

3505  LINE  TYPE  1! 


3510  Profile:!  COMPUTES  AND  PLOTS  PREDICTED  PROFILES 
3515  SCALE  -1.5,15,-1400,1200  !  Draw  vertical  lines  in 

lower  window 
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3520  LINE  TYPE  3 

3525  FOR  0-1.5  TO  13.5  STEP  3 

3530  MOVE  0,-1200 

3535  DRAW  0,0 

3540  NEXT  D 

3545  LINE  TYPE  1  !Get  set  for  upper  window 

3550  CLIP  0,15,-1200,1200 
3555  AXES  1,100,15,0,3,2,7 
3560  GRAPHICS 

35o5  Sbd-0  !  FLAG  FOR  SURFACE  BASED  QUCT  WARNING  TO  PRINT 
OUT 

3570  CSIZE  2. 7, .6 

3575  FOR  Cv-1  TO  5  i  label  left  side 

3 5d0  MOVE  -, 5,Cv*  200 

3505  LABEL  USING  " 40" ;Cv* 200 

3590  NEXT  Cv 

3595  Zg(l) -0 

3600  Zg(4 )*1200  ! 


3605  FOR  Pr-0  TO  4 
3610  Pm-Pr*150 

3615  Hrg-Nrll+6+12*Pr  i  Pull  out  nours  for  profiles 

3620  FOR  1-1  TO  4  !  Determine  multiples  to  use  for 

altitude 

3625  Zg(2)*G(Hrg  ,1) 

3630  Zg(3)»G(Hrg ,1) *1.1  !  TO  SET  INVERSION  THICKNESS  -.1* 

HEIGHT 

3635  !  COMPUTES  Q  AND  THETA 

3640  Ppl-G (Hrg ,4) + ( I>-3 ) *G (Hrg , 6) + ( I- 4) *Dqdzu* (1200-Zg (3) ) 
3645  Pp2-G (Hrg ,5 ) + ( I>*3 ) *G(Hrg , 7) + ( 1-4 ) *0 tdzu* ( 1200-Zg ( 3 ) ) 
3650  !  COMPUTES  PRESSURE  , VAPOR  PRESSURE  AND  TEMP ( K) 

3655  P«Press-. 120*Zg ( I) 

3660  E-P*  Ppl/622 
3665  T-Pp2- .0098*  Zg ( I) 

3670  I  COMPUTES  MODIFIED  INDEX  OF  REFRACTION 
3675  Pp(I) -77.6*P/T+373000*E/(T*T)+.157*Zg(I) 

3680  NEXT  I 

3685  LORG  5  1  Label  bottom  of  'm'  profile 
3690  SCALE  225,1050,-1400,1200 
3695  MOVE  Pm+375,-30 
3700  LABEL  USING  "K";"350  400" 

3705  MOVE  Pp{l) +Pm,Zg (1)  !  Draw  eacn  profile 

3710  FOR  1-2  TO  4 

3715  DRAW  Pp( I)+Pm,Zg ( I) 

3720  NEXT  I 

3725  IF  (Pp( 3) <Pp(l) )  OR  (ABS (Pp ( 3 ) -Pp (1 ) ) <5)  THEN  GOSU3 
Sfc  duct 
3730  LINE  TYPE  3 
3735  MOVE  Pp ( 3 ) +  Pm ,Zg ( 3 ) 

3740  DRAW  Pp  (  3)  +Pm, MAX  ( 0  ,  (Pp(  3) -Pp  (1 )  )  /  (  Pp  (2  ) -P?  (1 )  )  *Zg  (  2 ) 
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!  MOVE  TO  EACH  SIX 


3745 

3750 

3755 

37b0 

3765 

3770 

3775 

3780 

3785 

3790 

3795 

3800 

3805 

3810 

3815 

3820 

3825 

3830 

3835 

3840 

3845 

3850 

3855 

3860 

3865 

3870 

3875 

3380 

3885 

3890 

3895 

3900 


) 

LINE  TYPE  1 
MOVE  3  75  +  Pm,-l 1 50 
HOUR  SPOT 
A*FNRelhum((G(Hrg  ,4)  )  ,  (3(Hrg ,5) )  ) 

LABEL  USING  "3D.DM;A 

NEXT  Pr 

SEEP 

PAUSE 

EXIT  GRAPHICS 
Hd$*"N" 

INPUT  “DO  YOU  WANT  HARD  COPY  OUTPUT  (N/Y) ” fHd$ 

IF  Hd${l,lJ*"N"  THEN  SUSEXIT 
PRINTER  IS  0 
PRINT  LIN( 2) 

DUMP  GRAPHICS 
IF  NOT  Sod  THEN  3840 

PRINT  LIN  ( 2)  , "  The  possibility  of  a  surface 

based  duct(SSD)  is  indicated  and  only" 

PRINT  "due  to  tne  'A '  value  greater  at  the  surface 
than  at  the  inversion" 

PRINT  "height.  Variation  in  tne  vertical  structure 
must  be  empnasized" 

PRINT  "when  briefing  this  forecast  in  that  a  SBD  may 
not  exist." 

PRINT  LIN  ( 1)  , "  This  is  a  SIMPLIFICATION  of  the  real 
structure. " 

PRINT  LIN (1) ,  "  This  display  is  divided  into  UPPER 
and  LOWER  windows:" 

PRINT  "LOWER  WINDOW  displays  the  top  and  middle 
of  the  elevated  layer  and  its" 

PRINT  "  forecast  continued  movement  for  a  thirty 
hour  period  after  beginning." 

PRINT  "  At  the  bottom  is  the  Relative  Humidity  at 
each  six  hour  period." 

PRINT  "  The  lightly  dotted  line  is  the  lifting 
condensation  level." 

PRINT  "  {With  enough  moisture  near  the  inversion, 

this  can  De  used" 

PRINT  "  as  a  flag  for  possible  cloud  formation, 

or  for  nign  mixed-" 

PRINT  "  layer  humidity,  near  tne  surface,  to 

forecast  fog.]" 

PRINT  "UPPER  WINDOW  picks  out  M-value  structure  ( 
using  only  4  points) " 

PRINT  "  at  each  six  nour  period  and  will  indicate  a 
surface  based  duct" 

PRINT  "  only  if  the  elevated  M-value  is  less  tnan 
or  witnin  5  of" 

PRINT  "  the  surface  M-value.  The  sampled  times  are 
displayed  at  top." 
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3905 
3910 
3915 
3  920 
3925 
3930 
3935 
3940 
3945 
3950 
3955 
3960 


3965 

3970 

3975 

3980 

3985 

3990 

3995 


4000 


4005 

4010 

4015 

4020 

4025 

4030 

4035 

4040 

4045 

4050 

4055 

4060 

4065 

4070 

4075 

4080 

4085 

4090 

4095 


PRINT  PAGE 
SUOEX IT 

Sfc  duct:  LORG  1 
beIp 

LDIR  90 

MOVE  Pp  (  3)  +Pm-5 ,  Zg  (  3 )  +  20 
LABEL  USING  "K " ? "Surf ace  Duct" 
Sbd®  1 


LDIR  0 
LORG  5 
RETURN 
SUBEND i 


DEF  FNRelhum(Q, T } 

T1=T- .0098*10 

C®2  3.  84--2348/Tl-5.03*LGT(Tl) 
P®1012-.  120*  10 

Relhum®MIN (100, J*10 0/(625*1 0~C/P) ) 

RETURN  Relnum 

SUBEND! 


SUB  We2 ( Zi, Zlcl , Dte , Dqw, Qs , Thr 1 ,0liq,Thl ,Tsfc, Ust,Tst, 
Qst ,T  stv , We2 ,Delr2,Tsky)  !  NEED  TO  CALL  PARCEL 
FIRST 

Dth»Dte-2500*Dqw 
Thr=Thr 1+273. 16 
Th-Thl+273.16 

Epsilon®. 622  1RATIO  OF  WATER  AND  AIR 
Sigma®4.61E-ll  !  STEFANS  CONSTANT  RELATING  TEMP  TO 
BLACK  BODY  RADIATION 
Wl®Qliq*1250 

E  miss®  MAX  (  0,l-EXP(-(Zi-Zlcl)*.138*Wl/2)  )  !  EllISSIVITY 
Tc*Th-.0098*Zlcl  !  TEMP  AT  BOTTOM  OF  CLOUD 
T b®Thr- . 0 0 98 * Z i  !  TEMP  AT  TOP  OF  CLOUD 

Tbar*(Tb+Th)  /2 

Dqsdt® .  622*Qs*  (  2500/  ( .  28  7*Tbar  *Tbar)  )  !  Thr  IS  THETA 
AT  THE  CLOUD  TOP 

Beta*( 1+Th* (1+Epsilon) *Dqsdt) / ( l+2500*Dqsd t) !  L&L/Cp<- 
—2500  Cp»l 

D2*Beta* Dte-Tbar*Dqw  I  IF  NEGATIVE  ,  THEN  CLOUD 

TOP  INSTABILITY 
Dthl»Dth- (Thr-Th ) 

Dl*Dthl* ( 1+EDsilon* 0s) +Epsilon*Tbar*Dqw  !  GV*QS  IN  A 
CLOUD 

See*MIN (1 ,MAX( 0,Zlcl/ Zi) )  !  TRAP  SEE  BETWEEN  0  AND  1 

DEF  FNH vi  (A) ®MAX { 0 ,A )  !  HEAVISIDE  FUNCTION 

Rb»Eraiss*Sigma* (Tb"4-TsXy"4) !  SHOULD  SE  POSITIVE 
UNDER  MOST 

Rc»Emiss*Sigraa*  (Tsfc"4-Tc"4)  !  NORMAL  CONDITIONS 
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4100  It’  See»l  THEN  Rb-Rc*0 
4105  Wth»-Ust*Tst 

4110  Mq»-U3t*Qst/1000  !  USING  GRAMS  PER  GRAM 

4115  .4  thv*— Ust  *Tst v 

4120  Mthe*W  th+2500*Wq 

4125  Ml*{ 2-See)  *See*FNHvi  ( uthv)  +  ( 1-See)  *  ( l-3ee)  * FNHvi  ( Seta* 
Wthe-Th*Wq) 

4130  Ml  *(  2— See)  *See*  FNHvi  (-Nthv)  +  ( 1  -See)  *  (1-See)  *  FNHvi  (- 
Be  t  a*  w  t  he +Th  * Wq ) 

4135  M2»5ee*See*( FNHvi (Rd) -Rc)+(l-3ee) * (1-See) *3eta*Rc+(l- 
3 ee*See)  *Beta*FNHvi  (Rb) 

4140  M2»See*See*FNHvi  (-Rb)  +  (l-5ee*See)  *Beta*FNHvi  (-Rb) 

414  5  .4  3*  (1-See *See)  *  FNHvi  (-D2) 

4150  N3*See*See*01+ { l-3ee*See) *FNHvi(D2) 

4155  Ne2*( .2* (M1+M2) -(N1+N2) )/(N3-.2*M3)  !A*.2  RATIO  OF 

ENTRAINMENT  ENERGY 

4150  IF  Ne2>0  THEM  GOTO  4135  !  TO  TOTAL  TKE 

4165  IF  (D2 <0 )  AMO  (3ee<l)  THEM  Me 2* . 005* ( 1 -See* See )! CLOUD 
TOP  INSTABILITY 

4170  IF  (O2<0)  AND  (3ee*l)  AMD  (f.4tnv>0)  THEM  Me 2* 

.001! ENCROACHMENT 

4175  IF  { D2 <0 )  AND  (See*l)  AMD  (Wthv<0)  THEM  Ne2*.0005 
! STABLE  ENCROACHMENT 
4130  IF  Ne2<0  THEM  Ne2*.0001 

4135  Delr 2*Rb-RC  !  DOME  HERE  CALCULATED  FROM  THE  MATER 

CONTENT,  TEMP  OF  SKY,  CLOUD  EMISSIVITY  (THICKNESS 
OF  CLOUD)  TOP  MINUS  BOTTOM 

4190  SUBEND! 


4195  SUB  ParcelO (Tne ,gl0 , Z ,Rh,Tn  ,Qv,Ql , zs) 

4200  Gam*. 0098 
4205  Gamaew*. 00304 
4210  Lcp»2 . 430 
4215  0q*Qt=wl 0 
4220  Zz*0 

4225  GOSUB  Td_£ind 
4230  TdlO-Td 
4235  Zz*Z 

4240  Thr»The-Lcp*Qt 

4245  GOTO  4345 

4250  Theth : Th*Thez-Lcp*l 0 

4255  Cycle: Tdo* Th-Gam*Z 

4260  T*Tx*TdO 

4265  Zz*0 

4270  GOSUB  gvcal 

4275  0v*gvo 

4280  GOSUB  Gqdtcal 

428  5  D the* (The z- (Tb  +273 . 15 ) * ( I+Lcp* Gv/ ( T+27  3 . 1 6 ) ) +273 . 16 ) / ( 
l+Lcp*gv/(T+273.16)+(Tb+273.16)/(T+273.15)*Lcp*Dqdt) 
4290  IF  AB3 ( D the) <.0001  THEM  RETURN 
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4295 
4  300 
4  305 
4310 
431  5 
4320 
4325 

4330 
4335 
4340 
4345 
4350 
4355 
4360 
4  365 
4370 
4375 
4380 
4385 
4390 
4395 
4400 
4405 
4410 
4415 
4420 
4425 
4  4  30 
4435 
4440 
4445 
4450 


4455 

4460 

4465 

4470 

4475 

4480 

4485 

4490 

4495 

4500 

4505 

4510 

4515 

4520 

4525 


Th*Th+Otne 
30 TO  Cycle 
Qqdtcal :'fx*Tdo+. 00 1 
G0SU3  Qvcal 
Dadt*{  Cvb-0v)  /  .001 
RETURN 

Qvcai :  0vb*5.3542*EX?(-5399.286* (l/(Tx+273.I 5)-3.5  905E- 
3)) 

3ig*EXP(-1.1334E-4*Zz) 

sivb*Qvo/Sig 

RETURN 

Zs* (Thr-rdlO) /Gamdew 

IF  ZS<0  THEN  Zs*0 

IF  Z>*Zs  THEN  GOTO  4400 

Th=Thr 

Tx=Th-Gam*Z 

Zz*Z 

GOSUB  Gvcal 
Rn*Qt/Qvb*100 
aq«Qv*at 
GOSUB  Td  find 
GOTO  442(J 
Thez*The 
GOSUB  Theth 
Rh-100 
Td=Th-Gam*Z 
Ql*Qt-Qv 
SUBEXIT 

Td_f ind ;  Sig *EXP ( -1 . 1 384E-4*Zz) 

Qx*Oq*Sig 

Td»5 399. 286/ (21. 064-LOG (Qx) ) -273. 16 
RETURN 
SUBEND! 


SUB  Shrtwv(ThetafQlbm,ZcfZb,QtrDtemp) !  ADDED  ay 
CHRIS  FAIRALL  9FEB81 
IF  (Zc>Zb)  OR  (Theta=PI/2)  THEN  Dtemp*0 
IF  (Zc>*Zb)  OR  (Tneta=PI/2)  THEN  SUBEXIT 
Qlbm*MAX( .00005 ,Qlbm) 

0q*Uloni/2 

DIM  Xnd {4 ) , Dr (4) 

C  ldthk*Z  b-Zc 
Delz*Cld  thk 
Dr(1.0)  *2.5 
Dr (2) *5 
D  r  (  3 )  *  1  0 
Dr(4) *20 

3h*(Qt-aq/CldtnK’»Delz/2)*I.19  !  NAT  SR  VAPOR 
(JO* Qq*  1.19  1LI0UID  WATER 
Rra* (4774. 548* *q) ".333 


1  0  3 


0-  2  3 


4510  FOR  1-1.0  TO  4  !  XND  IS  THE  XLOUD  DROPLET  3PECT 

4  515  Xnd  ( I)  * . 3  29  96 1* EX? (-LOG  (Dr  ( I)  /Rm)  *LOG  (Dr  ( I)  /Rm)  * 

4.4164) *  1E3 
4540  NEXT  I 

4545  CALL  Nannrad <Cldthk,Sh ,Xnd ( *) , Theta ,0 ,Fheat, Ftop, 
Fbotdir , Fbotdif ) 

4550  !  Fbotdir  DOWNWARD  DIRECT  RAD  AT  CLOUD  BOTTOM 

4555  !  Fbotdif  DOWNWARD  DIFFUSE  RAD  AT  CLOUD  BOTTOM 
4560  1  Ftoo  DOWNWARD  DIRECT  RAD  AT  CLOUDTOP 

4565  !  Theta  ZENITH  ANGLE 

4570  Dterap=Delz* Fheat/Zb/ 60 
4575  SUBEND 

4580  SUB  Int779  (Nxy,  X  ,Xarray  (  *)  ,Yar r ay  ( *)  ,  if  ,Slooe , lx) 

4585  OPTION  BASE  1 

4590  Ix*MIN (MAX(1 ,Ix) ,Nxy-l) 

4595  IF  Xarray ( lx) <X  THEN  4615 
4600  IF  Ix<=l  THEN  5020 
4605  Ix=Ix-l 
4610  GOTO  4595 

46  2  5  IF  XOXarray (IX+1 )  THEN  4615 
4620  IF  Ix+l>*Nxy  THEN  5020 
4625  I X  *1 x+1 
4610  GOTO  4595 

4615  Slope® ( Yarray ( Ix+1 ) -Yarr ay (Ix))/(Xarray(Ix+l) -Xarray ( 
lx)) 

4640  Y=Slope* (X-Xarray ( lx) ) +Yarray ( I x ) 

4645  GOTO  4655 
4650  Ix—999 

4655  SUBEND 

4660  I  SUB  A52  -  INITIALIZE  XKV  ARRAY 
4665  SUB  A52 (Xkv( *) ,U) 

4670  !  DISP  * A52" 

4675  !  U  IS  PATH  IN  G/CM~2 
4680  OPTION  BASE  1 
4685  MAT  Xkv-ZER 
4690  X  kv ( 6 )  = . 0  5/U  “ .  1 9 !  .  0  2 2 
4695  Xkv(  7)  ».09/U“.  19!  .094 
4700  Xkv(8)*.26/U“.19J .169 

4705  Xkv(9) *.28/U~. 26 ! . 145 
4710  Xkv(10)  *.  16/U  '*,261.212 

4715  Xkv(  1 1)  ® .  66/U  “.511.591 
4720  Xkv(  1 2)  ® .  97/U  511.941 

4725  Xkv(ll)*l.  06/U\ 51!  .941 
4710  Xkv(14)*.19/U~. 681.76 

4715  Xkv(15)».ll/U“.681.585 

4740  SUBEND 

4745  !  SUB  A20  -  COMPUTE  NEEDED  ELEMENTS  FOR  AUGMENTS 

MATRIX  A 

4750  SUB  A20  (A(*)  , T , Al , Be , P ,X k, Fi  (*)  ,0ext  (*)  ,Dr  (*)  ,D,Sa(*)  , 
Af ( * ) ,Sh, Xkv ( *) , Xnd (* ) , Ag , I w, xmu ,Tau (* ) ) 

4755  OPTION  BASE  1 
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4760 

4765 

4770 
4775 
4780 
4785 
4790 
4795 
4800 
4805 
4810 
4315 
48  20 
4825 
48  30 
48  35 
4840 
48  4  5 
4850 
4855 

4860 
4865 
48  70 

4875 

4880 

4885 

4890 

4895 

4900 

4905 

4910 

4915 

4920 

4925 

4930 

4935 

4940 

4945 

4950 

4955 

4960 

4965 

4970 


!  DISP  "A20" 

CALL  A28 (G,r,Om, Dr (*) ,Xnd(*) ,Jext(*) ,D,Sa ( *) ,Af (* ) ,Sh, 
Xkv(  *)  ,Iw,Tau(*)) 

Fo-Fi(Iw) 

Fq=G*G 

G»G/(1+G) 

T*  { 1-Om*  Fq)  *T 

Om» ( 1-Fq) *Om/(l-Om*Fq) 

Dl*l-Om 
D2*l-0  m*G 
Xk*(  3*D1*D2)  *  .5 
P* ( 3* 01/D2 )  “.5 
D3*l-Xk*X  K*Xmu*X;nu 

Al=3*Om*  (  Fo/PI)  *Xmu*Xmu*  (1+G*D1 )  /(4*D3) 

3e=3*Om*  (Fo/PI)  *Xmu*  ( 1  +  3*G*D1*X.tiu*X;tiu )  /  ( 4 *D3 ) 

A  ( 1 , 1)  =  1  +  2*  P/3 
A(1 , 2)  =l-2*P/3 
A  ( 1 , 3)  =Al+2*Be/3 

A  (  2 , 1 )  =  ( l-Ag-2  * ( 1+Ag ) *P/3 ) *EXP(-Xk*T) 
A(2,2)=(l-Ag+2*(1+Ag) *  P/3)* EXP (Xk*T) 

A(2, 3) »( ( 1-Ag)  *Al-2*  (1+Ag)  *Be/3+Ag*X.Ttu*  (Fo/PI)  )  *EXP(- 
T/Xmu ) 

SUBEND 

I  SUB  A30  -  COMPUTE  NET  FLUXES  -  FTOP  AND  ARRAY  FL 
SUB  A 30 (Fl, Fbotdir ,Fbotdif ,Ftop,Ftopd,X {* )  ,Xk,Fi ( *)  ,T  , 
Al,3e , P, Iw,Xmu ) 

OPTION  BASE  1 
!  DISP  "A30" 

T0=»0 

I0-X(  1) *EXP(-Xk*T0) +X  { 2) *EXP (Xk*T0) -Al* EXP (-T0/Xmu ) 
I1«P*(X( 1) *EXP(-Xk*T0)-X(2) *EXP ( Xk*T0) ) -Be*EXP(-T0/ 
Xmu) 

Fd*Fi( Iw) *Xmu  1  DOWNWARD  COMP  AT  CLOUD  TOP  IW  BAND, 

PERPIND  TO  HORIZONTAL 

Ftop-Ftop-P I* (10-2/3*11) +Fd  i NET  (D-U)  AT  CLOUD  TOP, 
TOTAL 

Ftopd-Ftopd+Fd  J TOTAL  DOWNWARD  COMP 
I0=*X  (1)  *EXP(-Xk*T)  +X  (2)  *EXP  (Xk*T)-Al*  EXP(-T/Xrau) 
I1«P*(X( 1) *EXP ( -Xk*T) -X ( 2) *  EXP ( Xk*T ) ) -Be*EX  P  (-T/Xnui ) 
Fl»Fl+4/3*PI*Il  !NET  AT  LEVEL 

F1»F1+Fd* EXP (-T/Xmu)  (TOTAL  NET  (DIF+DIRECT)  LEVEL  I 
Fbotdi£*Fbotdif+PI* (10+2/3*1 1) 

Fbotdir»Fbotdir+Fd*  EXP (-T/ Xmu) 

SUBEND 

1  SUB  A 28  -  AVERAGE  3A  ,AF  ARRAYS  AND  STORE  IN  OM,G; 
COMPUTE  T 

SUB  A28  (G  ,T  ,Om,Dr  (*)  ,xnd(*)  ,Jext  (*)  ,  D  ,Sa  (  * )  , Af  ( * )  ,Sh, 
Xkv  ( * )  ,  I w  ,T  au  (  * ) ) 

!  DISP  "A 28  " 

OPTION  BASE  1 
Om*G=Ts*0 
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! 

li 
*  \ 
'  i 


4975  FOR  J-l  TO  4 

4980  Tau(J)  -Cext  (Iw,  J)  *PI*  (Dr  ( J)  ME-6)  ~2*Xnd  (J)  *D 
4985  Ts-Ts+Tau(J) 

4990  NEXT  J 

4995  T-Ts 

5000  FOR  J-l  TO  4 

5005  Om»Om+Sa( Iw,J) * (Tau(J) /Ts) 

5010  G-G+Af ( I w, J ) *Tau ( J) /Ts 

5015  NEXT  J 

5020  Vp-S  h*D* 1 E-4 

5025  On-Om/(l+Xicv(Iw)  *V?/T) 

5030  3UBEND 

5035  !  SUB  A100  -  COMPUTE  JEXT,SA,AF  ARRAYS  -  ROUTINE  FROM 

A.  GO ROC rt 

5040  SUB  A100 (Qext( * ) ,  Sa  ( * ) ,Af  ( *) ,  Dr ( *) ) 

5045  i  DISP  " A1 00  M 

5050  OPTION  BASE  1 

5055  DIM  X 1 9 ( 1 5 ) ,  Xn  9 ( 1 5 )  ,  X  K  9  ( !  5 ) 

5060  DATA  0.25,0.35,0.45,0.55,0.65,0.75,0.85,0.95,1.05,1.15 
5055  DATA  1.25,1.35,1.45,1.55,1.65 

5070  DATA  1.362,1.345,1.335,1.333,1.332,1.3,1.329,1.327, 
1.326  1.325 

5075  DATA* 1.32 3, 1.321, 1.320, 1.318, 1.316 

5080  DATA  2.5,5,10,20 

5085  RESTORE  5060 

5090  MAT  READ  Xl9,Xn9,Dr 

5095  MAT  XX9-ZER 

5100  REM  EVALUATE  EXTINCTION  USING  DEIRMENDJIAN 
APPROXIMATION 
5105  FOR  1x9-6  TO  15 
5110  X10-X19 (1x9) 

5115  XnO»Xn9 (1x9) 

5120  FOR  1-1  TO  4 
5125  R2-Dr ( I ) 

5130  RQ-4*PI*R2/X10* (Xn0-1) 

5135  Qsca*Qext (1x9,1 ) =2 

5140  Sa( 1x9 , I ) -Qsca/Qext (1x9,1) 

5145  IF  Sa( 1x9, I) =1  THEN  Sa ( 1x9 ,1 )=. 9999999999 
5150  Af  ( I  x9 ,  I )  -1.  8*  ( .  5+EXP  (**4*R0)  /R0+ (EXP(-4*R0)  -1)  /(4*R0)  “ 
2) 

5155  NEXT  I 
5160  NEXT  1x9 
5165  i  OISP  mm 
5170  SUBENO 

5175  SUB  Kannrad  (Cldthk,Sb  ,Xnd  (  *)  ,Tneta,I  test,  Fheat,  IincO, 
Fbotdir , Fbotdic ) 

5180  l  DISP  " KAHN RAD" 

5185  OPTION  BASE  1 

5190  DIM  A { 2 , 3 ) ,R(3) ,X(2) 

5195  DIM  Qext(15,4)  ,Sa( 15,4) ,Af (15,4) 

5200  DIM  Dr  (4)  , Tan (4)  ,Alam(15) 
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5205  DIM  F i  (15)  ,Xkv(15) 

5210  !  I.  INITIALIZATION 

5215  DATA  .0165, .1580, .2839, .2765, .2318, .18 22  , .2435,. 1149, 
.3948,  .0792 

5220  DATA  .  064  3  ,  .  0  51 3  ,  .  0 4  2  4  ,  . 0  3 48  ,  .  0  288 

5225  DATA  .  2  , . 5 35 5 , . 80 3  , .  909 , .  95 3 , .  960  , . 96 0  , .  920  , .  9 3 3  ,  .  940 , 
.850, ,730, .770, .830, .910 
5230  RESTORE  5215 
5235  HAT  READ  Fi,Alam 
5240  U*MAX(.05,Sh*Cldtttk*lE-4) 

5245  Fbotair»Fbotdif =Ftop=Ftopd=Ftopdir=0 
5250  IF  Itest*l  THEN  Tneta=PI/4 
5255  Xmu=COS (Theta) 

5260  IF  Xmu  =  0  THEN  5285 

5265  3ecmu=l/Xmu 

5270  FOR  1  =  1  TO  15 

5275  Fi(  I)  =  Fi(  I) *Alam( I) "Secmu 

5280  NEXT  I 

52b5  Ag= . 1  !  SURFACE  AL3ED0 

5290  IF  X.nu> .  7  THEN  Ag=.05 
5295  Delz=Cldthk 
5300  F1=0 

5305  D*Delz 

5310  CALL  A100  (Qext  (* ) ,3a (*) , Af ( * ) , Dr { * ) ) 

5315  CALL  A5  2 (Xkv ( * ) ,U) 

5320  !  II.  SOLUTION  OF  EQUATIONS 

5325  FOR  Iw-6  TO  15 

5  3  30  CALL  A20(A(*) ,T  ,A1 ,3e , P  ,Xk, F i  (*)  ,Qext(*) ,Dr(*)  ,D,Sa(*) 
,Af (*) ,Sh,Xkv ( *) ,Xnd( *) ,Ag  ,1  w  ,xmu  ,r  au  (*  ) ) 

5335  !  PERFORM  GAUSS  ELIMINATION  ON  MATRIX  A 

5340  FOR  L=1  TO  2 
5345  I1=L 
5350  Dd=A(L,L) 

5355  IF  ABS (Dd) >.000001  THEN  5405 

5360  11=11+1 

5365  IF  1 1 > 2  THEN  5545 

5370  FOR  Iwi=l  TO  3 

5375  R(Iwi)«A(L,Iwi) 

5380  A ( L  , I  wi)  =A  (  1 1  ,  Iwi  ) 

5385  A( II ,1 wi) =P.( Iwi) 

5390  NEXT  Iwi 
5395  Dd=A ( L , L ) 

5400  GOTO  5355 

5405  FOR  Ibi  =  l  TO  3 

5410  A(L,Ioi)=A(L,Ibi)/Dd 

5415  NEXT  Ibi 

5420  FOR  J=L+1  TO  2 

5425  Dd=A(J,L) 

5430  FOR  X*1  TO  3 

5435  A(J,X)=A{J,X)-Dd*A(L,K) 

5440  NEXT  K 
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544  5  NEXT  J 
5450  J  EXT  L 

5455  FOR  K*2  TO  2  STEP  -1 
5460  FOR  J*K-1  TO  1  STEP  -1 
5465  A(  J,  3)*A(  J,  3) -A(J,K)  *A(K,3) 

5470  A  ( J ,  X )  *0 
54  75  NEXT  J 
5480  NEXT  X 
5485  FOR  1*1  TO  2 
5490  X  ( I ) *A  ( I  ,  3 ) 

5495  NEXT  I 

5500  CALL  A30 ( FI, Fbotdir ,  F  bo  tdif  , Ft op , Ftopd ,X ( * ) , Xk , Fi ( * )  , 

T  ,A1 ,  Be ,  P, I w, Xmu ) 

5505  NEXT  Iw 

5510  !  III.  OUTPUT  MODULE 

5515  Fheat*(Ftop-Fl) /D*  35. 014 
5520  Iinc0*Ftopd 
5525  F top* I incO* I incO/1 . 7 1 36 
5530  Fbotdir=Pbotdir/l .7! 36 
5535  Fbotd ir*Fbotdif /I .7136 
5540  SUBEND 

5545  PRINT  "INDETERMINATE” 

5550  SUBEND 

5555  SUB  S ky2 ( Zi,Nsnd, Pot_temp ( *) , Zsnd( *) ,Spec_num( * ) ,Tsky , 
F  tOD } 

5560  DATA  0 ~,  .0 0001 ,.  000021 5  ,.  000046 4  ,.  0001 ,.  00 02 1 5 ,.  000464  , 
. 001,. 00 215,. 00464, .01,. 021  5, .0464  ,. 1  ,.2 15,. 464,1  , 
2. 15,4. 64, 10, 21. 5, 46. 4, -9 99 

5565  DATA  0,1.86,2.58,4.11,5.72,7.81,11.4,14.6,13.3,23.6, 
27.7,31.9,37.  4,41.7,46.2,52.  9,59,66.6,7  8.8,  88.1  , 
95.1,98.3,-999 

5570  DATA  . 0093 , 1 . 22 5 , 10094 , 4. 61E-1 1 , 0 , 1 , 0 
5575  OPTION  BASE  1 

5530  DIM  Tu(23) ,Temis(23) ,Smit(2) ,Tnrsnd(30) ,0tsnd(30) 

5585  RESTORE  5560 

5590  MAT  Qtsnd*5pec_num 

5595  MAT  Thrsnd*Pot”temp 

5600  MAT  READ  Tu,Temis 

5605  READ  Gad, RhoQ , Zscale ,Sig  ,U, Ii  ,3 

5610  Nemis*22 

5615  !  TRRSNDO  IS  TH  RADIOSND 

5620  !  ZSND ( )  IS  ALTI  IN  M , RAQSN 

5625  l  MEM I 3  IS  #  OF  EMISSIVITY  POINTS  IN  TEMI35 

5630  1  T U ( 2 3 ) ,U  DATA 

5635  !  TEMIS ( 23) , EMISSIVITY  DATA 

5640  !  ITOP  IS  INDEX  FOR  Z3ND=ZI  ,  IE  WHAT  LEVEL  IN  RAOSMD 

IS  AT  INVERSION 

5645  !  NSND  IS  NUMBER  IS  POINTS  IN  RAD IOSOND 

5650  l  ICLOUD  IS  INDEX  OF  UPPER  LEVEL  CLOUDS 
5655  Icloud*N snd  ! ASSUME  NO  UPPER  LEVEL  CLOUDS 
5660  Nsna=Nsnd-l 
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5665  IF  Msnd<*0  THEN  SUBEXIT 

5670  FOR  Itop*l  TO  Nsnd 

5675  IF  Zsnd { I  top) >Zi  THEM  GOTO  5685 

5680  NEXT  I  top 

568  5  I  top* I  top-1 

5590  DI3P  HSKY2H 

559  5  Tts*Thrsnd{ I  top) -Gad*  Zsnd ( I top) 

5700  Emit(li)*0 

5705  FOR  1*1 top+1  TO  Nsnd-1 

5710  1 1*3-1 i 

5715  IF  I>Icloud  THEN  GOTO  5735 

57  20  U=U+Qtsnd(  I)  *  ( .00  00  5*  ( Zsnd ( 1+1) -Zsnd ( 1-1)  )  *Rno0*£XP(- 
Zsnd(I) /Zscala) ) 

5725  CALL  I nt 779 ( Nemis , U , Tu ( * ) , Temis ( * ) , Emi t ( I i ) , Es , I e) 
5730  Tt*TPr snd ( I) -Gad* Zsnd ( I ) 

573  5  3=S-Smit (Ii) *Sig* ((Tt+273.16) *4-(Tts+273 . 1  6) “4) /I  00 
5740  Tts*Tt 
5745  NEXT  I 

5750  Tt*Thrsnd (Msnd )-Gad* Zsnd (Nsnd) 

5755  U*U+Qtsnd(Nsnd)*,0001* ( Zsnd (Nsnd) -Zsnd (Nsnd-1 ) ) *Rho0* 
EXP (-Zsnd (Nsnd ) /Z scale) 

5760  CALL  Int 779 ( Nemis , U ,Tu (*), Temis (*), Emit ( I i ), Es, Ie) 

5765  3*S+Sig* (Tt+273.16) *4*Emit(Ii) /100 

5770  Icloud=Nsnd 

5775  F  too*  0 

5780  GOTO  5795 

5735  Tt=Tiirsnd(  Icloud) -Gad* Zsnd  ( Icloud) 

5790  Ftop=Sig* (Tt+273. 16) "4 

5795  Tsky* ( ( Ftop* ( 1-Emit  (Ii) ) +S) /Sig) 25 

5800  DISP  "" 

5805  SUBEND 

5310  SUB  Zenith(Time,Julday,Lat,Tneta) 

5815  RAD 

5820  Hr *T ime/3600  MOD  24-12 

5825  Latr*Lat*. 017453 

5830  H=Julday* .017214 

5835  Decl*.006913-. 39991 2 *COS (H) +. 0 70257*3 IN( H) -. 006 753* 

CO S  (  2 * H ) +.000  907*SIN (2*H) -.002  697*005 (3* H)  + 
.000143*3IN(3*H) 

5340  Theta*ACS (SIN ( La  tr ) *3 IN (Decl) +C0S (Latr ) *C03 (Decl) * 

COS (Hr  *PI/1 2) ) 

5345  IF  Tneta>?I/2  THEN  Theta=PI/2 
5850  3U3SND 


APPENDIX  E: 


MAJOR  PROGRAM  VARIABLES  IN  PREDICTION  CODE: 


Alpha  t  ...  ratio  of  heat  transfer  to  momentum  transfer 
excnange  coefficient  for  neutral  conditions  (Psi*0) 

Aws  ...  subsidence  rate  value  at  initial  inversion 
heignt  [m/s] 

Cdn  ...  neutral  drag  coefficient 

Ctn_sqrt  ...  square  root  of  the  temperature  structure 
function  parameter 

Ddqwdt  ...  rate  of  cnange  of  specific  humidity  within  tne 
layer  [gAg-s] 

Ddtedt  ...  rate  of  change  of  equivalent  potential 
temperature  within  the  layer  [K/s] 

Delr  ...  radiation  flux  divergence  at  cloud  top  [km/s] 
Deltim  ...  prediction  time  step  [1800  sec] 

Dhdt  ...  rate  of  change  of  the  mixed  layer  height  and  is 
the  sum  of  entrainment  and  subsidence  effects  [m/s] 
Dqdzu  ...  lapse  rate  of  specific  humidity  above  the 
layer  [g/kg-km] 

Dqw  ...  jump  in  specific  humidity  at  the  inversion  [g/kg] 
Dqwdt  ...  rate  of  cnange  of  specific  humidity  jump  [g/kg-s] 
Dtdzu  ...  lapse  rate  of  potential  temperature  above  tne 
layer  [K/km] 

Dte  ...  jump  in  equivalent  potential  temperature  at  the 
layer  [ K ] 


Otedt  ...  rate  of  change  of  equivalent  potential  temperature 
Dtedzu  ...  lapse  rate  of  equivalent  potential  temperature 
above  tne  layer  (k/km) 

Oth  ...  jump  in  potential  temperature  at  the  layer  (k] 

Kkk  ...  Von  Karman  constant  (.35) 

Lop  ...  liquid  water  content 

Mixq  ...  initial  value  of  mixed  layer  specific  humidity 
Mixtemp  ...  initial  value  of  mixed  layer  potential 
temperature 

Nrll  ...  beginning  in  half  hour  period  [e.g.  1200s*24] 

Nrad  ...  100  times  the  ratio  of  buoyancy  in  clouds 
Press  ...  surface  pressure,  default  to  1012  mb  unless 
entered  from  IREPS  [mb] 

Psil  ...  stability  correction  to  logarithmic  surface  layer 
wind  profile.  Determined  from  surface  stability. 
(Businger  et  al ,  1971) 

Psi2  ...  stability  correction  to  logarithmic  surface  layer 
temperature  and  humidity  profiles,  determined  from 
surface  stability.  (Businger  et  al,  1971) 

Q10  ...  Op  ...  computed  well-mixed  specific  humidity  wnich 
is  assigned  to  the  10  meter  level  [g /kg ] 

Qst  ...  surface  layer  specific  numidity  scaling  parameter, 
-Wq/Ust,  [g/kgj 

Qw  ...  total  (liquid  +  vapor)  specific  humidity  [g/kg] 

5  ...surface  layer  stability  (S*Z/L,  L=Monin-Obuhxov  lengtn) 
no  ...  well-mixed  potential  temperature  wnich  is  assigned 
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to  the  10  meter  level  [K] 

Te  ...  equivalent  potential  temperature 

Th  ...  well-mixed  potential  temperature  in  degrees  Kelvin 
fsfc  ...  surface  temperature  [K] 

Tst  ...  surface  layer  temperature  scaling  parameter  [K] 

Tstv  ...  surface  layer  virtual  potential  temperature  scaling 
parameter  [K] 

U10  ...  well-mixed  wind  speed  which  is  assigned  to  the  10 
meter  level  [m/s] 

Ust  ...  surface  layer  momentum  scaling  parameter  [m/s] 

Vdte  ...  advection  above  tne  layer 
we  ...  entrainment  rate  [m/s] 
ws  ...  subsidence  rate  [m/s] 

wq  ...  surface  flux  of  specific  humidity  [g/g] 

wst  ...  mixed  layer  convective  scaling  velocity  [m/s] 

Wto  ...  surface  flux  of  temperature  [K-m/s] 

Wtv  ...  surface  flux  of  vertical  potential  temperature 
(buoyancy  flux)  [K-m/s] 

Zi  ...  heignt  of  the  layer  [m] 

Zlcl  ...  lifting  condensation  level  [m] 

Zo  . . .  surface  roughness  parameter  for  momentum  [m] 

Zot  ...  surface  roughness  parameter  for  temperature  and 
humidity  [m] 

Zws  ...  initial  height  used  for  subsidence  rate  referencefm] 


VARIABLE  CROSS  REFERENCE 


MAIN  PROGRAM 
COMMON  VARIABLES: 


Air 

115 

1070 

1850 

Change 

115 

1070 

Corel  ( *) 

120 

Delta 

115 

1070 

1850 

EnvnamS (*) 

125 

945 

1005 

1055 

Envaq( *) 

120 

930 

965 

995 

1055 

1060 

EvapS 

95 

1060 

Heights 

95 

Height(* ) 

110 

1065 

1080 

1090 

1095 

1455 

1465 

1480 

1620 

1630 

1365 

Htzero 

115 

1070 

LocS 

95 

440 

490 

530 

1060 

1335 

1845 

1940 

2630 

2900 

LosnamS (* ) 

125 

Lossq ( * ) 

120 

930 

Munits  (* ) 

110 

1065 

1865 

Names 

95 

1060 

Nmax 

115 

1070 

1075 

1080 

1090 

10  95 

1455 

1460 

1465 

1480 

1620 

1625 

1630 

1865 

1955 

Nunits(*) 

110 

1065 

Presur ( *) 

105 

1065 

1080 

1865 

Przero 

115 

1070 

1110 

1850 

Recenv 

115 

Relhm 

115 

1070 

1850 

Relhum(* ) 

105 

1065 

1085 

1865 

Sea 

115 

445 

340 

845 

8  50 

1070 

1850 

2370 

SysnamS  (* ) 

125 

Syssq( *) 

120 

930 

Temper ( *) 

105 

1065 

1085 

1090 

1365 

TimeS 

95 

440 

490 

530 

1060 

1335 

1845 

1945 

2630 

2900 

TypeS 

95 

1060 

Wind 

115 

575 

580 

1070 

1850 

WmoS (*) 

100 

1060 

WmohtS 

100  ' 

1060 

VARIABLES: 

Aa 

265 

285 

290 

295 

300 

Alpha_t 

190 

2330 

2430 

2470 

2490 

Aws 

2130 

2280 

2285 

Cc 

3000 

3005 

3010 

3015 

Cdn 

2350 

2355 

2375 

2385 

2430 

2470 

Ctn  sqct 

2380 

2385 

D  " 

1780 

1785 

Date? 

220 

Ddqwd  t 

2235 

2265 

Ddtedt 

2235 

2275 

Dele 

2235 

Deltim 

190 

2245 

Dhdt 

2235 

2245 

Dig  it_top 

1100 

1460 

Dqdzu” 

1920 

2115 

Dqp 

1915 

2110 

Dqw 

2120 

2235 

Dqwdt 

2235 

2250 

Dtdzu 

1930 

2115 

Dte 

2110 

2235 

Dtedt 

2235 

2260 

Dtedzu 

2115 

2205 

Dth 

1925 

2110 

Dummy 

1035 

1040 

Flag_data 

360 

1130 

Flag~rev 

310 

365 

Frame_date$ 

155 

1945 

Frame_loc$ 

155 

1940 

Ftop  “ 

1955 

G(*) 

130 

550 

1120 

2240 

2345 

2805 

G1 

2645 

2650 

G2 

2650 

2670 

G3 

2655 

2670 

G4 

2660 

26  70 

G5 

2665 

2670 

Hd$ 

430 

470 

2625 

3040 

Hight$ 

106  0 

I 

940 

945 

1010 

1090 

1870 

2030 

2945 

Idayl 

755 

765 

Ier 

910 

915 

Inll$ 

535 

540 

625 

Inl2$ 

645 

650 

J 

160 

960 

2075 

2430 

1790 

2470 

2250 

2260 

2265 

2275 

2285 

1625 

1860 

2205 

2220 

2565 

2120 

2270 

2265 

2270 

2310 

2205 

2565 

2275 

2310 

2235 

2310 

1045 

2525 

1125 

2565 

2565 

1050 

630 

675 

745 

770 

780 

2025 

2050 

2060 

2070 

2215 

2245 

2255 

2270 

2  30  5 

2310 

2565 

2655 

2660 

2665 

2670 

2825 

2855 

2870 

2940 

2950 

2655 

2660 

2665 

2670 

2675 

1835 

2585 

2590 

2595 

2600 

2690 

2395 

2905 

3020 

3030 

950 

960 

990 

995 

1005 

1015 

1020 

1075 

1080 

1085 

1095 

1100 

1105 

1860 

1865 

2040 

2045 

2050 

2055 

2060 

2085 

2810 

2825 

2855 

2875 

2950 

2955 

775 

920 

925 

545 

550 

565 

595 

610 

655 

6  60 

665. 

670 

675 

965 

970 

2020 

2065 

2070 

2175 

2135 

2210 

2215 

2225 

114 
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2240 

2245 

2255 

2270 

2  30  5 

2310 

2325 

2345 

2850 

2855 

2365 

Juld ay 

170 

445 

490 

495 

2230 

K 

2795 

2805 

2825 

2830 

2845 

2855 

2870 

2885 

Kkk 

190 

2  37  5 

2380 

2385 

2430 

2470 

2490 

2495 

L 

1235 

1240 

1250 

1255 

1270 

1275 

1280 

1285 

1295 

1300 

1305 

1310 

1430 

1435 

1440 

1445 

1460 

1465 

1470 

1625 

1630 

1635 

2825 

2830 

2835 

2840 

L$ 

165 

2630 

2635 

Lat 

17  5 

445 

500 

505 

2230 

Lcp 

2095 

2105 

2110 

2115 

2205 

2300 

2  30  5 

2310 

LI 

1265 

1275 

1280 

1290 

Mass$ 

75 

910 

920 

Mhr 

1960 

1980 

2175 

2565 

2645 

Mixq 

1650 

1660 

1665 

1685 

1725 

1730 

1910 

1915 

Mix  temp 

1490 

1515 

1525 

1540 

1575 

1580 

1900 

1925 

Mmm 

1530 

1650 

1675 

1780 

Mult 

2490 

2500 

2505 

Nd 

2055 

2060 

2065 

Nrl 

2035 

2040 

2055 

2060 

2070 

2080 

Nrll 

565 

620 

630 

675 

745 

755 

770 

780 

1980 

1990 

2020 

2025 

2035 

2145 

2150 

2155 

2175 

2210 

2215 

'565 

2645 

2735 

2740 

2745 

2755 

2765 

2775 

2785 

2805 

2825 

2830 

28  50 

2870 

2915 

2925 

2935 

2940 

2945 

2965 

2980 

2990 

3000 

Nrec 

1020 

1030 

1045 

Numhr 

1980 

2085 

02 

2220 

2235 

P 

2815 

2855 

2860 

Pot^temp (*) 

165 

1090 

1455 

1465 

1480 

1865 

1955 

Press 

1110 

1115 

2565 

Ps  il 

2405 

2420 

2430 

2465 

2470 

2495 

Ps  i2 

2410 

2425 

2430 

2460 

2470 

2490 

Q$ 

570 

585 

590 

635 

680 

695 

720 

7  25 

730 

735 

Q10 

2195 

2235 

2255 

Cdelta 

2370 

2385 

2505 

Qp 

1910 

2105 

2125 

2195 

2255 

2370 

Qstar 

2235 

2505 

2510 

Qw 

2125 

2250 

2255 

2305 

Rc$ 

355 

390 

400 

405 

410 

415 

Rec 

995 

1000 

1005 

1050 

1055 

1060 

w 


115 
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Run$ 

2170 

2210 

S 

2385 

2400 

2450 

SI 

2385 

2435 

So 

2385 

2390 

Spec  hum(*) 

165 

1085 

SticE$  (*) 

155 

190 

T 

20  90 

2095 

T10 

1900 

2015 

Tday$ 

220 

225 

2565 

Tdelta 

2365 

2385 

Te 

2105 

2260 

Th 

2015 

2090 

2305 

Theta 

2230 

2235 

Timesec 

2225 

2230 

Top_inv_q 

1675 

1685 

Top_inv_temp 
Top”q  ~ 

1530 

1700 

1540 

1710 

Top_q_alt 

1700 

1710 

Top_temp 

1555 

1565 

Top” temp  alt 

1555 

1565 

Ts  f  c 

350 

2235 

Tsky 

1955 

2235 

Tstar 

2235 

2500 

Tstv 

2235 

2510 

Tt 

2100 

2105 

2305 

U10 

2345 

2350 

Os  tar 

2235 

2495 

Vdt 

1995 

2005 

Vdte 

2005 

2275 

Vss 

2555 

2565 

W  in 

455 

460 

635 

765 

300 

Winl 

750 

765 

Wind  beginl 

575 

600 

WindTn(*) 

155 

460 

710 

1250 

Wn 

455 

56  0 

Ws 

20  5 

1935 

Wsl 

205 

445 

Wtv 

2235 

XI 

2450 

2465 

X2 

2455 

2460 

X5 

2060 

2070 

Z 

190 

2375 

2215 

2420 

2425 

2430 

2435 

2440 

2455 

2470 

2475 

2480 

2440 

2475 

2480 

2395 

2430 

2470 

1620 

1440 

1630 

1645 

1865 

1955 

2100 

2295 

2300 

2190 

2235 

235 

240 

445 

1335 

1975 

2630 

2500 

2305 

2105 

2110 

2115 

2190 

2205 

2310 

2365 

2385 

2510 

1690 

1710 

1735 

1915 

1920 

1545 

1565 

1585 

1925 

1930 

1740 

1920 

1740 

1920 

1590 

1930 

1590 

2365 

1930 

2510 

2110 

2310 

2115 

2205 

2290 

2295 

2355 

2335 

2495 

2570 

465 

640 

665 

670 

675 

690 

705 

710 

715 

760 

770 

775 

780 

785 

795 

770 

605 

610 

625 

630 

520 

565 

630 

670 

675 

745 

765 

770 

775 

780 

685 

2130 

705 

2235 

760 

2280 

800 

1235 

825 

830 

1415 

1935 

2380  2385  2490  2495 


116 


F-  4 


Zi 

1475 

1490 

1645 

1905 

2090 

2235 

Zil 

2130 

2290 

Zi3 

1475 

1715 

Zlcl 

2235 

2240 

ZO 

2375 

2495 

Zot 

190 

2330 

Zsnd( 

*) 

165 

1095 

Zws 

2135 

2285 

USER  DEFINED  FUNCTIONS: 

FNQvalue 

1085 

2370 

FNRelhum 

2665 

2940 

SU3  PROGRAMS: 

Frame 

2565 

M2 

2235 

Sky2 

1955 

Zenith 

2230 

JlJKP  TARGETS  : 

155 

80 

190 

185 

220 

Start: 

875 

255 

Option: 

260 

Menu: 

290 

295 

1040 

310 

Start_out: 

300 

315 

Start”out_ag 

405 

350 

RevcftangeT 

300 

400 

425 

List^cur  r : 

415 

490 

Jullat: 

320 

415 

515 

Enter_winds : 

325 

415 

530 

525 

540 

610 

590 

645 

One_wind: 

655 

660 

695 

Print_wind: 

650 

740 

745 

730 

795 

Fixwind: 

735 

300 

815 

Subsidence: 

330 

415 

840 

Sea  temp: 

335 

415 

860 

Assxg_error : 

905 

880 

No^data: 

860 

915 

900 

Ireps^data: 

245 

960 

All  out: 

935 

975 

All~full: 

955 

980 

Disp^data: 

315 

415 

1020 

Pick"”data: 

1000 

1505 

1525 

1565 

1580 

1585 

1665 

1710 

1715 

1730 

173  5 

1920 

1930 

1955 

1990 

2025 

2100 

2130 

2135 

2130 

2185 

2245 

2230 

2290 

2295 

1905 


2490 

1955 


2950 


300  305  360  410  415 

2335  2525  2570 


420 


305 

810 


925 


104  5 


1030 


117 


F-  5 


1110 

1140 

Plt_str uct : 

1770 

1885 

1900 

Run : 

2045 

2175 

Loope  r : 

2340 

Bulk: 

2420 

2450 

2490 

S  tars: 

2520 

Picture : 

2585 

Hard_output 

2610 

GvalT 

2710 

Mix_plot : 

3065 

Bye? 

3085 

30  95 

3110 

En: 

1080 

340 

415 

1790 

1785 

1835 

310 

345 

2050 

2085 

2030 

2320 

2395 

2445 

2390 

2435 

2415 

2435 

2475 

300 

2550 

2570 

425 

1830 

2605 

2570 
2570 
30  0 
3075 
3085 

2580  2705  3040 

2615  2890 


SOB  PROGRAM  M2 
COMMON  VARIABLES: 


VARIABLES: 


Cl 

3130 

3150 

Ddqwdt 

3115 

3205 

Ddtedt 

3115 

3200 

Delr 

3115 

3180 

3200 

3210 

Dhdt 

3115 

3195 

Dqdzu 

3115 

320  5 

Dqw 

3115 

3180 

3205 

3215 

Dqwdt 

3115 

3215 

Qte 

3115 

3180 

3200 

3210 

Dtedt 

3115 

3210 

Dtedzu 

3115 

3200 

Otemp 

3185 

3200 

3210 

Q10 

3115 

3155 

3160 

3135 

Q1 

3160 

3170 

3185 

Qliq 

3170 

3130 

Qs 

3165 

3180 

Qst 

3115 

3145 

3180 

Qv 

3160 

3165 

Rh 

3160 

3165 

T1Q 

3115 

3120 

3155 

3180 

Te 

3155 

3150 

Thet 

3125 

3150 

3200 

3210 

Theta 

3115 

3185 

Thr 

3160 

3130 

Tsfc 

3115 

3130 

F 


3060 


113 


o 


Tsky 

3115 

3130 

Tst 

3115 

3135 

3180 

Tstv 

3115 

3140 

3130 

Tt 

3120 

3125 

3130 

3150 

Ust 

3115 

3135 

3140 

3145 

we 

3190 

3195 

3200 

3205 

we  2 

3180 

3190 

Wq 

3145 

3150 

3205 

3215 

ws 

3115 

3195 

Wt 

3135 

3150 

wte 

3150 

3200 

3210 

wtv 

3115 

3140 

Zi 

3115 

3120 

3125 

3160 

3205 

3210 

3215 

Zlc  1 

3115 

3160 

3130 

3185 

3185 


USER  DEFINED  FUNCTIONS: 


SUB  PROGRAMS: 
ParcelO 
Shrtwv 
We  2 


3160 

3135 

3180 


JUMP  TARGETS: 


FNQvalue 

COMMON  VARIABLES: 
VARIABLES: 


Qvb 

3230 

3235 

R 

3225 

3230 

Tx 

3225 

3230 

USER  DEFINED  FUNCTIONS: 
'  SUB  PROGRAMS: 

JUMP  TARGETS: 


SUB  PROGRAM  Frame: 

COMMON  VARIABLES: 

VARIABLES: 

A 

Bb 

Bbb 

Cc 


3755 

3760 

3370 

3380 

3385 

3375 

3385 

3400 

3360 

3365 

3370 

3390  3400 

3375  3395  3405 


3200 


119 


F-  7 


Cv 

3410 

3415 

3590 

3420 

3425 

3575 

3580 

D 

3525 

3530 

3535 

3540 

Dqdzu 

3245 

3640 

Dtdzu 

3245 

3645 

E 

3660 

3675 

Frame  (* ) 

3245 

Frame_date$ 

3245 

3310 

Frame“loc$ 

3245 

3310 

G(*>  ” 

3245 

3430 

3445 

3460 

3470 

3485 

3625 

3630 

3640 

3645 

3755 

dd$ 

3735 

3790 

3795 

Hrg 

3615 

3625 

3630 

3640 

3645 

3755 

I 

3620 

3640 

3645 

3655 

3665 

3675 

3710 

3715 

3720 

Mhr 

3245 

3360 

3440 

3465 

3490 

Nr  11 

3245 

3270 

3275 

3280 

3285 

3305 

3330 

3350 

3360 

3415 

3430 

3460 

3465 

3  48  5 

8490 

3615 

P 

3655 

3660 

3675 

Pm 

3610 

3695 

3930 

3705 

3715 

3735 

3740 

Pp(*) 

3675 

3705 

3715 

3725 

3735 

3740 

Ppl 

3640 

3660 

Pp2 

3645 

3665 

Pr 

3605 

3610 

3615 

3765 

Press 

3245 

3655 

Sbd 

3565 

3315 

3940 

T 

3665 

3675 

TdayS 

3245 

3300 

V 

3440 

3445 

3450 

3465 

3470 

3475 

3495 

3500 

Zg  (*) 

3255 

3595 

3600 

3625 

3630 

3640 

3655 

3665 

3675 

3705 

3715 

3740 

3930 

USER  DEFINED  FUNCTIONS: 
FNRelhum  3755 

SUB  PROGRAMS : 

JUMP  TARGETS: 

3510  Profile: 

3840  3815 

3915  Sfc  duct:  3725 


FNRelhum 

COMMON  VARIABLES: 
VARIABLES: 


3585 


3495 

3680 

3320 

3440 

3750 

3930 


3490 

3645 

3735 


120 


F-  3 


F 


V* 

3975 

3985 

p 

3930 

3985 

0 

3965 

3985 

Relhum 

3985 

3990 

T 

3965 

3970 

T1 

3970 

3975 

USER  DEFINED  FUNCTIONS: 


SUB  PROGRAMS: 
JUMP  TARGETS : 


SUB  PROGRAM  We2 
COMMON  VARIABLES: 


VARIABLES: 


A 

4085 

3eta 

4060 

4065 

Dl 

4075 

4150 

D2 

4065 

4145 

Delr2 

4000 

4135 

Dqsdt 

4055 

4060 

□qw 

4000 

4005 

Dte 

4000 

4005 

Dth 

4005 

4070 

Dthl 

4070 

4075 

Emiss 

4035 

4090 

Epsilon 

4020 

4060 

Ml 

4125 

4155 

M2 

4135 

4155 

M3 

4145 

4155 

NX 

4130 

4155 

N2 

4140 

4155 

N3 

4150 

4155 

Qliq 

4000 

4030 

Qs 

4000 

4055 

Qst 

4000 

4110 

Rb 

4090 

4100 

Rc 

4095 

4100 

See 

4080 

4100 

4150 

3  igma 

4025 

4090 

Tb 

404  5 

4050 

Tbar 

4050 

4055 

Tc 

4040 

4095 

Th 

4015 

4040 

Till 

4000 

4015 

Thr 

4010 

404  5 

4125 

4130 

4135 

4140 

4150 

4165 

4170 

4175 

4065 

4075 

4065 


4095 

4075 


4075 

4135 

4140 

4135 

4135 

4135 

4125 

4130 

4135 

4140  4145 

4165 

4095 

4090 

4170 

4175 

4055 

4075 

4050 

4060 

4070 

4125  4130 

4070 

121 


F-  9 


Thcl 

4000 

4010 

Tsf  c 

4000 

4095 

Ts  ky 

4000 

4090 

rst 

4000 

4105 

rstv 

4000 

4115 

Ust 

4000 

4105 

4110 

4115 

We2 

4000 

4155 

4160 

4165 

4170 

4175 

Nl 

4030 

4035 

Wq 

4110 

4120 

4125 

4130 

Nth 

4105 

4120 

Wthe 

4120 

4125 

4130 

Nthv 

4115 

4125 

4130 

4170 

4175 

Zi 

4000 

4035 

404  5 

4030 

Zlcl 

4000 

4035 

4040 

4030 

USER  DEFINED 

FUNCTIONS: 

FNrfvi 

4035 

4125 

4130 

4135 

4140 

4145 

SUB  PROGRAMS : 

JUMP  TARGETS: 

4135 

4160 

SUB  PROGRAM  ParcelO 

COMMON  VARIABLES: 

VARIABLES: 

Dqdt 

4285 

4315 

Othe 

4285 

4290 

4295 

Gam 

4200 

4255 

4365 

4415 

Gamdew 

4205 

4345 

Lcp 

4210 

4240 

4250 

4235 

gio 

4195 

4215 

31 

4195 

4420 

Qq 

4215 

4385 

4435 

3t 

4215 

4240 

4380 

4385 

4420 

gv 

4195 

4275 

4235 

4315 

4385 

4420 

gvb 

4275 

4315 

4325 

4335 

4380 

3x 

44  3  5 

4440 

Rh 

4195 

4330 

4410 

Sig 

4330 

4335 

4430 

44  35 

T 

4260 

4235 

Td 

4230 

4415 

4440 

TdlO 

4230 

4345 

Tdo 

4255 

4260 

4305 

Th 

4195 

4250 

4255 

423  5 

4295 

4360 

4415 

The 

4195 

4240 

4400 

fhez 

4250 

4235 

4400 

4130 


4150 


4365 


122 


F-  10 


Thr 

4240 

4345 

4360 

Tx 

4260 

4305 

4325 

4365 

Z 

4195 

4235 

4255 

4355 

4365 

Zs 

4195 

4345 

4350 

4355 

Zz 

4220 

4235 

4265 

4330 

4370 

USER  DEFINED  FUNCTIONS: 

SUB  PROGRAMS: 

JUMP  TARGETS: 

4250  Theth: 

440  5 

4255  Cycle: 

4300 

4305  Dqdtcal: 

4230 

4325  Qvcal : 

4270 

4310 

4375 

4345 

4245 

4400 

4355 

4420 

4395 

4430  TdjEind: 

4225 

4390 

SUB  PROGRAM  Shrtwv 

COMMON  VARIABLES: 

VARIABLES : 

Cldthk 

4485 

4490 

4515 

4545 

Delz 

4490 

4515 

4570 

Dr  (*) 

4480 

44  95 

4500 

4505 

4510 

Dtemp 

4455 

4460 

4570 

Fbotdif 

4545 

Fbotdir 

4545 

Fheat 

4545 

4570 

Ftop 

4545 

I 

4530 

4535 

4540 

Qlbm 

4455 

4470 

4475 

Qq 

4475 

4515 

4520 

4525 

Q  t 

4455 

4515 

Rm 

4  525 

4535 

Sh 

4515 

4545 

Theta 

4455 

4460 

4465 

4545 

Xnd( *) 

4480 

4535 

4545 

Zb 

4455 

4460 

4465 

4435 

4570 

Zc 

4455 

4460 

4465 

4435 

USER  DEFINED  FUNCTIONS: 

SUB  PROGRAMS: 

Kahnrad  4  545 

JUMP  TARGETS: 


4370 
44  30 


4535 


123 


SUB  PROGRAM  Int779 


COMMON  VARIABLES: 
VARIABLES: 


lx 

4580 

4590 

4595 

4600 

4605  4615  4620 

4625 

4635 

4640 

4650 

Nxy 

4580 

4590 

4620 

Slope 

4580 

4635 

4640 

X 

4530 

4595 

4615 

4640 

Xarray( *) 

4  580 

4595 

4615 

4635 

4640 

Y 

4580 

4640 

Yarray( *) 

4580 

4635 

4640 

USER  DEFINED 

FUNCTIONS: 

SUB  PROGRAMS: 

JUMP  TARGETS: 

4595 

4610 

4630 

4615 

4595 

4635 

4615 

4655 

4645 

5020 

4600 

4620 

SUB  PROGRAM  A52 

COMMON  VARIABLES: 

VARIABLES: 

U 

4665 

4690 

4720 

4695 

4725 

4700 

4730 

4705 

4735 

4710 

4715 

X  kv ( * ) 

4665 

4685 

4715 

4690 

4720 

4695 

4725 

4700 

4730 

4705 

4735 

4710 

USER  DEFINED  FUNCTIONS: 
SUB  PROGRAMS: 

JUMP  TARGETS: 


SUB  PROGRAM  A 20 
COMMON  VARIABLES: 

VARIABLES: 

A(*)  4750  4830  4835  4840  4345  4350  4855 


124 


F-  12 


At  ( *) 

4750 

4765 

Ag 

4750 

4845 

4850 

4855 

A1 

4750 

48  20 

4840 

4855 

Be 

4750 

4825 

4840 

4855 

D 

4750 

4765 

01 

4795 

4805 

4810 

4820 

4325 

02 

4300 

4805 

4810 

03 

4815 

4820 

4825 

Dr  (*) 

4750 

4765 

Fi  (*) 

4750 

4770 

Fo 

4770 

43  20 

4825 

4855 

Fq 

4775 

4785 

4790 

G 

4765 

4775 

4730 

4800 

43  20 

4825 

I  w 

4750 

4765 

4  770 

Om 

4765 

4735 

4790 

4795 

4300 

4320 

p 

4750 

4810 

4830 

4835 

4845 

4850 

Qext(  *) 

4750 

4765 

Sa(  *) 

4750 

4765 

Sh 

4750 

4765 

T 

4750 

4765 

4785 

4345 

4350 

43  55 

Tau  (* ) 

4750 

4765 

Xk 

4750 

4805 

4815 

4845 

4350 

Xkv(* ) 

4750 

4765 

Xmu 

4750 

4815 

48  20 

4325 

48  5  5 

Xnd(  *) 

4750 

4765 

USER  DEFINED 

FUNCTIONS: 

SUB  PROGRAMS 

• 

• 

A28 

4765 

JUMP  TARGETS 

• 

• 

SUB  PROGRAM  A30 

COMMON  VARIABLES: 

VARIABLES: 

A1 

4870 

4890 

4915 

Be 

48  70 

4895 

4920 

Fbotdif 

43  70 

4935 

Fbotdir 

4870 

4940 

Fd 

4900 

4905 

4910 

4930 

4940 

F 1  (  * ) 

48  70 

4900 

FI 

48  70 

4925 

4930 

Ftop 

43  70 

4905 

Ftopd 

43  70 

4910 

10 

4890 

4905 

4915 

4935 

11 

4895 

'4905 

4920 

4925 

4935 

Iw 

48  70 

4900 

4825 


125 


F-  13 


p 

4870 

4895 

4920 

r 

4870 

4915 

4920 

4930 

4940 

TO 

4885 

48  90 

4895 

X(* 

)  4870 

4890 

4895 

4915 

4920 

Xk 

48  70 

4890 

4895 

4915 

4920 

Xmu 

4870 

43  90 
4940 

4895 

4900 

4915  4920  4930 

USER 

DEFINED  FUNCTIONS: 

SUB 

PROGRAMS : 

JUMP 

TARGETS : 

SUB  PROGRAM  A28 
COMMON  VARIABLES: 


VARIABLES: 


Af  (  * ) 

4955 

5010 

D 

4955 

4980 

Dr  (*) 

4955 

4980 

VJ 

4955 

4970 

Iw 

4955 

4980 

J 

4975 

4980 

5015 

Ora 

4955 

4970 

Qext(*) 

4955 

4980 

Sa  (*) 

4955 

5005 

Sh 

4955 

5020 

T 

4955 

4995 

Tau(*) 

4955 

4980 

Ts 

4970 

4985 

Vp 

5020 

5025 

Xkv(  *) 

4955 

5025 

Xnd  (*) 

4955 

4980 

5020 

5010 

5005 

5010 

5025 

4985 

4990 

5000  5005  5010 

5005 

5025 

5025 

4985 

5005 

5010 

4995 

5005 

5010 

USER  DEFINED  FUNCTIONS: 
SUB  PROGRAMS: 

JUMP  TARGETS: 


SUB  PROGRAM  A 100 
COMMON  VARIABLES: 


VARIABLES: 
Af  (*) 


5150 


F-  14 


«SEKB wate 


5040 


126 


Or(*> 

5040 

5090 

5125 

I 

5120 

5125 

5135 

1x9 

5105 

5110 

5160 

5115 

Cext( *) 

5040 

5135 

5140 

Qsca 

5135 

5140 

R0 

5130 

5150 

R2 

5125 

5130 

3a  (*) 

5040 

5140 

5145 

Xk9  (*) 

5055 

5095 

X10 

5110 

5130 

X19  (* ) 

5055 

5090 

5110 

XnO 

5115 

5130 

Xn9  ( *) 

5055 

5090 

5115 

USER  DEFINED  FUNCTIONS: 


SUB  PROGRAMS: 

JUMP  TARGETS: 

5060  5035 


SUB  PROGRAM  Kahnrad 
COMMON  VARIABLES: 
VARIABLES: 


A  (* ) 

5190 

5330 

5350 

5410 

5425 

Af  (  * ) 

5195 

5310 

5330 

Ag 

5285 

5290 

5330 

A1 

5330 

5  500 

Alam( *) 

5200 

5235 

5275 

Be 

5330 

5  500 

Cldthk 

5175 

5240 

5295 

D 

5  30  5 

5330 

5515 

Dd 

5350 

5355 

5395 

Delz 

5295 

5305 

Dr  (*) 

5200 

5310 

5330 

Fbotdif 

5175 

5245 

5500 

Fbotdir 

5175 

5245 

5500 

Fheat 

5175 

5515 

Fi(*> 

5205 

5235 

5275 

FI 

5300 

5  500 

5515 

Ftop 

5245 

5500 

5515 

Ftopd 

5245 

5500 

5520 

Ftopdir 

5245 

I 

5270 

5275 

5280 

11 

5345 

5360 

5365 

Ibi 

5405 

5410 

5415 

127 


5140  5145  5150  5155 

5135  5140  5145  5150 


5375  5380  5385  5395  j 

5435  5465  5470  5490  ] 


5410  5425  5435 


5535 

5530 

5330  5500 

5525 

5485  5490  5495 

5380  5385 


F-  15 


S^3SEBfe 


IincO 

5175 

5520 

5525 

I  test 

5175 

5250 

Iw 

5325 

5330 

5500 

5505 

Iwi 

5370 

5375 

5330 

5335 

5390 

J 

5420 

5425 

5475 

5435 

5445 

5460 

5465 

5470 

K 

5430 

5435 

5430 

5440 

5455 

5460 

5465 

5470 

L 

5340 

5345 

5350 

5375 

5380 

5395 

5410 

5420 

5425 

5435 

5450 

P 

5330 

5500 

Cext( *) 

5195 

5310 

5330 

R(*) 

5190 

5375 

5385 

Sa(*) 

5195 

5310 

5330 

Secmu 

5265 

5275 

Sh 

5175 

5240 

5330 

T 

5330 

5500 

Ta  u ( * ) 

5200 

5330 

Theta 

5175 

5250 

5255 

U 

5240 

5315 

X{*) 

5190 

5490 

5500 

Xk 

5330 

5500 

Xkv(  * ) 

5  20  5 

5315 

5330 

Xmu 

5255 

5260 

5265 

5290 

5330 

5500 

X  nd  (  * ) 

5175 

5330 

USER  DEFINED  FUNCTIONS: 


SUB  PROGRAMS: 


A100 

5310 

A20 

5330 

A30 

5500 

A52 

5315 

JUMP  TARGETS: 

5215 

5230 

5285 

5260 

5355 

5400 

5405 

5355 

5545 

5365 

SUB  PROGRAM  Sky2 
COMMON  VARIABLES: 
VARIABLES: 


Emit( *) 

5580 

5700 

5725 

5735 

5760 

£s 

5725 

5750 

Ftop 

5555 

5775 

5790 

5795 

Gad 

5605 

5695 

5730 

5750 

5735 

123 


F-  15 


I 

5705 

5715 

I  cloud 

5655 

5715 

Ie 

5725 

5760 

Ii 

5605 

5700 

I  top 

5670 

5795 

5675 

Nemis 

5610 

5725 

Nsnd 

5555 

5655 

Pot_temp(*) 

5555 

5755 

5595 

Qtsnd{ *  j 

5580 

5590 

RhoO 

5605 

57  20 

S 

5605 

5735 

Sig 

5605 

5735 

Spec_hum(*) 

5555 

5590 

Temis( *) 

5580 

5600 

Thrsnd( *) 

5530 

5595 

Tsky 

5555 

5795 

Tt 

5730 

5735 

Tts 

5695 

5735 

Tu(*) 

5580 

5600 

U 

5605 

57  20 

Zi 

5555 

5675 

Zscale 

5605 

5720 

Zsnd (*) 

5555 

5675 

USER  DEFINED 

FUNCTIONS: 

5785 

SUB  PROGRAMS: 

Int779 

5725 

5760 

JUMP  TARGETS: 

5560 

5585 

5685 

5675 

5785 

5715 

5795 

5780 

SUB  PROGRAM  Zenith 
COMMON  VARIABLES: 


VARIABLES: 


Decl 

5835 

5840 

a 

5830 

5835 

Hr 

5820 

5340 

Julday 

5810 

5330 

Lat 

5310 

5325 

Latr 

5825 

5840 

Theta 

5810 

5340 

5720 

5730 

5745 

5770 

5785 

5710 

5725 

5735 

5760 

5765 

5680 

5760 

5685 

5695 

5705 

5660 

5770 

5665 

5670 

5705 

5750 

5720 

5755 

5755 

5765 

5795 

5765 

5790 

5795 

5725 

5760 

5695 

5730 

5750 

5785 

5740 

5750 

5765 

5735 

5790 

5740 

5725 

5760 

5725 

5755 

5760 

5755 

5695 

5720 

5730 

57  50 

5755 

5845 

F-  17 


.  (9  ’ 


Time 


5310  5320 


USER  DEFINED  FUNCTIONS: 
SUB  PROGRAMS: 

JUMP  TARGETS: 
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